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solution.  The  existence  theory,  which  is  in  the  setting  of  semigroups  and  < 
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solution  are  new. 
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depends  on  the  solution.  This  introduces  a  nonlinearity  to  the  problem  which  makes  it  of  the  form  frequently 
called  quasilinear.  This  is  a  stronger  form  of  nonlinearity  than  semilinear.  Results  of  a  numerical  example  are 
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Preface 


My  research  proposal  was  based  on  my  search  for  an  unsolved  problem  concerning 
the  dynamics  of  large  space  structures.  My  strategy  was  to  identify  some  feature  of  the 
dynamics  for  which  adequate  mathematical  tools  were  not  available.  In  my  reading,  on 
several  occasions,  I  encountered  phrases  about  nonlinearities  not  being  adequately  mod¬ 
elled.  I  determined  that  some  of  the  simple  traditional  models  for  beam  behavior  did 
not  have  nonlinear  counterparts.  That  is  to  say,  while  many  sophisticated  models  had 
been  developed,  some  rather  straightforward  models  apparently  had  not.  I  was  curious 
whether  the  mathematics  was  available  for  such  models.  Finding  none,  and  obtaining  the 
concurrence  of  my  committee,  I  started  my  research. 

I  do  not  know  how  well  my  model  will  serve  for  analysis  of  space  structures.  I  do 
know  that  the  analysis  of  mathematical  tools  has  been  fascinating  and  enjoyable. 

I  would  like  to  extend  many  thanks  to  the  courteous  and  professional  staff  of  the 
AFIT  Library.  They  have  efficiently  handled  countless  Intejrlibrary  Loan  requests  and 
have  helped  in  a  variety  of  local  searches.  I  thank  them  all. 

I  would  also  like  to  mention  that  when  my  proposed  research  topic  was  still  too 
broad,  Dr.  Bagley  posed  several  questions  which  were  very  effective  in  helping  me  to  focus 
more  rapidly  on  a  specific  equation.  I  thank  him. 

Many  thanks  are  also  due  to  Dr.  Lair,  who  answered  a  variety  of  technical  ques¬ 
tions  and  pointed  out  several  refinements.  His  countless  hours,  in  my  behalf,  are  greatly 
appreciated. 

Finally,  I  thank  Dr.  Quinn.  He  has  been  a  constant  source  of  encouragement  and 
guidance.  He  was  always  eager  and  helpful  in  guiding  me  in  the  pursuit  of  my  research. 
He  went  to  great  effort  to  guide  my  efforts  as  a  researcher  so  that  they  would  be  most 
effective.  Yet,  he  was  careful  not  to  redirect  me  along  any  lines  other  than  those  I  wished 
to  pursue.  He  taught  me  to  do  research. 


Carl  Edwin  Crockett 
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Abstract 

A  nonlinear  partial  differential  equation,  motivated  by  the  transverse  vibration  of  a 
beam,  is  shown  to  have  a  unique  solution.  The  existence  theory,  which  is  in  the  setting  of 
semigroups  and  evolution  operators,  is  a  composite  and  synthesis  of  theorems  of  Kato.  The 
formulation  of  the  problem  and  the  verification  that  the  formulation  leads  to  a  solution 
are  new. 

The  introductory  chapter  provides  background  on  the  topic  generally.  Chapter 
provides  detailed  formulations  for  the  constant  coefficient  case  Chapter  3  describes  nonau- 
tonomouc  cases.  The  most  general  theorem  is  presented  hero.  In  Chapter  4,  a  move  general 
case  is  considered.  Namely,  Kelvin- Voigt  damping  with  a  coefficient  which  depends  on  the 
solution.  This  introduces  a  nonlinearity  to  the  problem  which  makes  it  of  the  form  fre¬ 
quently  called  quasilinear.  This  is  a  stronger  form  of  nonlinearity  than  semilinear.  Results 
of  a  numerical  example  are  presented. 
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AN  EVOLUTION  OPERATOR  SOLUTION  FOR 
A  NONLINEAR  BEAM  EQUATION 


I.  Introduction 

The  central  problem  of  this  dissertation  concerns  a  nonlinear  partial  differential  equa¬ 
tion  for  the  modelling  of  the  transverse  vibration  of  a  beam.  The  particular  form  of  non¬ 
linearity  that  arises  gives  the  equation  a  form  that  is  frequently  referred  to  as  quasilinear. 
This  is  a  stronger  nonlinearity  than  the  form  known  as  semilinear.  The  equation  is  shown 
to  have  a  unique  solution.  The  applicable  existence  theory  is  a  composition  of  some  the¬ 
orems  of  Kato.  The  formulation  of  the  problem  and  the  verification  that  the  formulation 
leads  to  a  solution  are  new. 

This  introductory  chapter  provides  background  on  the  topic  generally.  Chapter  2 
provides  detailed  formulations  for  the  constant  coefficient  cases.  Chapters  3  and  4  provide 
the  generalizations. 

1.1  Some  general  background 

Partial  differential  equations  are  an  essential  element  in  the  mathematical  modelling 
of  the  physical  world.  The  topics  that  arise  in  the  study  of  natural  phenomena  with 
partial  differential  equation  models  are  numerous  and  varied.  Indeed,  the  fact  that  two 
physical  phenomena  are  very  similar  does  not  guarantee  any  particular  connection  between 
the  corresponding  mathematical  models.  In  fact,  two  models  of  the  same  behavior,  with 
one  slightly  more  sophisiticated  than  the  other,  can  lead  to  entirely  distinct  mathematical 
entities.  Even  changes  which  appear  to  be  very  minor  or  superficial  can  radically  alter  the 
nature  of  the  solution,  and  even  whether  there  is  a  solution.  In  cases  where  a  solution  still 
exists,  it  may  be  necessary  to  use  entirely  different  methods  to  find  it  (or  them). 

In  the  classical  theory  of  partial  differential  equations  the  three  basic  second  order 
forms  are 

1.  the  heat  equation  ut  =  uxx 
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2.  the  wave  equation  Utt  =  uxx 

3.  Laplace’s  equation  ««  =  —uxx. 

It  is  well  known  that  these  have  quite  distinct  solution  characteristics.  The  same  holds 
true  for  higher  order  equations.  When  problems  are  formulated  abstractly,  it  remains  true 
that  small  changes  in  appearance  can  lead  to  significant  changes  in  solution  properties. 

The  equations  above  have  been  given  names  to  distinguish  their  types.  The  heat 
equation  is  known  as  an  example  of  a  parabolic  equation.  It  has  a  unique  real  characteristic 
curve.  The  wave  equation  is  called  hyperbolic.  It  has  two  distinct  characteristic  curves. 
The  Laplace  equation  is  elliptic  and  has  no  real  characteristic  curves.  The  classifications 
are  used  to  distinguish  equations  whose  solutions  have  certain  distinctive  behavior. 

As  an  example  of  a  model  refinement  which  may  not  look  too  different,  Friedman 
and  Hu  [17:pp  249,252]  discuss  an  equation  referred  to  as  the  hyperbolic  heat  equation.  It 
arises  when  a  slight  modification  to  a  basic  mathematical  assumption  avoids  the  physically 
unrealistic  property,  of  the  usual  model,  that  global  temperature  changes  instantaneously. 
The  change  is  to  replace  the  Fourier  Law  q  =  -k6x  with  q  +  rqt  =  -kOx  for  some  time 
delay  r,  probably  quite  small.  The  resulting  equation  is 

T$tt  +  &t  ~  Oxx  =  0.  (1) 

In  the  case  of  abstract  equations  the  affect  of  a  negative  sign  or  of  a  change  in  the 
order  of  the  equation  can  be  pointed  out  by  observing  that  the  abstract  equations 

Ut  =  ^xx 

Utt  ~  uxx 

Utt  ~  ~Uxxxx  (2) 

all  have  classical  solutions.  But,  in  their  abstract  formulation  with  homogeneous  Dirichlet 
boundary  conditions,  none  of  the  equations 


Ut  —  uxx 


2 


utt 


(3) 


—  uxx 
Utt  —  UXxxx 


has  a  nontrivial  classical  solution. 

The  terms  parabolic  and  hyperbolic  do  not  have  the  same  meaning  in  the  abstract 
theory  as  they  do  for  classical  second  order  equations.  (The  term  elliptic  is  not  used  for 
abstract  equations.)  An  abstract  equation  of  the  form 


ut  =  Au  (4) 

is  called  parabolic  if  the  operator  A  is  the  generator  of  an  analytic  semigroup.  It  is  called 
hyperbolic  if  A  is  the  generator  of  only  a  Cq  semigroup  (see  Definition  2). 

The  present  work  addresses  an  Euler-Bernoulli  beam  equation  with  Kelvin- Voigt 
damping.  The  abstract  formulation  results  in  an  equation  which  is  hyperbolic.  Linear  and 
nonlinear  versions  will  be  considered  as  well  as  autonomous  and  nonautonomous  versions. 


1.2  An  introduction  ti  the  basic  equation 

A  fundamental  element  of  most  large  structures  is  the  beam.  Over  a  period  of 
many  years,  as  the  need  for  more  and  more  accurate  modelling  has  become  apparent, 
the  mathematical  equations  used  to  describe  beam  behavior  have  become  more  and  more 
sophisticated.  Similarly,  more  advanced  mathematical  tools  are  used  in  the  development 
and  analysis  of  such  equations.  Herein,  attention  will  be  focused  on  equations  for  the 
transverse  vibration  of  a  beam. 


In  particular,  an  equation  with  nonlinear  damping  will  be  considered.  This  will 
be  preceeded  by  the  treatment  of  simpler  equations  in  order  to  build  up  the  appropriate 
mathematical  machinery.  Existence  and  uniqueness  of  a  solution,  in  the  simple  cases, 
will  be  established  from  the  viewpoint  of  semigroup  theory.  The  more  complex  cases  will 
require  the  mild  generalization  of  a  semigroup  known  as  an  evolution  system. 

In  this  section  a  review  of  an  early  model  of  beam  motion  will  be  presented.  In 
the  next  chapter  the  model  will  be  considered  from  a  mathematically  sophisticated  ( ie 


3 


Figure  1.  Cantilevered  beam 

semigroups)  point  of  view.  This  will  provide  a  point  of  departure  for  the  investigation 
of  more  complex  and  potentially  more  precise  mathematical  models  of  beam  behavior. 
Development  of  the  mathematics  related  to  these  models  is  the  pm  pose  of  the  present 
work.  It  is  stated  at  the  outset  that  while  beam  equations  motivated  this  mathematical 
problem,  the  engineering  is  left  to  others.  That  is  to  say,  while  the  equations  developed 
may  be  of  engineering  value,  it  will  not  be  proved  here.  The  current  work  is  limited  to 
establishing  of  the  validity  of  certain  mathematical  operations. 

1.2.1  The  traditional  starting  place 

Consider  the  cantilevered  beam.  This  is  a  beam  which  is  fixed  at  one  end,  as  to  a 
wall,  and  free  at  the  other  end  (see  Figure  1).  Assume  a  thin  beam  whose  density  and 
other  physical  characteristics  are  uniform.  It  is  customary  to  establish  a  coordinate  system, 
assumed  inertial,  whose  origin  is  at  the  intersection  of  the  beam’s  centerline  and  the  wall. 
The  a;- axis  is  chosen  along  the  beam’s  centerline.  It  is  supposed  that  the  centerline  is 
straight  and  represents  the  natural,  rest,  or  unperturbed  position  of  the  beam.  Let  L  be 
the  length  of  the  beam. 

Vertical  motion  of  each  point  x  along  the  beam,  with  vertical  displacement  at  time 
t  represented  by  u(t,x),  is  considered.  Notice  that  the  physical  configuration  requires  (for 
consistency)  the  following  mathematical  assumptions: 

u(t,  0)  =  0 

vx(t,  0)  =  0 

^ara;(t,  jL)  —  0 
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(5) 

(6) 
(7) 


(8) 


(Subscripts  indicate  partial  differentiation  with  respect  to  the  variable  appearing  as  the 
subscript.)  The  physical  interpretations  of  the  boundary  conditions  are  as  follows:  (5)  says 
that  the  left  end  of  the  beam  is  held  in  a  fixed  position,  (6)  says  that  the  beam  centerline  is 
perpendicular  to  the  wall,  (7)  says  that  there  is  no  curvature  at  the  right  end  of  the  beam, 
and  (8)  says  that  the  curvature  at  the  right  end  is  not  changing  (ie  there  is  no  external 
torque). 

A  large  variety  of  physical  situations  can  be  described  with  only  slight  changes  in 
(5)-(8).  For  example  if  both  ends  were  free,  conditions  like  (7)  and  (8)  would  apply  at  the 
left  end.  If  a  control  mechanism  were  attached  at  the  right  end,  then  conditions  (7)  and 
(8)  would  have  nonzero  right  hand  sides.  Clearly  there  are  many  other  physical  situations 
described  by  closely  related  sets  of  conditions.  Of  course,  u(t,  x)  must  be  sufficiently 
differentiable  for  the  boundary  conditions  to  make  sense.  An  assumption  to  this  affect  is 
consistent  with  typical  physical  situations. 

The  simplest  model  for  the  transverse  motion  of  a  beam  is  called  the  Euler-Bernoulli 
beam  model.  It  has  been  in  use  for  over  one  hundred  years.  (See  Russell  [67:pp  177-216] 
for  an  excellent  historical  review  of  models.)  The  following  assumptions  are  made: 


1.  The  density  of  the  beam,  m(x),  is  known  and  appropriately  smooth.  It  is  assumed 
to  be  constant  with  respect  to  t. 

2.  Vertical  displacements  are  small  compared  to  the  length  of  the  beam  and  any  hori¬ 
zontal  component  of  motion  is  neglected. 


3.  Young’s  modulus,  E,  (an  experimentally  determined,  material  dependent  quantity 
that  relates  the  amount  of  stretching  to  the  amount  of  force  applied,  see  Kolsky 

faV’ncr  Ql  Cryr  ovamnloi  prop  mnmont  ryf  inortia  nlnno-  Kon^incr  pyic  T  pro 

known  and  sufficiently  smooth  (often  taken  as  constant).  They  are  assumed  to  be 
constant  with  respect  to  time. 


4.  Mechanical  energy  is  conserved. 
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Details  of  the  derivation  are  provided  in  Appendix  A.  The  resulting  equation  is 


utt  +  (auxx)xx  =  0;  0  <x  <  L;  0  <t 


(9) 


where  a  = 

This  is  the  simplest  mathematical  model  for  a  beam.  Transverse  vibration  is  allowed 
and  elasticitv  of  bending  is  considered.  Longitudinal  stretching  and  vibration  are  ignored, 
as  are  time  ndent  changes  in  material  properties.  Also  ignored  are  changes  in  material 
properties  due  to  the  history  and  current  state  of  the  beams  position  and  motion.  Damping 
is  not  modelled  and  other  imperfections  are  sure  to  be  present.  Nevertheless,  some  starting 
place  is  needed  and  the  Euler-Bernoulli  model  is  the  usual  place.  It  provides  a  specific 
equation  for  use  in  demonstrating  the  mathematical  tools  to  be  used  in  further  analysis. 

1.2.2  Proposed  generalizations  and  analysis  of  the  equation 

The  linear  partial  differential  equation  above  is  reduced  to  a  system  which  is  first 
order  in  time.  Then  an  appropriate  Banach  space  is  chosen  and  the  system  is  formulated 
as  an  abstract  Cauchy  problem.  Existence  and  uniqueness  of  a  solution  is  then  established 
from  standard  theorems.  A  key  challenge  is  to  construct  the  Banach  space  wisely  so  that 
the  operator  in  the  abstract  problem  is  densely  defined  and  has  the  properties  necessary 
to  justify  the  application  of  the  theorems. 

The  quantity  a  will  then  be  allowed  to  vary,  thus  generalizing  the  problem.  First,  a 
will  be  allowed  to  vary  with  the  space  variable.  This  case  does  not  require  any  additional 
theory. 

The  next  generalization  will  be  the  introduction  of  a  Kelvin- Voigt  damping  term.  In 
this  case  the  equation  has  the  form 

utt  -b  ( fi'U'txx)xx  -b  QV'xxxx  =  9.  (19) 

The  underlying  assumption  of  Kelvin- Voigt  damping  is  that  the  damping  depends  on 
velocity  in  the  same  way  the  basic  equations  depends  on  position.  A  positive  value  of 
/?  represents  damping  in  the  system.  Negative  values  of  (3,  which  represent  energy  input 
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or  energy  generation,  will  not  be  considered.  To  maintain  the  step  by  step  approach 
to  generalization,  a  is  held  constant  while  the  damping  term  is  considered.  Separate 
treatments  will  be  given  to  the  cases  of 

1.  constant  coefficient  of  damping, 

2.  spatially  varying  coefficient  of  damping, 

3.  temporally  varying  coefficient  of  damping, 

4.  combined  space  and  time  variation  of  the  coefficient  of  damping,  and 

5.  solution  dependent  coefficient  of  damping. 

It  is  in  the  last  of  these  that  a  true  nonlinearity  appears.  The  form  of  the  equation 
with  this  nonlinearity  is  such  that  the  term  quasi-linear  is  appropriate.  Interest  in  this 
particular  equation  motivated  the  current  research. 

In  another  generalization,  a  will  be  allowed  to  vary  with  time.  This  case  requires 
a  more  powerful  existence  theorem.  A  more  general  theorem,  which  is  a  synthesis  of  two 
theorems  of  Kato  and  covers  this  case,  is  presented. 


1.3  A  review  of  the  literature 

The  review  of  the  literature  can  be  divided  into  three  basic  areas.  First,  a  general 
overview  of  some  problems  which  use  the  same  general  type  of  theory  as  that  which  will 
be  used  in  this  work  is  given.  Second,  some  problems  which  are  closely  related  (at  least 
in  general  appearance)  to  the  current  problem  are  presented.  Finally,  a  review  of  the 
literature  which  provides  the  theoretical  framework  for  this  work  is  presented. 


1.3.1  A  sampling  of  problems  from  the  literature 

Higher  order  equations  are  used  to  describe  beam  behavior.  Ball  [3:pg  399]  describes 
a  model  which  accounts  for  several  affects.  The  equation  is 


Utt  +  OiUxxxx  -  \P  +  k 


/ 


u\dx 


Uxx  +  '  'U'txxxx 


-a  f  uxutx  dx  uxx  +  6ut 
Jo 


0 


(11) 
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He  discusses  existence,  uniqueness,  and  regularity  for  the  constant  coefficient  problem. 
The  approach  is  more  classical  than  modern  (te  not  a  semigroup  approach). 

Pivovarchik  [59:pg  647]  considers  the  spectrum  of  a  similar  looking  equation,  namely 
Ot'U'tXXXX  "1"  'U'XXXX  +  ( g(x)ux)x  +  k(x)ut  +  utt  =  o.  (12) 

More  specific  examples  from  the  literature  will  be  presented  shortly.  A  brief  pause 
to  make  some  general  comments  is  appropriate. 

The  current  literature  contains  a  wide  variety  of  problems  related  to  fluid  flow,  plate 
vibration,  beam  vibration,  and  the  control  of  these  kinds  of  phenomena.  Since  small 
alterations  often  make  problems  which  are  mathematically  distinct,  this  literature  is  quite 
voluminous.  No  complete  review  will  be  attempted.  However,  it  is  certainly  appropriate 
to  present  a  sample  of  the  kind  of  work  that  is  being  done. 

A  brief  list  of  the  kinds  of  items  involved  in  the  set  up  of  a  problem  is  now  given. 
A  small  change  in  any  of  these  items  can  lead  to  entirely  different  solution  behavior.  The 
domain  may  be  bounded  or  unbounded,  fixed  or  changing,  and  the  conditions  specified  at 
the  boundary  can  have  a  variety  of  forms.  The  operators  which  appear  may  be  bounded  or 
unbounded,  linear  or  nonlinear,  autonomous  or  nonautonomous,  they  may  be  self-adjoint, 
skew-adjoint,  compact,  or  closed.  The  underlying  space  may  be  Hilbert,  Banach,  normed, 
linear,  metric,  or  just  a  set.  The  order  of  the  equation  is  significant. 

For  all  of  these  factors,  small  changes  can  affect  whether  the  problem  is  well-posed, 
whether  it  has  a  solution,  whether  solutions  are  unique,  what  kind  of  algorithm  can  be 
used  to  obtain  the  solution,  or  indeed,  what  notion  of  solution  is  appropriate  to  consider. 

A  few  samples  of  equations  in  the  recent  literature  are  given. 

Bernis  [5:pg  227]  establishes  existence  and  uniqueness  for  the  equation 

(~&)mu  + g(x,u)  =  f  (13) 


8 


on  all  of  %tN  with  certain  limitations  on  /  and  g.  (The  symbol  A  is  used  for  the  Laplace 
operator.)  He  also  establishes  existence  and  uniqueness  for  the  parabolic  problem 

ut  +  (-A)mu+  |  u  |p_1  u  =  /  (14) 

on  3lN  x  (0, 1).  Again,  certain  restrictions  apply. 

The  Schrodinger  equation,  [68:pg  823], 

Q 

— u  +  iAu  +  |u|p_1u  =  0  (15) 

where  i  —  (— 1)?,  also  receives  attention. 

Cannarsa  et.  al.  [7:pg  2]  treat  a  damped  wave  equation 

utt(i,x)=:  A(u(t,x)  + Kut(t,x))  + f(t,x)  (16) 

on  a  bounded  domain  with  appropriate  boundary  and  initial  conditions.  They  note  (pg 
3)  that  even  a  slight  change  in  boundary  conditions  sends  the  problem  into  the  realm  of 
the  unknown.  Weak  solutions  are  obtained  following  a  transformation  which  changes  the 
problem  to  one  where  the  operator  is  the  generator  of  an  analytic  semigroup.  They  claim 
this  approach  can  be  extended  to  higher  order  (eg  beam)  equations. 

Kuttler  and  Hicks  [42:pg  1]  address  existence  and  uniqueness  of  a  global  weak  solution 
to 

mtt  +  (P(m«))*  -  (a(mx)mxt)x  =  f(t,  x ).  (17) 

Their  emphasis  is  on  time  dependent  boundary  conditions. 

Other  uses  of  similar  theory  are  found  in  Euler  Equations  [2:pp  367-382],  porous 
medium  systems  [21:pg  86],  and  Navier-Stokes  equations  [36:pg  891],  But,  it  is  time  to 
turn  to  equations  more  closely  related  to  the  current  work. 

1.3.2  Review  of  literature  related  to  the  current  problem 
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Fitzgibbon  [15:pg  536]  addresses  (11)  as  a  specific  example  that  comes  under  an  ab¬ 
stract  formulation  that  he  gives  for  a  class  of  quasilinear  evolution  equations.  He  examines 

utt  +  cxAut  +  Au  =  F(t,  u,  ut)  (18) 

with  appropriate  initial  conditions.  The  operator  A  is  allowed  to  be  unbounded,  but  it 
does  not  depend  on  t  or  u.  Also,  a  is  a  constant.  When  A  is  chosen  to  correspond  to  the 
problem  of  Ball,  it  turns  out  that  A  is  self-adjoint  and  —A  is  the  generator  of  an  analytic 
semigroup,  see  also  [71:pp  631-633].  Existence  and  uniqueness  results  are  obtained. 

Huang  [23:pg  714]  discusses  the  closely  related  equation 

utt  +  But  +  Au  =  0  (19) 

where  B  is  related  to  A  in  a  certain  way. 

Authors  concerned  with  control  theory  still  tend  to  use  the  more  basic  forms  of  the 
equation  while  they  concentrate  on  progress  in  the  area  of  controls.  For  example,  Lasiecka 
and  Triggiani  [44:pg  330]  and  [45:pg  1]  use  the  relatively  mple  form 


Utt  +  A  2u  =’  0.  (20) 

(See  also  [39:pg  288],  for  a  classical  treatment) 

Some  recent  uses  of  Kelvin- Voigt  damping  appear  in  [4:pg  l]  and  [6:pg  1391].  The 
first  uses  it  in  solving  for  a  material  property  as  a  distributed  parameter.  The  second 
solves  for  a  displacement  and  stress  field  in  a  solid. 

Standard  treatments  of  semigroup  theory  are  in  [22],  [19],  [20],  and  [58], 

Several  authors  have  addressed  numerical  issues  related  to  the  implementation  of 
iterative  schemes.  An  early  paper  discussing  the  convergence  of  difference  schemes  is  [55:pg 
321].  For  a  semigroup  style  numerical  analysis  book,  see  as  an  example,  [64].  Another 
reference  on  discrete  schemes,  with  emphasis  on  time  dependent  operators,  is  [57].  Another, 
with  emphasis  on  nonlinearity,  is  [26]. 
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1.S.3  Some  methods  of  analysis 

Two  lines  of  argument,  for  analysis  of  the  types  of  equations  that  have  been  men¬ 
tioned,  will  now  be  briefly  discussed.  Many  authors  have  written  on  these  issues,  but  it  is 
not  necessary  to  give  a  comprehensive  review.  The  two  lines  of  argument  to  be  mentioned 
are  due  principally  to  Kato  and  Crandall.  A  review  of  the  bibliographies  of  the  articles 
cited  will  quickly  lead  to  broad  coverage  of  the  topic. 

The  development  of  the  theory  of  semigroups  goes  back  several  decades.  It  is  appro¬ 
priate  to  pick  up  the  story  with  an  article  by  Kato  in  1953  [28].  The  equation  of  interest, 
in  his  paper,  is 

^u(t)  =  A(*)u(t)  +  /(i).  (21) 

This  is  called  an  evolution  equation  because  of  the  time  dependence  in  the  definition  of 
the  operator  A.  He  identifies  conditions  for  the  existence  of  a  unique  solution. 

In  the  years  following  this  publication,  several  generalizations  were  obtained.  Most 
of  the  extensions  were  in  the  direction  of  weakening  one  hypothesis  or  another.  Of  course, 
the  goal  of  such  work  was  to  widen  the  range  of  applicability  of  the  theory. 

The  method  of  proof,  of  the  existence  of  a  solution,  involves  construction  of  a  sequence 
of  operators  which  are  shown  to  converge  to  the  operator  in  the  original  problem.  The 
modified  problems  corresponding  to  each  of  the  new  operators  are  solvable.  Results  which 
confirm  that  the  function,  to  which  the  sequence  of  solutions  converges,  is  a  solution  to  the 
original  problem  are  of  interest.  In  this  regard  the  paper  of  Trotter  [70]  and  a  correction 
to  the  proof  of  one  of  its  theorems  by  Kato  [29]  are  applicable. 

By  the  mid  60’s,  nonlinear  versions  were  receiving  appreciable  attention.  Of  course, 
in  this  more  complex  setting,  many  distinct  variations  of  a  problem  can  cause  it  to  fail  to  be 
linear.  Hence,  the  term  nonlinear  is  used  with  several  different  meanings  in  the  literature. 
In  some  cases  it  simply  means  that  a,  function  is  set- valued  rather  than  single- valued.  If  the 
forcing  function  depends  on  the  solution  the  term  semilinear  is  used.  When  the  operator 
depends  on  the  solution,  the  term  quasilinear  is  used.  The  terms  genuinely  nonlinear 
and  fully  nonlinear  are  also  encountered.  Any  of  the  above  terms  may  be  abbreviated  to 
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nonlinear  in  the  literature.  Hence,  it  is  important  to  be  wary  of  one’s  own  preconceived 
notions  of  what  the  word  suggests. 

The  1971  paper  of  Crandall  and  Liggett  [8]  is  considered  definitive  in  bringing  the 
nonlinear  problem  well  in  hand.  It  establishes  convergence  of  an  exponential  type  limit 
which  converges  to  the  solution,  if  one  exists. 

Also  in  the  early  70’s,  Kato  was  publishing  in  the  area  of  linear  evolution  equations 
of  hyperbolic  type  [31]  [32].  Recall  that  if  the  operator  of  an  equation  generates  a  strongly 
continuous  semigroup,  then  the  equation  is  called  hyperbolic.  An  equation  whose  operator 
is  the  generator  of  an  analytic  semigroup  is  called  parabolic. 

By  the  mid  70’s,  Kato  was  working  on  quasilinear  evolution  equations  of  hyperbolic 
type  [34]  [35],  In  these  papers,  conditions  are  given  for  existence  and  uniqueness  of  solutions 
and  several  applications  are  discussed.  Specific  applications  are  discussed  in  [24]  and  [33]. 

Through  the  80’s,  Crandall  and  Souganidis  have  published  [9],  [10],  and  [11]  in  the 
area  of  nonlinear  equations.  Their  approach  is  to  start  with  a  difference  scheme  as  an  alter¬ 
nate  formulation  of  the  problem.  Then  the  convergence  of  the  scheme  is  addressed.  This 
approach  appears  to  be  more  independent  of  the  equation’s  type.  Further,  the  nonlinear 
theory  is  developed  directly  as  opposed  to  being  an  extension  of  the  linear  theory.  Special 
effort  is  made  to  show  that  situations  which  satisfy  Kato’s  hypotheses  also  satisfy  those  of 
Crandall.  But  only  in  [11]  does  Crandall  claim  to  have  convergence  results  comparable  to 
those  of  Kato. 

In  the  present  work,  emphasis  is  on  the  theorems  of  Kato  from  [24]  and  [35].  An 
excellent  text  by  Pazy  [58]  incorporates,  in  well  summarized  form,  a  large  portion  of  Kato’s 
work.  For  convenience,  the  text  of  Pazy  will  be  cited  for  the  introduction  of  terminology 
and  most  of  the  preliminary  results.  Another  text,  which  includes  a  summary  of  Kato’s 

■4 

work  is  [27:pp  237-247]. 

Other  works  to  broaden  perspective  include  [14]  and  [53]. 
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1.4  Outline  of  what  is  to  follow 

The  second  chapter  provides  a  discussion  of  the  abstract  version  of  (9).  The  discussion 
includes  basic  concepts  of  the  abstract  theory  and  theorems  which  are  sufficient  to  establish 
the  existence  and  uniqueness  of  a  solution.  Much  of  the  work  done  in  this  setting  is  key  to 
the  solution  of  the  more  general  problem.  A  linear  space  and  some  operators  are  carefully 
chosen  and  shown  to  have  desirable  properties.  The  abstract  version  of  the  problem  that 
has  been  formulated  is  shown  to  satisfy  the  hypotheses  of  appropriate  theorems.  Kelvin- 
Voigt  damping  is  then  introduced,  leading  to  the  equation 

Htt  4"  (/3lltxx)x x  4*  (QUxx)xx  —  0  (22) 

and  the  individual  cases  of  (3  a  constant  and  /3  dependent  on  x  are  considered. 

In  the  third  chapter,  a  will  be  held  constant  and  time  dependent  ft  will  be  considered. 
This  requires  additional  terminology  and  theory,  which  is  presented.  The  cases  (1  depends 
on  f,  /?  depends  on  x  and  t,  and  finally,  (3  depends  on  u  are  also  considered. 

The  fourth  chapter  considers  a  case  with  a  dependent  on  t.  This  requires  more 
theory,  which  will  be  presented. 

The  final  chapter  summarizes  conclusions  and  recommends  specific  further  work. 

Citations  to  the  literature  are  in  brackets  and  can  thus  be  distinguished  from  refer¬ 
ences  to  equations.  Throughout  the  paper  the  end  of  a  proof  is  indicated  by  a  box,  like 
this.  □ 
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II.  Abstract  formulations 


In  this  chapter  the  basic  cases  are  considered.  There  is  quite  a  bit  of  terminology  to 
introduce.  This  is  done  in  the  context  of  solving  the  constant  coefficient  problem.  Several 
basic  theorems  are  also  presented. 


2.1  Preliminaries  for  the  constant  coefficient  case  with  no  damping 

The  differential  equation  derived  in  Appendix  A  to  approximate  the  unforced  trans¬ 
verse  vibration  of  a  beam  with  no  damping  is 


uti  +  DliaDlu)  =  0 


(23) 


where  Dx  =  ^  and  a  >  0  represents  ET/m.  The  case  for  constant  a  is  not  difficult. 
Indeed,  it  is  easily  solved  directly  by  product  separation  of  variables,  eg  [49:pg  117-120]. 
Nevertheless,  this  simple  case  provides  an  opportunity  to  introduce  terminology,  notation, 
and  a  strategy  for  solution  which  will  also  be  applicable  to  more  complex  cases. 


The  equation  will  be  reformulated  as  a  system  of  first  order  equations  in  t.  Let 

\ 


u{t)  = 


where  the  x  dependence  is  suppressed.  Explicitly,  the  values  of  u(t) 


are  elements  of  a  function  space  which  carries  the  x  dependence.  For  convenience,  the 


usual  notation  will  be  u  = 


l  „  \ 

til 

\“2/ 


.  The  functions  tii  and  u^,  when  evaluated  at  any  point 


in  their  domain,  are  required  to  satisfy  appropriate  boundary  conditions. 


The  single  higher  order  equation  is  now  replaced  with  the  first  order  system 


tij  +  Au  =  0;  t  >  0,  ti(0)  =  tio  (24) 

where  the  vector  variable  u  depends  only  on  t.  In  order  to  keep  the  independent  variables 
clear,  it  is  reemphasized  that  while  U\  and  are  explicit  functions  of  t,  their  values  at  a 
particular  value  of  t  are  functions  of  x.  The  vector  variable  u  and  the  boundary  conditions 


14 


on  its  components  are  given  by 


u  = 


\U2J 

The  operator  A  is  given  by 


;  «i(0)  =  «i'(0)  =  ua(l)  =  ui'(l)  =  0. 


(25) 


A  =  Ai  = 


0  -1 
ctD 4  0 


(26) 


The  subscript  on  A  is  used  to  distinguish  this  particular  operator  from  others  that  will  be 
introduced  later.  When  A  appears  without  a  subscript,  no  specific  operator  is  intended. 
The  operator  D  is  The  symbol  0  is  used  as  a  vector  where  appropriate  without  any 
special  notation.  Also,  the  symbol  1  is  used  for  the  identity  operator  without  special 
notation.  The  domain  of  A\  will  be  denoted  by  D(A\).  The  beam  is  assumed  to  be  finite, 
with  length  1.  The  boundary  conditions  are  as  indicated. 


2.1.1  Some  definitions 

Basic  terminology  is  now  reviewed  to  establish  a  foundation  for  discussion  of  the 
abstract  problem.  The  first  definition  is  that  of  a  semigroup  [58:pg  1]. 

Definition  1  Let  X  be  a  Banach  space .  A  one  parameter  family  S(t),  0  <  t  <  oo,  of 
bounded  linear  operators  from  X  into  X  is  a  semigroup  of  bounded  linear  operators  on  X  if 

1.  S(0)  =  I,  (the  identity  operator) 

2.  S(t+s)  =  S(t)S(s)  for  every  t,s  >  0. 

There  are  several  notions  of  continuity  for  semigroups.  In  this  work,  strongly  continuous 
semigroups  [58:pg  4]  are  used  almost  exclusively. 

Definition  2  A  semigroup  S(t),  0  <  t  <  oo,  of  bounded  linear  operators  on  X  is  a  strongly 
continuous  semigroup  (abbreviated  Cq  semigroup),  of  bounded  linear  operators  if 


Urn  S(t)x  —  x  for  every  x  E  X.  (27) 
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The  terminology  refers  freely  to  semigroups  of  bounded  linear  operators  (even  when  the 
generators  are  not).  This  is  appropriate,  as  the  following  theorem  shows. 

Theorem  3  Let  S(t)  be  a  Cq  semigroup.  There  exist  constants  a>  >  0  and  M  >  1  such 
that 

||S(f)||  ^  Meut  for  0  <  t  <  oo.  (28) 


Proof:  See  [58:pg  4].  □ 

If  w  <  0  the  semigroup  is  said  to  be  uniformly  bounded.  If,  in  addition,  M  =  1  then 
the  semigroup  is  called  a  Co  semigroup  of  contractions. 

It  may  be  possible  to  identify  a  generator  (conceptually  very  much  like  a  derivative) 
for  a  given  semigroup  [58:pg  1]. 


Definition  4  The  infinitesimal  generator  of  a  semigroup  S(t)  is  the  linear  operator  A 
defined  by 


Ax  -  lim 

no 


S(t)x  -  3 


(29) 


whenever  this  limit  exists. 


Notice  the  conceptual  similarity  of  the  generator  to  a  derivative.  Not  all  linear  op¬ 
erators  are  the  generators  of  semigroups.  And,  some  semigroups  may  have  infinitesimal 
generators  which  are  only  defined  on  a  portion  of  their  domain.  The  notion  of  an  infinites¬ 
imal  generator  (often  referred  to  simply  as  a  generator)  is  useful  when  its  domain  is  dense 
in  the  domain  of  definition  of  the  semigroup;  ie,  the  space  X  mentioned  in  the  first  two 
definitions.  The  collection  of  all  operators  A  such  that  —A  is  the  infinitesimal  generator 
of  a  semigroup  on  X  bounded  by  a  particular  /  r,u  pair  is  denoted  by  G(X,  M,u>). 

Later,  it  will  be  useful  to  invert  operators  of  the  form  I  -  \A  for  rational  A  >  0.  It 
is  appropriate  to  identify  the  values  of  A  such  that  this  inverse  exists.  Special  attention  is 
given  to  values  of  A  for  which  the  inverse  is  defined  on  all  of  X.  Some  appropriate  terms 
are  now  introduced  [58:pg  8]. 

Definition  5  The  resolvent  set,  p(A),  of  A  is  the  set  of  all  complex  numbers  A  for  which 
XI  -  A  is  invertible.  That  is,  (A I  -  A)-1  is  a  bounded  linear  operator  on  X. 
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Definition  6  The  family  R(X]A)  =  (XI  -  A)-1,  for  all  X  6  p(A),  of  bounded  linear 
operators  is  called  the  resolvent  of  A. 


2.1.2  Choosing  an  appropriate  space  X 

The  dependent  space  variable  is  suppressed  in  the  notation  of  (2d).  In  the  original 
formulation  u  =  u(t,x),  but  now  the  problem  is  formulated  in  a  linear  space,  X.  The 
elements  of  X  are  functions  of  x  which  satisfy  the  boundary  conditions  of  the  problem. 
Thus  u  is  a  function  of  t  and  assigns,  to  each  value  of  t,  a  unique  element  of  X.  The  x 
dependence  is  hidden  in  the  domain  and  is  not  explicit  in  the  function. 

It  is  reasonable  to  think  of  (24)  as  an  ordinary  differential  equation  over  a  linear 
space.  This  is  referred  to  as  the  abstract  formulation  of  the  differential  equation  and  is 
known  as  an  abstract  Cauchy  problem.  To  be  specific,  choose 


A'  =  iy  = 


H 

^ ,J2  j 


6  II2[ 0,  l]  x  //°[0,  l)  1 7/1  satisfies  boundary  conditions  (25)  >  .  (30) 


Here  IIP  represents  the  standard  Sobolev  space  (see  Definition  S).  Notice  that  X  is  a 
vector  space  with  each  component  coming  from  a  Sobolev  space.  T1  .  boundary  conditions 
are  as  specified  in  (25).  The  following  lemma  allows  an  alternate  description  of  A'. 


Lemma  7  If  y  e  II 2  and  7/(0)  =  y(l) -=  y'( 0)  =  2/'(l)  =  0,  then  y  €  7/q. 


Proof:  See  Appendix  D.  □ 

The  particular  boundary  conditions  under  consideration  make  it  appropiiale  to  apph 
Lemma.  7  and  describe  X  in  the  abbreviated  form 


( 

\ 

j 

Vi 

\y  = 

{ 

\  y2 

) 

S  /io[0,l]  x  7/°[0,l] 


(31) 


The  domain  [0,1]  will  often  be  omitted  in  the  sequel.  However,  it  is  intended  foi  the 
entirety  of  the  paper. 


17 


t 

For  convenience,  the  standard  Sobolev  norm  is  defined.  At  this  point  the  Sobolev 
norm  is  only  used  to  identify  functions  which  are  components  of  elements  in  the  set  X. 
Points  in  the  linear  space  X  will  have  a  norm  to  be  specified  momentarily. 

Definition  8  The  standard  (£%  style)  Sobolev  space  Hp  consists  of  those  functions  u  such 
that 

DP‘<C<co  (32) 

t=0 

where  D°u  =  u  and  the  norm  is  given  by 

ll«I»  =  (pit'll!,)  (33) 

The  points  (pairs  of  functions)  in  X  have  been  identified  and  it  is  easy  to  see  that 
they  form  a  linear  space.  It  is  desirable  to  have  a  Banach  space,  and  hence,  a  norm  must 
be  specified  for  the  elements  of  X.  Alternatively,  it  is  appropriate  to  specify  an  inner 
product  for  the  elements  of  X  and  let  the  norm  be  the  one  naturally  induced  by  it.  An 
inner  product  will  be  specified  and  then  it  will  be  shown  that  A  is  a  complete  space  under 
the  induced  norm.  A  preliminary  lemma  is  required. 

Lemma  9  If  z  is  absolutely  continuous  on  [0,1]  and  z'{x)  —  0  almost  everywhere,  then  z 
is  constant. 

Proof:  See  [65:pg  105].  □ 

An  inner  product  is  now  presented.  This  inner  product  has  been  used  previously  in 
[I2:pg  16]  and  [35:pp  144,  147]. 

Theorem  10  The  expression 


where  a  is  a  positive  constant,  defines  an  inner  product  on  X. 
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Proof:  This  is  shown  to  be  an  inner  product  in  a  straightforward  fashion  from  the  defini¬ 
tion  (to  review  the  definition,  see,  for  example  [54:pg  272],  [65:pg  210],  or  [73:pp  39-40]). 
The  necessity  of  having  the  boundary  conditions  in  the  definition  of  X  becomes  clear  in 
the  details. 

The  only  issue  in  doubt  is  that  of  positive  definiteness.  Suppose  that 


This  requires  each  integral  to  be  zero,  and  hence  u\"  and  u 2  to  be  zero,  at  least  almost 
everywhere  (a.e.).  Since  functions  which  only  differ  on  a  set  of  measure  zero  are  treated 
as  identical,  it  is  clear  that  u2  is  zero.  If  the  boundary  conditions  force  u\  =  0,  then  the 
proposed  inner  product  is  legitimate. 

Lemma  9  is  now  applied  to  show  that  if  z  £  H2  and  z"  =  0  a.e.,  then  z  =  0.  Indeed, 
for  z  to  be  an  element  of  H 2  requires  Jq(z")2  dx  to  exist.  Since  the  domain  of  integration 
is  bounded,  Jq  z"  dx  exists.  It  follows  that  z'  is  absolutely  continuous  (a.e.).  Since  z'  is 
a.e.  and  z"  —  0  a.e.,  the  lemma  says  z'  is  constant.  The  boundary  conditions  on  z'  force 
the  constant  to  be  zero.  Thus,  z  is  constant.  The  boundary  conditions  on  z  force  this 
constant  to  be  zero  also.  Apply  the  lemma  to  u\  and  positive  definiteness  is  clear. 

This  completes  the  proof  of  the  theorem.  □ 

The  norm  chosen  for  X  is  ||  •  ||^  =  (•,  -)1/2.  The  set  X  with  this  norm  is  the  Banach 
space  that  will  be  referred  to  in  the  sequel.  (Since  the  norm  comes  from  an  inner  product 
it  would  be  appropriate  to  call  X  Hilbert  space.  However,  the  more  general  term  Banach 
space  will  usually  be  used.) 

2.1.3  Completeness  of  X 

It  is  important  that  X  be  complete.  This  issue  will  now  be  addressed.  Several 
preliminaries  are  needed.  As  a  matter  of  notation,  ||  •  denotes  the  supremum  norm. 
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Lemma  11  If  z  e  H 2  then  \\z\\l2  <  ||z||oo- 


Proof:  It  is  important  for  the  proof  that  the  function  have  a  bounded  domain  of  definition. 


/  n  „  s 

1/2 

^  j  z2(x)  dx 

) 

If1  0  > 

1/2 

(/„ 

) 

IWI»-  D 

(36) 


Lemma  12  If  z  e  H 2  then,  pljoo  <  |^r(0)[  +  \W\\l2 


Proof:  Since  z  6  II2,  z'  is  square  integrable.  Since  the  domain  is  a  bounded  interval,  it 
follows  that  z'  is  integrable.  Then  [65:pg  101]  z  is  given  (a.e.)  by 


It  follows  that 


zfx)  =  2(0)  +  J  z'(x)di 


K*)l  = 


< 

< 

< 

< 


2(0)  +  f  z'(x)dx 
Jo 

|  ^(O)  |  +  |  jf  z'{x)dx 
|  *(°)  |  +  J  |  z'{x)  |  dx 
I  ,(0)  |  +  j  |  z'(x)  |  dx 

I  z(0)  I  +M*- 


(37) 


(38) 


The  Schwarz  inequality  has  been  used  in  the  last  step.  Since  the  right  hand  side  is  free  of 
x  it  follows  that 

Moo  <  I  *(o)  |  +  IMli,  (39) 

as  desired.  □ 

This  lemma  has  an  obvious  corollary. 

Corollary  13  If  z  £  II2  and  2(0)  =  0,  then  [|^||oo  <  IN'Ilz** 
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A  similar  lemma  and  corollary,  for  z '  instead  of  2,  follow  immediately.  When  2(0)  = 
z'( 0)  =  0,  these  combine  to  give 

Mh  <  IMI00  <  \\A\L2  <  Halloo  <  \\z"\\l>-  (40) 

When  the  above  lemmas  apply,  the  following  composite  result  is  available. 

Lemma  14  For  the  given  norm  on  X,  ||2j|#s  <  ^ 

Proof:  This  is  proved  by  using  the  previous  results. 

Ma?  = 

< 


This  concludes  the  preliminaries  for  an  argument  on  the  completeness  of  X. 


(MIL  +  ML  +  Ml!,) 1/2 
(3IN"IIL)i/2 


3 


(!) 


\°  / 
( 


1/2 

x) 


2 

l  0 


(41) 


\°/ 


Theorem  15  The  space  X,  with  the  inner  product  introduced  in  Theorem  10,  is  complete. 


Proof:  Let  {j/n}  be  a  Cauchy  sequence  in  X.  Recall  that  elements  of  X  have  two  compo¬ 
nents.  To  be  explicit, 


{fc}  = 


' 

(  M 

< 

* 

(42) 


This  will  converge  to  a  point  y  = 


V 

w 


6  X  if  and  only  if  the  first  components  of  the 


given  sequence  converge  to  y\  G  Hq,  and  the  second  components  converge  to  t/2  6 
Convergence  will  be  addressed  for  each  component  separately.  Keys  to  the  strategy  are 
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the  relationship 


/  \ 

(  \ 

(  n  \ 

Vln  I 

yin 

0 

= 

+ 

\ V2n  I 

{ 0  J 

V y2n  / 

and  the  fact  that  (for  the  specified  norm  on  X)  that 


/ 

\ 


/ 

(  n  \ 

yin 

+ 

l 0 

) 

X  . 

\y2n) 

(43) 


(44) 


Since  the  left  hand  side  forms  a  Cauchy  sequence,  it  is  clear  that  each  term  on  the  right 
hand  side  must  form  a  Cauchy  sequence.  Now  focus  attention  on  the  second  term. 

For  each  n,  y2n  €  H°.  Furthermore, 


=  lljfejtfo- 


(45) 


Hence,  {t/2ri}  is  Cauchy  in  ||  •  \\H0.  But,  H°  is  complete  by  definition,  eg  [l:pg  44].  Hence, 
there  exists  some  y2  G  H°,  y2n  -*■  t/2*  This  is  the  desired  y2. 


The  case  of  y\  is  nearly  as  simple.  Each  element  of  {yin}  is  in  II q,  and  the  strategy 
is  to  find  y\  in  the  complete  space  Hq  such  that  yin  ->  yi.  But,  this  depends  on  {y\n} 
being  a  Cauchy  sequence  in  the  norm  of  Hq.  This  follows  immediately  from  Lemma  14. 
Thus  {j/in}  is  Cauchy  as  a  sequence  in  the  complete  space  Hq  and  converges,  say  to  y\. 


It  is  easy  to  see  that 
the  proof  of  the  theorem. 


/  \ 


y  1 

U; 


G  X  and  it  follows  that  X  is  complete.  This  completes 

□ 


It  is  worth  noting  that  the  boundary  conditions  were  important  in  these  arguments. 
Extension  of  the  result  to  other  boundary  conditions  is  not  trivial.  Additional  comments 
on  this  point  are  in  Appendix  B. 


2.1.4  Identifying  the  domain  of  A\ 


It  is  important  to  identify  the  portion  of  X  on  which  A\  makes  sense  and  lias  its 
image  again  in  X.  Let 

£  H4  X  H$  |  y\  satisfies  boundary  conditions  (25)  j  .  (46) 

Notice  that  the  second  components  satisfy  the  boundary  conditions.  This  is  cxpicssod  in 
the  11$  notation  and  omitted  from  the  specifier.  Recall  that  the  boundan  conditions  weto 
specified  in  (25). 

Several  comments  are  appropriate  to  prepare  to  establish  that  the  domain  of  .-lt  is 
dense  in  _Y.  First,  note  that  it  is  sufficient  to  consider  the  components  individually.  Second, 
in  the  way  of  notation,  the  collection  of  all  functions  which  are  infinitely  difici callable  and 
have  compact  support  on  (0, 1)  will  be  denoted  .  (It  would  be  nioie  st  audit  id  to  wiile 
(0.1)  so  there  is  some  abbreviation  here.)  Third,  the  collection  of  I'iiik lions  in  //'’ 
which  satisfy  the  boundary  conditions  will  be  denoted  //$. 

Theorem  16  For  A\  in  (26),  D(A\)  is  dense  in  X . 

Proof:  The  skeleton  of  an  appropriate  argument  has  been  presented  in  [I2:pg  19].  For 
denseness  of  the  first  component  it  is  sufficient  to  establish  that  H £  D  II J.  Note  that 

C0~  C  Hi  c  III  C  III  =  Q?  (‘17) 

where  the  closure  is  with  respect  to  the  H 2  norm.  The  first  two  inclusions  are  clear.  The 
last  inclusion  follows  from  Lemma  7.  The  equality  holds  by  definition,  eg  [l:pg  45]. 

Note  that  Hq  =  11$  and  that  II q  is  the  first  component  of  D(A\).  Also,  11$  D  11 J 
and  hence  Hq  D  II J.  Denseness  of  the  first  components  is  now  clear. 

For  the  second  component  it-  is  sufficient  to  establish  that  //?  3  // 0  Consider 

C0~  C  III  C  II2  C  H°  =  u  C  C?  (48) 
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where  the  closure  is  with  respect  to  the  L2  norm.  The  last  containment  follows  from 
Theorem  2.19  of  [l:pg  31]  (see  also  [20:pg  253]  or  [54:pg  592]).  Now  the  denseness  is  clear. 

Since  each  component  of  D(A\)  is  dense  in  the  corresponding  component  of  X,  it  is 
clear  that  D(A\)  is  dense  in  X.  □ 

Corollary  17  The  inclusions  Hq  C  Hq  C  H°  are  dense. 

2.1.5  Special  properties  of  A\  and  its  adjoint 

Later,  there  will  be  a  need  for  the  adjoint,  AJ,  of  A\.  Initial  definitions  and  com¬ 
putations  are  presented  now.  See  [54:pp  352,  527],  [63:pp  201,  215],  or  [73:pg  196]  for 
details. 


Definition  18  The  operator  A*  is  called  the  adjoint  of  the  operator  A  if(Au,  v )  =  (u,  A*v ) 
for  all  u,v  €  D(A). 

The  following  derivation  identifies  the  operator  A*  for  A\  defined  by  (26). 


/ 

Ai 

\ 

\U2J 

v2JJ 


-a  f  D2u2D\  dx  +  a  f  v2DAu\  dx 
Jo  Jo 

-aDu2D2v\  |o  -fa  J  Du2D3vi  dx 

+ av2D3ui  |o  -a  f  Dv2D3u\  dx 
Jo 

au2D3v\  |o  -a  /  u2D4v  1  dx 
Jo 

-ocDv2D2ui  |o  +a  J  D2v2D2u\  dx 


«i 


V2 


\\U2J  \ 


( 

(  \ 

/ 

U\ 

J 

\ 

Kui) 

\ 

-aD4v  1  J  J 


0 

-aD4 


\ 

(  \\ 

'  ®1 

) 

l V2)> 

(49) 
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where  all  of  the  boundary  terms  are  zero  and  are  dropped  after  their  first  appearance. 

It  is  interesting  that  AJ  =  —A\.  When  this  occurs,  the  operator  A\  is  called  skew- 
adjoint.  Since  D(A\)  =  D(-A\)  it  follows  that  jD(AJ)  =  D(Ai).  Thus,  since  A\  is  densely 
defined,  so  is  A\.  (Actually,  it  is  not  quite  this  straightforward.  Special  care  is  required  in 
describing  the  domain  of  an  unbounded  operator.  For  A\  it  is  more  complete  to  argue  as 
follows:  Let  y  £  D{A\)  be  given.  The  goal  is  to  show  that  y  £  D(A\).  This  requires  the 
existence  of  some  z  £  X  such  that  for  every  a:  £  D(Ai),(Aix,y)  =  (*,2).  Let*  £  X>(Ai)be 
given.  By  direct  computation  ( Aix,y )  =  (s,-Aiy).  Choose  2  =  -A\y.  Since  y  £  D(A{) 
and  jD(Ai)  =  D(-A\ )  it  follows  that  2  £  X.  Thus  y  £  D(A\)  and  D(Ai)  C  -D(AJ).  The 
argument  here  is  reversible  and  the  containment  goes  both  ways.  The  equality  of  the  two 
domains  is  now  clear.)  The  issue  of  A\  actually  being  defined  on  the  dual  space  has  been 
ignored  since  X  is  a  Hilbert  space  and  is  identified  with  its  dual  (eg  [73:pg  91]). 

Another  useful  concept  is  that  of  the  dissipative  operator  [58:pg  14]. 

Definition  19  A  linear  operator  A  is  dissipative  if  ||(A/ -  A)u||  >  A||u||  for  all  A  >  0  and 
all  u  £  D(A ). 

The  idea  of  an  operator  being  dissipative  is  quite  simple.  If  Au  +  u  is  in  some  way 
less  than  u,  then  it  is  quite  reasonable  to  call  A  dissipative.  If  u  were  a  vector  it  would 
be  reasonable  to  think  of  Au  as  having  a  component  in  the  direction  of  -u.  This  could 
be  validated  algebraically  by  considering  the  inner  product  and  requiring  ( Au,u )  <  0.  For 
a  definition  that  is  acceptable  in  a  Banach  space  it  is  reasonable  to  require  u  -  Au  to  be 
greater,  in  some  sense,  than  u.  If,  for  every  A  >  0,  it  happens  that  Xu  -  Au  is  bigger  than 
Xu,  then  it  is  certainly  reasonable  to  call  A  dissipative.  This  is  just  what  the  definition 
does. 

The  following  lemma  shows  that,  in  a  Hilbert  space,  the  desired  implication  holds. 
Lemma  20  If  (Au,  u)  <  0  for  every  u  £  D(A),  then  A  is  dissipative. 

Proof:  Basic  definitions  are  sufficient  to  establish  the  implication. 

\\(XI  -  A)u||  =  ((XI  -  A)u,  (XI  -  A)u)1/2 
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=  (A  Iu  -  Au ,  A  In  -  Au)1/2 
=  ((Xu,  Xu)  -  (Xu,  Au)  -  (Au,  Xu)  +  (Au,  Au))1/2 
=  .  (X2(u,  u)  —  2X(Au,  it)  +  (Au,  Au))1/2 

>  (A2(u,  u)  -  2X(Au,  u))1/2 

>  (X2(u,u))1/2 

=  X\\u\\.  □  (50) 


It  is  useful  to  establish  that  the  operators  A\  and  A\  are  dissipative.  First,  consider 


Al 


Lemma  21  The  operator  A\  in  (26)  is  dissipative. 


(  /  \  /..  w 

Ui 


Proof:  Proof  is  by  direct  computation. 

=  -a  f  D2uiD2u\  dx  +  a  [  u2D'lui  dx 
Jo  Jo 

=  —aDu2D2u\  |o  +c*  f  Du2D3ui  dx 
Jo 


Ai 


\U 2  ) 


«i 

\u2J 


+au2D3ui  |o  -a  [  Du2D3U\  dx 
Jo 


=  0. 


(51) 


where  the  boundary  terms  are  again  zero.  This  is  sufficient,  by  Lemma  20,  to  establish 
that  A  is  dissipative.  □ 

It  is  immediately  clear  that  —Ai,  A\,  and  -A\  =  (-Aj)*  are  all  dissipative.  This  is 
a  convenience  of  having  an  equality  when  an  inequality  suffices  in  (51).  Of  course  this  was 
not  accidental,  the  choice  of  inner  product  made  it  happen.  It  is  interesting  to  note  how 
such  an  inner  product  is  chosen. 
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If  the  usual  L2  inner  product  is  used,  then  the  following  computations  are  sufficient 
to  check  whether  A\  is  dissipative. 


This  quantity  must  be  less  than  or  equal  to  0. 

One  strategy  is  to  redefine  the  inner  product  so  that  everything  cancels.  Possible 
replacements  for  —  Jq1  u^ui  dx  would  be 


-  [  av,2D4Ui  dx, 
Jo 

/  aDu2D3Ul  dx,  or 
Jo 

-  /  ctD2u2D2Ui  dx. 
Jo 


(53) 


Each  of  these,  after  the  appropriate  application  of  integration  by  parts  would  show  Ai 
to  be  dissipative  by  making  the  inner  product  exactly  zero.  The  last  candidate  in  (53) 
corresponds  to  taking  a  times  the  L2  inner  product  of  the  second  derivatives  of  the  first 
components.  It  has  been  verified  in  Theorem  10  that  such  a  replacement  leads  to  a  legiti¬ 
mate  inner  product. 


2.1.6  Closedness  of  A\ 

It  is  also  important  to  establish  that  the  operators  A\  and  Ai  are  closed.  A  definition 
is  first  given.  See  [54:pp  241,  529],  [63:pg  300],  or  [73:pg  77]. 

Definition  22  A  linear  operator  A  is  closed  if,  when  any  pair  of  sequences  {xn}  C  D(A)  C 
X  and  {yn}  =  {Aa;n}  C  X  both  converge,  say  xn  -*  x  G  X  and  yn  — >  y  G  X  then  x  G  D(A ) 
and  Ax  =  y. 
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Lemma  23  If  A  is  a  closed  operator,  then  so  is  —A. 


Proof:  This  is  clear  from  the  linearity  of  the  space  and  the  equality  D(A)  =  D(—A).  □ 

The  proof  of  the  theorem  that  will  be  used  to  establish  closedness  of  A\  requires  a 

lemma. 

% 

Lemma  24  Let  (•,  •)  denote  an  inner  product  on  a  linear  space  X.  For  fixed  x  6  X,  let 
f(y)  =  (x,y).  Then  f  is  continuous.  That  is,  the  inner  product  is  a  continuous  function 
of  its  arguments. 

Proof:  See  [20:pp  179-180]  or  Appendix  D.  □ 

There  will  be  several  occasions  to  use  the  next  result  [54:pg  529].  It  is  actually  a 
special  case  of  the  theorem  cited.  An  additional  hypothesis  that  A  be  closable  and  an 
additional  result  that  jD(A*)  is  dense  are  omitted.  (The  proof  of  the  omitted  portion  is 
relegated  to  guided  exercises  in  [54:pp  531-532].  A  different  proof  is  in  [25:pg  172].  See 
also  [63:pp  299-300].) 

Theorem  25  Let  A  be  a  densely  defined  linear  operator  on  a  Hilbert  space  H.  Then  the 
adjoint,  A*,  of  A  is  a  closed  operator. 

Proof:  Let  u *  — ►  u*  and  A*u*n  — >  w*  with  u*  G  D(A*)  for  all  n.  By  the  definition  of 
adjoint,  ( Au ,  «*)  =  (u,  A*u * )  for  all  u  G  D(A).  From  the  continuity  of  the  inner  product  it 
follows  that  ( u,w *)  =  (Au,m*).  But,  again  recalling  the  definition  of  adjoint,  u*  G  D(A*) 
if  there  is  some  z  e  X  such  that  (Au,u*)  =  ( u,z ).  When  this  occurs  2  =  A*u*.  Clearly, 
z  —  w*  is  the  element  needed  and  A*u*  =  w*.  Thus  u *  G  D(Am). 

This  establishes  that  A*  is  closed.  □ 

Lemma  26  The  operator  —Ai,  in  (26),  is  closed. 

Proof:  Recall  that  A\  =  —  A\  and  that  D(A\)  =  D(Ai).  Thus,  since  D(A\)  is  dense  in 
X  so  is  D(A\).  (A  more  general  approach  would  be  to  apply  Theorem  1.4.5.C  of  [58:pg  15] 
and  Theorem  7.10.3  of  [54:pg  529],  but  this  is  not  necessary  in  the  present  case.)  Then, 
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by  Theorem  25,  A\*  is  a  closed  operator.  But  AJ*  =  A\  (It  is  almost  trivial  to  verify  this 
directly.  It  is  not  surprising  in  the  Hilbert  space  setting,  for  example  see  [54:pg  353]  for  the 
case  of  bounded  A.)  so  A\  is  a  closed  operator.  It  is  immediate  that  -A\  is  also  closed.  □ 

2.1.7  The  generator  of  a  Co  semigroup  of  contractions 

It  is  important  that  the  operator  -A\  be  the  generator  of  a  Co  semigroup  of  con¬ 
tractions.  This  is  an  immediate  result  of  tha  next  theorem. 

It  is  useful  to  note  that  the  literature  is  not  uniform  in  the  use  of  A  and  —A.  This 
result'  from  the  form  in  which  the  original  differential  equation  is  written.  That  is  to  say 
ut  +  Au  =  0  and  ut  =  Au  correspond  to  opposite  sign  conventions.  Also,  some  authors 
use  accretiveness  instead  of  dissipativeness  which  leads  to  opposite  sign  conventions.  It  is 
important  to  be  self-consistent.  Nevertheless,  it  frequently  requires  a  conscious  effort  to 
keep  this,  somewhat  annoying  detail,  straight. 

The  next  theorem  provides  a  useful  tool  for  establishing  that  certain  operators  are 
the  infinitesimal  generators  of  Cq  semigroups. 

Theorem  27  Let  A  be  a  densely  defined  closed  linear  operator.  If  both  A  and  A*  are 
dissipative,  then  A  is  the  infinitesimal  generator  of  a  Co  semigroup  of  contractions  on  X. 

Proof:  See  [58:pg  15]  or  [40:pg  87,  Theorem  4.4].  □ 

Corollary  28  The  operator  -A\  tn  (26)  is  the  generator  of  a  Co  semigroup  of  contrac¬ 
tions. 

Proof:  It  has  been  established  in  the  preceeding  subsections  that  the  hypotheses  of  the 
theorem  are  satisfied.  Hence  -A\  is  the  generator  of  a  Co  semigroup  of  contractions.  □ 

2.2  Remarks 

It  is  interesting  to  note  that  A\  is,  in  fact,  the  generator  of  a  group.  (Use  Theorem 
27  on  A*  and  A**.  Also,  see  [19:pg  22,  2.16,  pg  32  Theorem  4.7],  and  [58:pg  22,  Theorem 
1.6.5  and  pg  41,  Theorem  1.10.8].)  This  fact  is  not  of  any  immediate  interest.  However, 
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it  is  appropriate  to  note  that  this  feature  can  be  useful.  For  an  example,  see  [43:pg  745] 
where  it  is  used  in  a  controls  problem. 

A  distinctive  feature  of  this  work  is  the  attention  to  hyperbolic  problems.  This  means 
that  —At  is  the  generator  of  a  Co  semigroup  but  not  necessarily  of  an  analytic  semigroup 
[35:pp  128-129].  If  it  were  the  generator  of  an  analytic  semigroup  then  the  theory  of 
parabolic  equations  would  apply,  eg  [16:pg  108].  It  is  appropriate  to  verify  that  the  current 
problem  does  not  fit  into  the  more  specialized  category.  This  is  proved  in  Anpendix  C  to 
avoid  too  much  of  a  distraction  at  this  point. 

2.3  Finishing  up  for  constant  a 

It  is  time  to  complete  the  constant  coefficient  problem.  The  preliminaries  have  been 
rather  complete  and  have  included  verification  of  the  hypotheses  of  the  existence  theorem 
to  be  presented  now.  The  next  theorem  guarantees  the  existence  of  a  unique  C1  solution 
to  (24)  for  any  uq  €  D{A).  A  solution  is  an  X  valued  function  u(t)  that  is  continuous  on 
[0,oo),  continuously  differentiable  on  (0,oo),  has  u{t )  €  D(A)  for  all  t  >  0,  and  satisfies 
the  differential  equation  (24)  for  all  t  >  0  [58:pp  105, 139]. 

2.3.1  The  existence  theorem 

Theorem  29  Let  the  operator  -A  be  the  infinitesimal  generator  of  a  Co  semigroup  S(t). 
Then  (24)  has  a  unique  solution  u(t),  which  is  continuously  differentiable  on  [0,  oo),  for 
every  initial  value  in  D(-A). 

Proof:  See  [58:pp  102-104].  □ 

Theorem  30  The  abstract  Cauchy  problem  (24),  with  A.  replaced  by  A\,  has  a  unique 
classical  solution. 

Proof:  To  see  that  the  hypotheses  of  the  preceeding  theorem  are  satisfied,  it  suffices 
to  review  Theorem  16  and  Theorem  27  (which  in  turn  depends  on  Lemma  21  and  the 
comments  following  it).  This  theorem  now  follows  from  Theorem  29.  □ 
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This  completes  the  existence  and  uniqueness  argument  for  constant  a  >  0.  The 
solution  is  given  by 

u(t)  =  S(t)u0  (54) 

where  S(t)  is  the  semigroup  generated  by  -A\. 

2.3.2  Continuous  dependence 

It  may  be  that  «o  is  not  known  exactly,  or  even  if  it  is  known,  it  may  be  that  it 
cannot  be  represented  exactly.  In  either  of  these  cases  it  becomes  important  to  know 
whether  small  changes  in  uo  lead  to  only  small  changes  in  u(t).  This  is  usually  referred  to 
as  continuous  dependence  of  the  sou  '‘ion  on  the  data. 

.  For  the  current  problem,  S(t)  is  a  Co  semigroup  of  contractions.  This  means  that 
||S||  <  1.  Hence,  the  following  theorem  is  straightforward. 

Theorem  31  The  unique  classical  solution  of  (24),  with  A  replaced  by  A\,  depends  con¬ 
tinuously  on  the  initial  data. 

Proof:  Suppose  «o  is  an  initial  condition,  possibly  different  from  wo-  Let  e  >  0  be  given. 
For  t  >  0,  let  u(t)  =  S(t)uo.  Then 

\\u(t)  -  u(t)|U  =  H-S'(t)t£o  —  ■S’CO^oIIa' 

=  ||5(t)(u0  -  fio)iU 

<  ||5,|U  ||«o  -  &o||Ar 

<  \\uo  -  uo|U  •  (55) 

Choose  6  =  e  and  the  continuity  of  the  dependence  is  established.  □ 

2.3.3  Constructing  a  solution 

Before  consideration  is  given  to  generalizations  of  the  problem,  it  is  appropriate 
to  consider  how  the  solution,  whose  existence  has  just  been  established,  can  actually  be 
obtained. 
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If  the  operator  —A  is  bounded,  and  the  infinitesimal  generator  of  a  Co  semigroup, 
then  the  solution  to  (24)  is  given  by  u(t)  =  e~tAuo.  Since  A\  in  (26)  need  not  be  bounded, 
limit  definitions  (with  the  same  flavor  as  that  of  e~iA )  must  be  carefully  analyzed  to  assure 
convergence.  But,  the  conditions  which  assure  that  -A\  is  the  generator  of  a  Co  semigroup 
do  just  that.  When  it  exists,  e~tA  is  the  semigroup  S(t).  This  motivates  the  terminology 
for  -A  as  the  generator  of  S(t).  The  solution  always  has  the  form  u(t)  =  S(t)u o.  The  key 
issue  in  applications  is  to  determine  when  5(f)  can  be  determined  from  A.  When  S(t)  can 
be  obtained,  the  first  component  of  u(t)  is  the  solution  to  the  original  problem. 

Here  is  one  way  of  obtaining  a  semigroup  [58:pg  33]. 


Theorem  32  Let  S(t)  be  a  Co  semigroup  on  X.  If  A  is  the  infinitesimal  generator  of  S(t) 
then 


\“n 

n  (n  Y 

x  =  hm 

-R  -f  :A) 

/  n— *co 

[t  \t  J 

(56) 


for  x  €  X,  and  the  limit  is  uniform  in  t  on  any  bounded  interval. 


Proof:  See  [58:pp  34-35],  A  more  general  development  is  presented  later  in  the  same 
reference,  [58:pp  89-92].  □ 

In  the  event  that  the  operator  A  is  either  not  known  precisely  or  cannot  be  represented 
precisely,  it  is  of  interest  to  know  whether  the  semigroup  generated  by  an  approximate  A  is 
close  to  the  semigroup  that  would  have  been  generated  by  an  exact  A.  This  has  the  same 
flavor  as  the  continuous  dependence  considered  above.  Indeed,  it  would  be  appropriate  to 
ask  whether  the  semigroup  depends  continuously  on  the  generator. 

Results  of  this  type  are  obtained  by  considering  a  sequence  An  which  is  assumed  to 
converge  to  A  in  an  appropriate  sense.  The  question  is  whether  the  corresponding  semi¬ 
groups  Sn{t)  converge  to  the  semigroup  generated  by  A.  Formal  results  in  this  direction 
often  consider  also  whether  A  depends  continuously  on  S.  Results  of  this  type  are  referred 
to  as  Trotter  or  Trotter-Kato  theorems.  (Trotter  presented  the  pioneering  work  in  the 
linear  case.  Kato  corrected  an  error  in  the  published  proof.) 
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For  results  applicable  to  the  current  problem,  see  [58:pp  35,  84-89],  [19:pp  44,  48], 
or  [73:pg  269].  The  results  are  that  three  distinct  convergences  axe  in  fact  equivalent:  the 
operators,  the  semigroups,  and  the  resolvents  of  the  operators. 

The  works  of  Crandall,  eg  [10]  and  [11],  provide  an  alternate  approach.  The  alternate 
approach  will  not  be  pursued  here.  However,  it  is  appropriate  to  note  that  it  started  with  a 
fundamental  paper  in  1971  [8].  The  key  theorem  has  a  somewhat  sketchy  proof.  A  detailed 
version  is  provided  as  Appendix  E. 


2.4  Spatially  dependent  coefficient 

Consider  utt  +  D2(a(x)D7u)  =  0.  Assume  that  a  is  continuously  differentiable  with 
respect  to  x  and  for  some  am,-„  >  0,  a(x)  >  amin  >  0  for  all  a;  G  [0, 1].  Choose  the  space  X 
and  its  inner  product  as  before.  For  the  validation  of  the  inner  product  to  go  as  before,  the 
strict  inequality  in  the  requirement  a(a:)  >  0  is  necessary.  This  is  because  (35)  must  hold 
with  a  inside  the  integral,  otherwise  the  argument  is  unchanged.  This  time  the  operator 
is  defined  as  follows: 

/ 

0  -1 
D\a{-)D 2)  0 

Clearly  Ai  is  a  linear  operator.  Note  that  D{Ai)  =  D(A\).  In  order  to  determine  let 
u  G  D(A2)  be  given,  then 


“«2 

D2(a(-)D2)Ul 


A2 


(57) 


=  -  /  a(x)D2v,2D2vi  dx  +  f  V2D2(a(x)D2U\)  dx 
Jo  Jo 

=  -a{x)Du2D2v\  |o  +  J  DuiD(a(x)D2vi)  dx 

+V2D(a(x)D2U\)  |o  -  /  Dv2D(a(x)D2Ui)dx 
Jo 

=  U2D(a(x)D2v{)  |o  -  /  U2D2(a(x)D2vi)dx 

■  Jo 

-Dv2d{x)D2u\  !o  +  /  ot{x)D2u\D2V2dx 
Jo 
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// 


til 

ti2 


0  1 
-D\a{x)D2)  0 


(  \\ 

Vi 

/ 

\v2h 

(58) 


where  all  of  the  boundary  terms  are  zero. 


Observe  that  A\  = 


\ 


and  hence,  as  before  A\  =  -A2.  Since  the 


0  1 
-D2(a(-)D2)  0  j 

steps  are  reversible,  it  is  easy  to  see  that  (j42)*  =  -A\  =  (-A2)*  =  A2.  By  Theorem  25, 
-A2  is  closed. 


Next,  it  is  necessary  to  show  that  -A2  is  dissipative.  For  u  G  £>(A2) 


( 


— A2 


ti2 


til 


-K2(a(-)D2tii)  j  \  ti2 


=  f  a(x)D2U2D2uidx  -  f  u,2D2(a(x)D2ui)dx 

Jo  Jo 

=  a(x)Dv,2D2ui  |o  -  /  Du2D(<x(x)D2u\)dx 
Jo 

-v,2D(a(x)D2 tii)  |o  +  J  Du2D{a{x)D2ui)dx 


=  0 


(59) 


where  again  the  boundary  terms  are  each  zero. 

Clearly  A2  is  also  dissipative  in  this  case.  But  A2  =  (-A2)*  and  thus  (-A)*  is 
dissipative.  Now  all  of  the  theorems  apply  as  in  the  case  of  constant  a  and  guarantee  the 
existence  of  a  unique  solution  to  the  differential  equation. 


2.5  Constant  coefficient  Kelvin-Voigt  damping 

Consider  utt  +  (Jutxxxx  +  oluxxxx  =  0  with  ti(0,z)  =  ti0,  tit(0,a;)  =  uio  with  the 
boundary  conditions  as  before.  Throughout  this  section  a  is  a  positive  constant.  To  begin, 
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- D2(a(x)D 2)  0 


u2 


\V2J 


where  al"  of  the  boundary  terms  are  zero. 

/  0  1 
-D2(a(-)D2)  0 

steps  are  reversible,  it  is  easy  to  see  that  (AJ)*  =  -A2  =  (—A2 )*  =  A2.  By  Theorem  25, 
— A2  is  closed. 


Observe  that  A £  = 


and  hence,  els  before  A2  =  —  A2.  Since  the 


V 


Next,  it  is  necessary  to  show  that  -A2  is  dissipative.  For  u  €  D(A2) 


( 

(  \ 
U\ 

\ 

—A2 

> 

{ 

J 

\U2) 

/ 

u; , 

^  ~D2(a(-)D2u,) 


U 1 
«2 


\ 

) 


f  a(x)D2u2D2u\  dx  -  f  u2D2(a(x)D2ui)dx 
Jo  Jo 

a(x)Du2D2Ui  |o  -  /  Du2D(a(x)D2ui)  dx 
Jo 

-u2D(a{x)D2u{)  |o  +  /  Du2D{a{x)D2u\)dx 
Jo 


(59) 


where  again  the  boundary  terms  are  each  zero. 

Clearly  A2  is  also  dissipative  in  this  case.  But  A2  =  (—A2)*  and  thus  (-^4)*  is 
dissipative.  Now  all  of  the  theorems  apply  as  in  the  case  of  constant  a  and  guarantee  the 
existence  of  a  unique  solution  to  the  differential  equation. 


2.5  Constant  coefficient  Kelvin-Voigt  damping 

Consider  utt  +  /3utxxx x  +  auxxxx  =  0  with  u(0,a;)  =  u0,  ut(0,a;)  =  «<o  with  the 
boundary  conditions  as  before.  Throughout  this  section  a  is  a  positive  constant.  To  begin, 


(64) 


=  pDu2D2u 2  \l  -p  J\d2u2)2  dx 
<  0 


since  P  >  0. 

The  next  task  is  to  identify  (-^3)*. 


( -M 


I  ux 


\ 


\  \U2J 


V2  /  / 


u2 

—aD4u\  —  PD4u2 


\ 

(  \\ 

Vi 

> 

/ 

kV2)j 

(65) 


=  a  [  D2u2D2V\dx  —  a.  [  v2D4u\dx  —  P  [  v2D4u2dx 
Jo  Jo  Jo 

=  aDu2D2v\  -a  f  Du2D3v\  dx  -  av2D3U\  |q 
Jo 

+a  J  ,Uv2D3Ui  dx  -  Pv2D3u2  |q  +/?  J  Dv2Dzu2  dx 
-  ~au2Dzv\  |o  +<*  J  u2D4V\  dx  +  <xDv2D2U\ 

rl  r\ 

-a  D2v2D2u\  dx  +/3Dv2D2u2  |J  -/?  /  D2v2D2u2dx 
Jo  Jo 

=  -a  f  D2v2D2ui  dx  +  a  [  u2D4V\  dx  -  /3Du2D2v2  |i  +p  f  D3v2Du2dx 
Jo  Jo  Jo 

=  -a  I  D2v2D2u\  dx  +  a  I  u2D4V\  dx  +  /3u2Dzv2  In  — /?  f  D4v2u2  dx 
Jo  Jo  Jo 


0 

aD4 


-1 

-pD4 


\ 

\ 

1>1 

V2 ) 

J 

(66) 


/ 


Thus  (— A3)*  = 


V 


0  -1 
aD4  - PD 4 

Notice  that  this  operator  is  not  skew-adjoint.  Hence,  the  operator  —A?,  is  not  the 
generator  of  a  group  as  in  the  previous  cases. 


It  is  necessary  to  check  whether  the  operator  (-A3)*  is  dissipative.  The  details  here 
follow  as  for  -^3  except  for  a  few  sign  changes  in  the  intermediate  steps  so  that  the  result 
is  the  same,  and  consequently,  the  details  are  omitted. 
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It  is  also  important  that  -A3  be  a  closed  operator.  But  the  adjoint  of  a  densely 
defined  linear  operator  is  always  closed  (see  Theorem  25).  Since,  in  the  present  case 
,  —A3  =  ((-.A3)*)*  and  (-.A3)*  is  densely  defined  (same  domain  as  .A3),  then  -A3  is  closed 
as  desired. 

The  existence  of  a  unique  Cl  solution,  u(t),  to  (60)  for  any  uq  6  D(Az)  is  now 
guaranteed  by  Theorem  27  and  Theorem  29. 


2.6  Spatial  dependence  in  the  damping  coefficient 

Consider  utt  +  D2((3(x)D2ut)  +  aD4u  =  0.  Assume  ffix)  >  0  for  all  x  e  [0, 1],  and  a 
a  positive  constant  as  before.  Formulate  the  abstract  system  as  before  with  the  operator 
defined  as  follows: 


A  =  A4  = 


/ 


V 


0  -1  | 

aD4  D2(/3(-)D2)  j 


(67) 


Note  that  D(A4)  =  D{Az).  It  will  now  be  established  that  -A4  is  dissipative.  The 
details  follow  the  pattern  of  the  case  for  constant  /?.  Boundary  terms  are  zero  as  before. 
The  symbol '  and  the  symbol  D  will  both  be  used  to  represent  differentiation  with  respect 
to  the  spatial  variable.  Let  u  e  -D(A4)  be  given.  Then, 


/ 

0  1  \ 

(  \ 
tii  ' 

(-A4ti,ti)  = 

\ 

\  -aD4  -D2(p(-)D2)  j 

\U2J 

\U*  )) 

u2 


^  -aD4Ul  -  D2(p(-)D2u 2)  j 


ui 


U2 


w 


a  j  j,2"ui"dx  +  J  (~aD4ui  -  D2(/3(x)D2  ti2))  «2  dx 

aujui"  |o  -a  f  U2Ui'"dx-au2Ui"\l+a[  uz'ui"  dx 
Jo  jo 

-u2D(P(x)D2u2)  |o  +  J  U2  D((3(x)D2u2)  dx 

U2(3(x)D2U2  lo  -  /  U2"P(x)U2"  dx 
Jo 

-  /  P{x)(u2")2  dx 
Jo 


(68) 
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and  this  is  clearly  not  positive  since  P(x)  >  0  for  each  x. 

Next  the  adjoint  (-A4)*  of -A4  will  be  computed.  Again  the  details  follow  as  in  the 
case  for  constant  /?.  For  u  €  D(A4), 


(~A4u,v)  = 


// 


0  1 
- aD 4  -D2(p(-)D2) 

«2 


\ 

(  \ 

(  \\ 

Ul  1 

J 

\U2) 
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[V2  j) 

\ 

/  ,  \\ 

f  *>1 

/ 

\V2J> 

-aD4U!  -  D2(p(-)D2u2) 

=  a  J  u2"vi"  dx  —  J  (aD4u\  +  D2(P(x)D2  «2  fj  v2  dx 

=  OiU2v\"  lo  -a  f  U2Vi"dx-av2Ui"\l+a[  v2u\ dx 
Jo  Jo 

-v2D(P(x)D2u2)  |5+/o  V2  D(P(x)D2u2)  dx 
=  -au2vi"  |o  +a  J  %2D4V\  dx  +  av2'u\"  |J  -a  J  v^'ux"  dx 

+v2'P(x)u2"\l-  [  vJ'PWul'dx 
Jo 

=  -a  [  v2"ui"  dx  +  a  [  u2J94vi  dx  -  v2"P(x)u2"  |o  +  f  u2  D(P(x)v2")  dx 
Jo  Jo  Jo 

=  -a  f  v2"ui"dx  +  a[  u2D4v\  dx  +  U2D(P(x)v2")  |o  -  [  U2D2(p{x)v2)dx 
Jo  Jo  Jo 


U\ 

«2 


0  -1 
aD4  -D2(P(-)D2) 


\ 

(  X\ 

f  Vi 

) 

l  V2)j 

=  (u,(-A4)*v). 


Therefoie, 


(-A4r  = 


0  -1 

aD4  -D2(P(-)D2) 


(69) 


(70) 


The  argument  to  show  that  (- A4)*  is  dissipative  follows  the  established  pattern  and 
has  the  desired  result. 


38 


The  domain,  ,D((— A4)*),  of  (— A4)*  is  the  same  as  D(A^).  This  is  clear  from  an 
inspection  of  the  operators.  Hence  (— A4)*  is  densely  defined.  Note  that  ((-Aj)*)*  =  —  A4 
and  thus  -A\  is  closed  by  Theorem  25. 

Existence  of  a  Cl  solution  in  now  assured  by  Theorem  27  and  Theorem  29  as  before. 
2. 7  Chapter  summary 

Basic  theorems  of  semigroup  theory  have  been  reviewed  and  the  concept  of  abstract 
formulation  of  a  differential  equation  has  been  discussed.  A  careful  description  of  a  linear 
space  and  a  careful  selection  of  operators  have  been  provided  to  demonstrate  the  termi¬ 
nology.  The  hypotheses  of  an  appropriate  existence  theorem  were  shown  to  be  satisfied. 

The  beam  equations  for  constant  and  spatially  varying  a  have  been  considered.  Also, 
for  constant  a,  the  cases  of  constant  and  spatially  varying  (3  have  been  considered.  Exis¬ 
tence  and  uniqueness  of  solutions  has  been  established. 

In  the  next  chapter,  consideration  is  given  to  a  nonautonomous  problem.  All  of 
the  present  chapter’s  concepts  will  be  needed  there.  Additional  concepts  will  also  be 
introduced. 
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III.  The  nonautonomous  problem 


In  this  chapter  the  coefficient  of  damping  is  allowed  to  vary  with  time.  A  suitable 
theory  for  this  case  will  be  presented.  Combined  temporal  and  spatial  dependence  is  not 
difficult  once  the  time  dependent  case  is  complete.  The  case  of  a  time  dependent  a  is  also 
treated.  This  will  require  the  introduction  of  an  additional  theorem. 

S.l  An  overview  of  the  simple  nonautonomous  case 

Consider  utt+D'2(0(t)D2Ut)+aD4u  =  0.  Assume  a  >  0and/3(i)  >  0  for  all  t  6  [0 ,T]. 
Furthermore,  /?  is  assumed  to  be  continuously  differentiable.  As  before,  formulate  the 
problem  as  the  abstract  system 


ut  +  Au  =  0;  ti(0)  =  u0 


(71) 


where  now, 


A  =  A(i)  =  As  = 


/ 

V 


0  -1 
aD4  (3(t)D4 


\ 

/ 


(72) 


where  -D(As)  =  D(Az).  Note  that  for  each  t  this  behaves  the  same  as  in  the  case  of 
constant  (3.  The  arguments  for  the  adjoint  and  dissipativity  are  not  repeated. 


Because  of  the  explicit  time  dependence,  (71)  is  called  an  evolution  equation.  For 
each  t  >  0,  the  previous  arguments  apply  to  establish  that  -As  is  the  generator  of  a  Co 
semigroup  of  contractions  which  provides  a  solution  as  before.  If  /3(t j  is  suitably  smooth 
it  is  reasonable  to  expect  to  be  able  to  piece  together,  from  the  solutions  for  individual 
values  of  t,  an  overall  solution.  However,  the  theorems  which  gave  solutions  for  individual 
values  of  t  do  not  guarantee  the  necessary  behavior  (eg  continuity,  differentiability)  for 
their  composite  to  be  a  solution  of  the  evolution  equation. 


The  solution  strategy  is  as  follows.  At  t  =  0  the  state  of  the  system  is  specified. 
The  solution  of  the  equation  for  /?(0)  is  some. surface,  as  in  Figure  2.  But,  since  (3  changes 
with  time,  this  solution  is  only  accurate  for  small  values  of  t.  Suppose  T,  the  largest  value 
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u 


Figure  2.  The  initial  solution 

of  interest,  is  larger  than  the  interval  for  which  the  solution  based  on  /?( 0)  is  sufficiently 
accurate.  Then  an  iterative  procedure  is  used.  This  will  be  briefly  outlined. 

Let  0  =  t0  <  <  t2  <  •••  <  tn  =  T  describe  a  partition  of  [O,!1].  Apply  the 

results  of  the  earlier  work  to  the  constant  coefficient  case  corresponding  to  j9(0).  Let 
u(t)  =  So(t  -  fo)«o  on  [*0)*i]>  where  Sq  is  the  semigroup  generated  by  -As(0).  Then 
u(ti)  is  an  estimate  of  u(ti).  Apply  the  earlier  work  again,  this  time  for  the  constant 
coefficient  P{i\).  This  gives  u(t)  =  S\(t  -ti)u(ti)  on  Take  £(*2)  as  an  approximate 

value  for  u(t2).  If  fi(t)  is  appropriately  smooth  then  a  solution  with  any  specified  accuracy 
is  possible  by  choosing  a  small  enough  upper  bound  on  max{tj  —  tj-\\  j  =  1,2, 
Conceptually,  this  is  like  Euler’s  method  in  numerical  analysis. 

In  this  section,  necessary  conditions  for  the  solution  strategy  to  make  sense  are 
identified.  This  will  be  preceeded  by  appropriate  introductions  of  additional  terminology 
and  theorems. 

In  the  earlier  solution  strategy  a  semigroup  S(t)  was  obtained  from  its  infinitesimal 
generator  -A3.  But  A5  =  A(t)  depends  on  t  and  there  is  some  risk  of  confusion  of 
parameters.  The  semigroup  generated  by  —A(t)  will  be  denoted  5j(s).  Suppose  now  that 
a  limiting  operator,  which  always  uses  the  current  element  of  the  semigroup,  is  desired. 
This  might  be  represented  by  something  like  Ss(s),  but  this  would  certainly  be  confusing. 
An  operator  is  needed  which  can  propogate  from  some  time  s  to  some  time  t  while  using 
the  appropriate  element  of  the  semigroup  at  each  instant  of  time. 

The  common  choice  of  notation  for  this  situation  is  U  with  two  parameters.  Namely, 
U(t,s )  is  used  for  the  desired  operator,  when  it  exists.  The  term  evolution  system  is  used 
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to  refer  to  U(t,s).  (Some  authors  use  the  term  evolution  operator.)  So,  while  U  and  S  play 
roles  which  are  conceptually  very  similar,  the  terminology  is  quite  different.  The  following 
standard  definition  is  taken  from  [58:pg  129]. 

Definition  33  A  two  parameter  family  of  bounded  linear  operators  U(t,s ),  0  <  s  <t  <T, 
on  X  is  called  an  evolution  system  if  the  following  two  conditions  ore  satisfied: 

1.  U(s,s )  =  I,  U(t,r)U(r,s)  =  U(t,s)  for  0  <s<r<t<T. 

2.  ( t,s )  ->  U(t,s )  is  strongly  continuous  for  0  <  s  <  t  <  T. 

When  there  is  an  evolution  system  corresponding  to  a  problem  with  a  time  dependent 
operator,  the  evolution  system  is  used  to  produce  a  solution  in  the  same  way  a  semigroup 
is  used  for  the  autonomous  case.  That  is,  u(t)  =  U(t,0)uo. 

Notice  that  property  1  of  an  evolution  system  (Definition  33)  says  that  small  steps 
and  large  steps  give  the  same  result  as  long  as  the  ultimate  end  points  are  the  same. 

With  As  allowed  to  depend  on  t,  the  same  basic  ideas  as  before  are  still  applicable.  It 
is  necessary,  however,  to  modify  the  requirements  on  A$  to  assure  existence  of  a  solution. 
Of  course  some  smoothness  of  the  map  t  -*  A(t)  is  required.  The  map  is  required  to  be 
continuous  and  also  the  concept  of  a  stable  family  {A(<)}(6[oj]  is  introduced.  Further, 
since  A5  is  allowed  to  change  as  t  changes,  there  is  the  possibility,  for  an  iterative  scheme, 
that  the  image  of  v4s(ti)  would  fail  to  be  in  the  domain  of  Asfa)  (see  Figure  3).  This 
would  cause  the  iterative  scheme  to  fail.  It  is  useful  to  identify  an  appropriate  subspace  Y, 
Y  C  X,  such  that  for  all  t  G  [0,T],  it  happens  that  Y  C  D(A(t)).  In  this  regard  the  notion 
of  7  being  A(t)-admissible  is  presented.  Such  a  Y  will  not  be  very  meaningful  unless  it  is 
dense  in  X.  This  will  be  required. 

When  a  solution  is  obtained  in  Y,  it  is  termed  a  Y-valued  solution  [58:pp  139-.140]. 
This  is  not  a  significant  restriction  for  the  application  under  consideration  since  Y  will  be 
all  of  D(A).  However,  the  properties  of  Y-valued  solutions  are  useful. 

Definition  34  A  function  u  €  C([0,T]  :  Y)  is  a  Y-valued  solution  of  the  initial  value 
problem  (71)  if  u  G  C1((0,CT] :  X)  and  (71)  is  satisfied  in  X. 
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Figure  3.  The  different  domains 

Theory  sufficient  to  establish  the  existence  of  a  unique  solution  is  presented  in  the 
next  section.  The  development  is  gradual  in  that  several  of  the  theorems  will  be  minor 
extensions  of  the  immediately  preceeding  theorem.  While  this  means  there  will  be  some 
overlap,  it  also  has  the  advantage  of  step  by  step  development.  Most  of  the  proofs  will  be 
by  citation  only.  That  is,  while  all  the  theorems  are  available,  is  seems  prudent  to  gather 
statements  of  the  theorems  together  here.  It  is  not  necessary  to  reproduce  all  the  proofs. 

3.2  The  simple  nonautonomous  case 

The  specific  problem  introduced  in  (71)  and  (72)  is  now  addressed. 

3.2.1  Some  technical  preliminaries 

A  fundamental  concept  is  that  of  a  stable  family  [58:pg  130]. 

Definition  35  A  family  {yl(<)}<e[0,T]  °f  infinitesimal  generators  of  Co  semigroups  on  a 
Banach  space  X  is  called  stable  if  there  are  constants  M  >  1  and  u  (called  stability  con¬ 
stants)  such  that 

p(A(t ))  D  (w.  oo)  for  all  t  6  [0.T]  (73) 

and 

k 

n  i2(A  :  A(tj))  <  M( A  -  u})~k  (74) 

5= i 

for  A  >  u  and  every  finite  sequence  0  <  t\  <  <2  <  •  •  *  <  tk  fs  T]  k  =  1, 2,  •  •  •. 
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For  k  =  1  this  is  just  the  resolvent  condition  for  an  individual  operator  A  to  be  the 
generator  of  a  Cq  semigroup  [58:pg  12,  Corollary  1.3.8],  Stability  of  the  family  essentially 
means  that  the  bounds  on  the  resolvents  of  individual  operators  A(tj)  can  be  combined 
to  give  a  bound  on  the  composition  of  resolvents  of  j4(tj)’s.  There  are  other  conditions 
equivalent  to  (74),  eg  [58:pg  131,  Theorem  5.2.2]. 

The  next  theorem  is  a  standard,  very  powerful,  result. 

Theorem  36  (Hille-Yosida)  A  linear  operator  A  is  the  infinitesimal  generator  of  a 
Co  semigroup  of  contractions  S(t),  t  >  0,  if  and  only  if  A  is  closed,  D(A)  =  X,  and 
the  resolvent  set  p(A)  of  A  contains  &+,  and  for  every  A  >  0 

||1£(A  :  .A)||  <  1/A.  (75) 

Proof:  See  [58:pp  8-11],  □ 

There  will  be  several  occasions  to  use  the  following  criterion  to  establish  that  a 
particular  family  is  stable. 

Lemma  37  If,  for  each  t  €  [0,T],  Aft)  generates  a  Cq  semigroup  of  contractions,  then 
the  family  is  stable. 

Proof:  This  is  a  straightforward  application  of  Definition  35.  For  contraction  semigroups 
it  is  appropriate  to  choose  w  =-  l .  Then,  for  each  tj,  appeal  to  Theorem  36  to  bound  the 
resolvent  by  1/A.  Now  it  is  clearly  suitable  to  choose  M  =  1.  □ 

Now  the  phrase  A-admissible  will  be  defined  [58:pg  122]. 

Definition  38  Let  S(t)  be  a  Co  semigroup  and  let  A  be  its  infinitesimal  generator.  A 
subspace  Y  of  X  is  called  A-admissible  if  it  is  an  invariant  subspace  of  S(t),  t>  0,  and  the 
restriction  of  S(t)  to  Y  is  a  C0  semigroup  in  Y  (ie,  it  is  strongly  continuous  in  the  norm 

II  •  II y). 
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3.2.2  Key  theorems 

A  first  step  toward  finding  the  evolution  system  for  (71)  is  given  in  the  following 
theorem  [58:pg  135].  The  space  Y  in  the  theorem  has  not  yet  been  chosen.  When  the 
choice  is  made,  it  will  be  with  these  hypotheses  in  mind.  The  points  of  Y  will  come  from 
X  and  the  restriction,  of  an  operator  A,  to  Y  will  be  called  the  part  of  A  in  Y  and  be 
denoted  by  A.  When  it  is  time  to  apply  these  rems,  the  place  of  A  will  be  taken  by 
-A5  defined  in  (72). 

Theorem  39  Let  A(t),  0  <  t  <  T,  be  the  infinitesimal  generator  of  a  Co  semigroup  St(s), 
s  >  0,  on  X.  Suppose  the  family  {A(t)}i€[o(T]  satisfies  the  conditions 

1.  {A(t)}tg[0,T]is  a  stable  family  with  stability  constants  M,u. 

2.  Y  is  A(t)-admissible  for  t  6  [0,T]  and  the  family  {A(t)}t6r0|T)  of  parts  A(t)  of 
A(t)  in  Y,  is  a  stable  family  in  Y  with  stability  constants  M,w. 

3.  For  t  €  [0,  T\,  D(A(t))  D  Y,  A(t)  is  a  bounded  operator  from  Y  into  X  and  t  -*•  A(t) 
is  continuous  in  the  B(Y,X )  norm  ||  •  jjy_+^. 

Then,  there  exists  a  unique  evolution  system  U(t,s),0  <  s  <  t  <T,  in  X  satisfying 

1.  ||  17(4,  s)||  <  Mcw(‘-S)  for  0  <  s  <  t  <  T. 

2.  7ftU(t,s)v  |<=s=  A(s)u  for  t;  G  Y,  0  <s<T. 

3.  -§^U(t,s)v  =  —U(t,s)A(s)v  for  v  G  Y,  0  <  s  <  t  <  T. 

The  derivative  from  the  right  in  the  second  item  and  the  derivative  in  the  third  are  in  the 
strong  sense  in  X. 

Proof:  See  [58:pg  135-138],  □ 

The  properties  of  U  established  in  the  theorem  are  useful  in  demonstrating  the 
uniqueness  of  the  candidate  solution  it  generates.  The  portion  of  the  proof  which  describes 
the  construction  of  the  evolution  system  is  repeated  here  for  convenience  and  because  of 
its  role  in  the  construction  of  a  solution. 
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Consider  an  approximation  of  the  family  by  piecewise  constant  families 

{J4n(t)}te[0,T])  n=  1)2,..*,  defined  as  follows:  Let  =  ( kjn)T ,  k  =  0, 1, . . . ,  n  and  let 


An(t)  =  A(t%)  for  fg  <  t  <  t£+1,  k  =  0,l,...,n-  1 
^  An(T)  =  A(T). 

Since  t  -*  A(t)  is  continuous  in  the  B(Y,X )  norm  it  follows  that 

||A(t)  -  yfn(t)||y-Oi:  ->  0  as  n  -»  oo 


(76) 


(77) 


uniformly  in  t  G  [0,T].  From  the  definition  of  An(t)  and  the  hypotheses  of  Theorem  39  it 
follows  readily  that  for  n  >  1,  {^n(0}i6(o,T]  is  a  stable  family  in  X  with  constants  M,u> 
while  {-4n(0}t€[o,T]  is  a  stable  family  in  Y  with  constants  M,u. 

Next,  for  each  n  define  a  two  parameter  family  of  operators  Un(t,  s ),  0  <  s  <  t  <  T 

by, 

Stn(t  -  s)  for  <  s  <  t  <  Vj+1 

Un(t,  s)  =  •  a»((  -  (J)  [nf;/+1  sf,  (I)]  s,r  W+1  - »)  (78) 

for  k  >  l,  <  t  <  t%+1 ,  t "  <  s  <  . 

It  is  straightforward  to  verify  that  Un(t,  s )  is  an  evolution  system.  Then,  let 


U(t,s)x  —  lim  Un(t,s)x  for  x  e  X,  0  <  s  <  t  <  T. 


(79) 


Details  of  the  proof  that  this  limit  exists  and  is  an  evolution  system  are  in  [58:pg  135-138]. 
Now  consider  the  associated  uniqueness  result. 


Theorem  40  Let  {-A(t)}46[0,T]  be  a  family  of  infinitesimal  generators  of  Co  semigroups 
on  X  satisfying  the  conditions  of  Theorem  39.  If  the  initial  value  problem  (71)  has  a 
Y-valued  solution  u,  then  this  solution  is  unique;  and  moreover,  it  is  given  by 


where  U(t,s)  is  the  evolution  system  provided  by  Theorem  39  and  v  is  the  value  of  v  at 
time  t  =  s. 

Proof:  See  [5S:pg  140].  □ 

The  next  theorem  completes  the  list  of  properties  required  of  U  in  older  to  guaiautec 
existence  of  a  solution.  In  its  proof,  the  need  for  results  so  far  obtained  becomes  clcai. 

Theorem  41  Let  {>l(0}:e[o,Tj  satisfy  the  conditions  of  Theorem  39  and  lei  U(t,$),  0  < 
s  <t  <T  be  the  evolution  system  given  in  Theorem  39.  If 

1.  U(i,s)Y  C  Y  for  0  <  s  <  t  <  T  and 

2.  For  v  e  U(t,s) v  is  continuous  in  Y  for  0  <  $  <  t  <  T 

then  for  every  v  €  V,  U(t,s)v  is  the  unique  Y -valued  solution  of  the  initial  value  problem 

du(t)/dt  =  A(t)u(t)  for  0  <  s  <t  <T  (SI ) 

w(s)  =  v.  (82) 

Proof:  See  (53:pg  141].  □ 

3.2.3  A  strengthened  hypothesis 

An  alternative  for  the  second  hypothesis  of  Theorem  39  is  how  presented.  It  will 
be  referred  to  as  condition  2+.  This  condition  appears  in  [35:pg  13S].  It  will  be  used  to 
establish  the  additional  hypotheses  of  Theorem  41. 

(2+)  There  is  a  family  {Q(0}te[o,T]  of  isomorphisms  ([54:pg  173]  or  [65:ppl$4,  199]) 
of  I'  onto  A'  such  that  for  every  v  £  Y.  Q{t)v  is  Lipschitz  continuous  in  .V  on  /  €  [0.T] 
and 

=  A(l)  +  lJ(i)  (33) 

where  B(i),  0  <  t  <  T,  is  a  strongly  continuous  family  of  bounded  operators  on  X. 

The  next  lemma  establishes  that  this  is  indeed  at  least  as  strong  a  condition  as 
hypothesis  2  of  Theorem  39. 
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Lemma  v  Hypothesis  1  of  Theorem  89  and  condition  2+  imply  hypothesis  2  of  Theorem 
89. 

Proof:  See  [58:pg  142].  The  condition  2+  is  not  as  strong  as  the  corresponding  condition 
in  the  reference.  Yet,  the  proof  carries  over  with  no  essential  modification.  Also,  see 
comments  in  Appendix  D.  □ 

Next  it  will  be  shown  that  the  condition  2+  is,  in  fact,  stronger  than  hypothesis  2  of 
Theorem  39.  A  preliminary  lemma  is  appropriate. 

Lemma  43  Let  U(t,s),  0  <  $  <  t  <  T  be  an  evolution  system  in  a  Banach  space  X 
satisfying  ||Z7(t,s)||  <  M  for  0  <  s  <  t  <  T.  If  H(t)  is  a  family  of  integrable  linear 
operators  in  X  such  that  for  almost  every  t,  ||/f  (t)\\  <  H  <  oo,  then  there  exists  a  unique 
family  of  bounded  linear  operators  V(t,s),  0  <  s  <t  <T  on  X  such  that 

V(t,s)x  =  U(t,s)x  +  J  V(t,r)H(r)U(r,s)xdr  for  x£X  (84) 

and  V(t,s)x  is  continuous  in  s,tfor  0  <  s  <t  <  T. 

Proof:  The  proof  is  standard  for  Volterra  integrals  of  the  second  kind.  For  details  see 
[58:pp  142-143].  The  hypothesis  of  integrability  is  replaced,  in  the  reference,  with  strong 
continuity.  But,  with  the  bound  II  described  in  the  hypothesis,  the  cited  proof  holds  with 
no  essential  modification.  □ 

Theorem  44  Let  A(t),  for  0  <  t  <  T,  be  the  infinitesimal  generator  of  a  Co  semigroup  on 
X.  If  the  family  {A(t)}jer0ir]  satisfies  the  conditions  1  and  8  of  Theorem  39  and  condition 
2+,  then  there  exists  a  unique  evolution  system  U(t,s ),  0  <s<t<T,inX  satisfying  the 
following  5  conditions: 

1.  ||tf(t,s)||  <  Me^~s)  for  0  <  s  <  t  <  T, 

2.  jj-U(t,s)v  | t=s=  A(s)v  for  v  €  Y,  0  <  s  <  T, 

3.  ys U(t,s)v  =  -U(t,s)A(s)v  for  v  €  Y,  0  <  s  <  t  <  T, 

U(t, s)Y  C  Y  for  0  <  s  <  t  <  T,  and 
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5.  For  v  €  Y,  U(t,s)v  is  continuous  in  7  for  0  <  s  <  t  <  T. 

Proof:  See  [58:pp  143-145].  Some  comment  on  the  cited  proof  is  appropriate  since  the 
condition  2+  in  this  work  is  weaker  than  the  corresponding  condition  in  the  reference. 
In  the  reference,  Q  is  required  to  be  continuously  differentiable.  Then  Q  represents  the 
derivative.  In  this  work,  Q  is  only  required  to  be  Lipschitz  continuous.  But  Lipschitz 
continuity  implies  absolute  continuity  which  implies  that  there  is  some  integrable  function, 
which  will  also  be  denoted  Q,  such  that 

Q(t)  =  Q(  0)+  fQ{r)dr.  (85) 

Jo 

A  similar  result  is  obtained  in  [38:pp  505-507], 

It  follows  that  Q  is  differentiable  a.e.  with  the  derivative  given  by  Q  wherever  the 
derivative  exists. 

The  boundedness  of  Q~l  and  the  differentiability  of  Q  a.e.  lead  to  the  differentiability 
of  Q -1  a.e..  Where  it  exists,  the  derivative  is  given  by 

^  (QW1*)  =  -Q{t)~lQ(t)Q(t)-lx.  (86) 

With  these  observations,  the  proof  goes  through  as  cited.  □ 

These  results  can  be  combined  as  a  corollary,  as  in  [58:pg  145]. 

Corollary  45  Let  {A(t)}t€[o,r]  a  family  of  infinitesimal  generators  of  Co  semigroups 
on  X .  If  {A(f)}ie[o,T]  satisfies  the  hypotheses  of  Theorem  44  far  every  v  €  Y  then  the 
initial  value  problem 

du(t)/dt  =  A(t)u(t)  for  s  <  t  <  T,  (87) 

u(s)  =  v,  (88) 

has  a  unique  Y -valued  solution,  u,  on  s  <t  <T. 

3.2.4  The  existence  theorem  for  solving  (71) 
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The  problem  of  solving  (71)  will  fall  under  the  following  special  case  of  Corollary 
45.  The  theorem  is  to  be  applied  to  -As,  as  discussed  previously.  The  next  theorem  is  a 
special  case  of  Theorem  5.4.8  in  [58:pg  145].  See  also  [31:pg  252]. 

Theorem  46  Let  {-4(t)}(g(o,Tj  be  a  stable  family  of  infinitesimal  generators  of  Co  semigroups 
on  X.  If  D(A(t))  =  D  is  independent  of  t  and  for  v  G  D,  A(t)v  is  Lipschitz  continuous  in 
X  then  there  exists  a  unique  evolution  system  U(t,s),  0  <  $  <  t  <  T,  satisfying  the  five  re¬ 
sults  of  Theorem  44  where  Y  is  the  set  D  equipped  with  the  norm  [|u||y  =  ||u||x  +  ||A(0)u||jf . 

Proof:  See  [58:pp  145-146].  Also  see  Appendix  D  which  elaborates  on  one  portion  of  the 
proof  cited.  E 

It  may  be  noted  that  Theorem  46  made  no  mention  of  condition  2+.  The  hypotheses 
of  Theorem  46  are  such  that  Q(t )  =  I  +  A(t)  is  an  acceptable  choice.  Thus,  whenever 
Theorem  46  applies,  there  is  no  need  to  explicitly  identify  the  isomorphism  Q.  In  the 
following  application,  the  choice  Q  =  I  -  (-As)  =  I  +  As  will  be  appropriate. 

8.2.5  An  application  of  Theorem  46  to  the  case  of  nonautonomous  damping 

The  hypotheses  of  Theorem  46  are  satisfied  for  the  problem  represented  by  (71)  and 
(72).  Theorem  46  is  to  be  applied  to  -As-  The  satisfaction  of  the  hypotheses  will  be 
considered  in  some  detail.  Note  that  D(As )  =  D(Az). 

Let  7  be  the  same  set  of  points  as  D(As)  but  with  the  graph  norm  of  ^5(0) 

II*  =  WU  +  IM5(0)*.  (89) 

Lemma  47  The  linear  space  Y  is  a  Banach  space  under  ||  •  ||y. 


Proof:  Since  it  is  easy  to  see  that  Y  is  a  normed  linear  space,  all  but  completeness  is 
clear.  Let  { yn }  be  a  Cauchy  sequence  in  Y.  It  is  immediate  that  {y„}  and  {As(0)j/n}  are 
Cauchy  sequences  in  X.  Hence,  there  exist  y,  y  6  X  such  that  yn  —>  y  and  As(0)yn  — ^  V- 


Since  -  As(0)  is  the  generator  of  a  Co  semigroup  of  contractions  (as  shown  in  the  case 
of  /?  a  constant)  it  follows  that  A5(0)  is  a  closed  operator  (see  Theorem  36).  (The  closedness 
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of  --As(O),  as  an  operator  in  X ,  was  established  in  showing  that  it  is  the  generator  of  a 
Co  semigroup.)  From  the  closedness  of  -As(O)  it  follows  that  y  £  D(A$)>  and  hence  y  £  Y. 
Thus  Y  is  complete.  □ 

For  each  t  £  [O.T]  it  is  clear  that  -As(t)  is  the  generator  of  a  Co  semigroup  of 
contractions  (from  the  case  for  constant  /?).  Hence,  by  Lemma  37,  the  family  is  stable  as 
required.  The  domain  has  no  time  dependence  as  is  easily  seen  from  its  definition. 

Finally,  it  is  necessary  to  verify  the  Lipschitz  continuity  of  A5  with  respect  to  t.  Let 
v  £  D  be  given.  First  consider  a  straightforward  approach  to  continuity.  Let  e  >  0  be 
given.  Identify  6  >  0  such  that 

|  -  <2  |<  6  =>  ||As(ii)t>  -  A5(t2)t>|U  <  e.  (90) 

From  the  definitions 

(  0  0  \(Vl\ 

V  0  (/?(<!)  -  /?(i2))X>4  )  \V2  ) 

0 

mi)-m))D4v2 

f  rl  \ 

=  \m-m\  \}ADiv2)Hx)  (91) 

Now,  V2  £  HA  and  /?  £  C°  makes  clear  the  existence  of  suitable  6  >  0.  Specifically, 
v2  £  H4  =>•  3 M  such  that  ||u2||tf4  <  M.  Also,  continuity  of  (3  guarantees  that  there  is 
some  6  such  that 

\ti-h\<8*\0(ti)-P(&)\<e/M.  (92) 

This  is  the  required  6. 
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In  fact,  whatever  smoothness  ft  has  will  carry  over  to  As(t)v.  For  example,  to  identify 
$-tAs(t)v,  consider 


lim  =  lim 

/i— >0  h  /i— t0 


=  lim 
h-*0 


l/o  0 

^  1  0  (/ 3(t  +  h)-(3(t))D 4 


\ 


/  \ 

\ 

/ 

\ 


/ 


[(P(t  +  h)-p(i))D1 2 3 4V2  jj 

\ 


mi^m D4V2 


'  o  N 

P'D\  2 


(93) 


Since  ft1  is  continuous  and  D4v 2  does  not  depend  on  t  it  is  now  clear  that  A5(t)t>  is  contin¬ 
uously  differentiable  with  respect  to  t  for  t  €  (0,  T).  Thus,  in  this  particular  application,  a 
stronger  condition  than  Lipschitz  continuity  is  satisfied. 


This  completes  the  verification  of  hypotheses  for  the  application  of  Theorem  46  to 

(71). 

A  brief  review  of  what  is  known  about  solving  (71),  as  a  result  of  satisfying  the 
hypotheses  of  Theorem  46,  is  appropriate. 


1.  Theorem  46  guarantees  the  existence  of  U(t,s)  satisfying  the  five  results  of  Theo¬ 
rem  44. 

2.  Corollary  45  says  the  problem  has  a  unique  7-valued  solution. 

3.  Theorem  41  says  I7(f,0)uo  is  the  unique  7-valued  solution  of  the  equation. 

4.  Theorem  39  (proof)  describes  the  construction  of  a  sequence  { Un }  whose  limit  is  U. 

Since  U(t,s)  is  bounded,  the  continuous  dependence  on  the  initial  condition  vector 
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follows  as  in  Theorem  31. 


8.8  Combined  t,x  dependence  in  the  damping  coefficient 

Consider  utt  +  D2(P(t,x)D2ui)  +  a D4u  =  0.  Require  the  constant  a  >  0  and 
j 6(t,x)  >  0  as  before.  Furthermore,  require  P  and  to  be  twice  continuously  differ¬ 
entiable  with  respect  to  a:,  and  continuously  differentiable  with  respect  to  t.  The  problem 
can  be  formulated  as  a  system  as  in  (71)  with 


0  -1  \ 
aD 4  D\p(t,-)D2)  j 


(94) 


Note  that  D(Ae)  =  jD^). 

It  is  not  difficult  to  see  that  the  hypotheses  of  Theorem  46  are  again  satisfied.  Indeed, 
for  each  t,  —  A<$  generates  a  Co  semigroup  of  contractions  as  in  the  P(x)  case.  Hence,  the 
family  {-i46(t)}te[o)rj  is  stable.  Furthermore,  note  that  D(-Aq)  does  not  depend  on  t. 
Also,  -A&(t,  •)  £  C1  as  long  as  P  is  appropriately  smooth,  as  is  required.  Some  detail  is 
appropriate  for  this  last  point. 

The  computation  which  parallels  (91)  is  now  presented. 


Me(<i)t>-  AG(t2)v\\x  = 


0  0 
0  D2(P(tux)-p(t2,x))D2 

0 

D2  {{P(h,x)-P(t2,x)\D2v2) 


vi 

V2 


X 


=  (jf  \d2  {[P(i  1 , x)  - P(t2 , x))D2v^j ] 2  dx'j 
=  Qf  \D  ([P(.h,x)~  P(t2,x)}D3v2 

/  Q  Q  \  \-|2  \*/2 

+  [dift*1  ’  “  dx *7 1,2 V2J  dx) 

=  [(P{h,x)~  p(t2,x))D4v2 

+2  (J-p(h,x)-^t2,x^D% 


53 


\ 

2  \ 

j  -D21>2 

dx 

1/2 


<  31/2  (jf  [(/?(*i,*)-/?(i2,*))2  {d*v2) 

+4  (j^  (/5(<i.*) "  P(h,x))j  (f>3v2) 

+  (j^Wi,x)-/3(h,x))j  (D\) 

<  31/2  [Q  [/?(*i,x)  -  /?(t2>x)]2  (-E>S)2  dx^ 

+2  (/?(*i,x)  -  /3(*2,x))]  (-D3^)2  dx 

+  (jo  [£»(^(<1’*)‘"^(<2’*))  l1^2)2 


\  1/2 

da;  1 
1/2 


1/2 

1/2' 


(95) 


The  key  to  continue  is  that  ^  is  twice  continuously  differentiable  with  respect  to  x  on  a 
closed  and  bounded  set.  Thus,  there  exists  some  Mp  such  that 


sup 

(*,z)€[0,T]x(0,l} 


(JW  s  M > 


(d  (  8*  > 
SUP  o7  a~2P 

(t,r)e(o,T]x[o,i]  ) 


<  Mp. 


(96) 


Then, 


P(h,x)~  P(t2,x)  <  Mp\ti-t2\ 
■g^Wi,x)-ff(ta,x))  0 


Txe(.h,x)  -  Tjoa,x) 


< 


fd  (  d  r\\.  . 

~‘2' 

<  Mp\h-h\ 


d 2 


,*)-/?(<2,*))  < 


(97) 
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Then,  to  continue, 


peMv-AjteMU  <  31/2  Mp\ti  -  t2\  (J  (d4v2)2  da) 

+2Mfi\ti  -  h |  Qf  (£>3u2)2  da) 

+Mfi\ti  -h\  Qf  (d2v2)2  da)  j 
r  i  2/2 

<  31/2AfJJ|t1  -  i2j  (7)4u2)2  da)  ' 

+2  ( Dzv2 )2  dx^  +  Qf*  (p2vtf  da)  '  .  (98) 

From  this  point  the  argument  which  follows  (91)  applies  as  before. 

The  computation  which  parallels  equation  (93)  is  presented  next. 


lim  -46(t  +  h,x)v~  Ae(t,x) 

h — >0  }i 


lim 
h-+  0 


0  0 
0  D7(p(t  +  h,x)-P(t,x))D 2 


=  lim  I  - 


D2  {{P{t  +  h,x)~  0(t,  a))  D2v2) 


Continuity  with  respect  to  t  follows  from  an  inequality  similar  to  (98). 

Thus,  Theorem  46  applies  to  give  the  evolution  system  in  terms  of  which  a  solution 
to  the  problem  is  given  as  described  in  the  P(t)  case. 
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3,4  A  more  general  theorem 

The  most  general  theorem  used  in  this  work  is  Theorem  48.  It  is  a  specialized 
combination  of  Theorem  3.1  of  [35:pp  169-170]  and  The.  rem  1  of  [24:pp  275-276].  The 
proof  is  patterned  after  their  proofs. 

The  problem  of  interest  is 

+  A(t,u(t))u  =  0;  u(0)  =  uo.  (100) 

The  theorem  identifies  sufficient  conditions  to  guarantee  the  existence  of  a  unique  solution. 
Numerous  hypotheses  are  required  to  describe  the  setting.  They  are  presented  first.  The 
proof  will  be  presented  immediately  afterward. 

Throughout  the  remainder  of  this  chapter,  the  following  hypotheses  will  be  assumed. 

1.  Let  X  be  a  real,  reflexive,  separable,  Banach  space  with  norm  ||  •  ||x.  Let  Y  be  a 
subset  of  X  which,  when  endowed  with  an  appropriate  norm,  ||  •  ||y,  is  itself  a  real, 
reflexive,  separable,  Banach  space.  Let  ||  •  ||y>  ||  •  ||jy. 

2.  Let  W  be  an  open  set  in  Y. 

3.  Let  Q(t,w)  be  a  collection  of  isomorphisms  of  Y  onto  X  for  ( t,w )  £  [0,T]  X  W. 
Assume  T  >  0. 

4.  There  are  real  constants  A q,  Xq,  and  HQ  such  that  ||Q(i, w)||y_*x  <  Aq, 
||<9(^,w)-1|U_fy  <  Xq,  and  ||<?(i,w)  -  Q(i,w)\\Y->x  <  Hq{\  t-i  \  +||w  -  u>||x) 
for  arbitrary  (t,w),  (t,w)  £  [0,T]  X  W. 

5.  Let  N(X)  be  the  collection  of  all  norms  on  X,  equivalent  to  the  given  one,  ||  •  ||x* 
That  is,  N(X)  =  {||  •  ||M  :  ||  •  ||M  is  equivalent  to  ||  •  ||x>  where  n  comes  from  some 
index  set}.  Let  a  metric  (this  will  be  validated  shortly)  for  N(X)  be  given  by 

411-iui  •«,)='<*  js,— {fit*  fife}'  (101) 

6.  Let  N(t,w )  :  [0,T]  x  W  -*  N(X)  be  a  function  satisfying  d(N(t,w),  j|  •  ||x)  < 
A/y/,  d(N(t,w),N(i,w ))  <  ^(|  t  —  t  |  +||u;  -  u)||x)  for  fixed  nonnegative  A jv,  hn 
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G  §£  and  arbitrary  (t,tu),  (t,  w)  G  [0 ,T]  X  W.  Let  denote  the  space  X  with 

the  norm  N(t,w). 

7.  Let  {A.(t,tu) :  (i,tu)  G  f  G  [0,T]  x  W}  be  a  family  of  operators,  A(t,w)  G  B(Y,X). 
There  are  real  positive  constants  Xa  a?  i  M4  such  that  for  each  (t,  w)  G  [0,  T]  x  W, 
the  following  hold: 

A(t,w)  G  G(XN(t,w)>  1»  0) 

M(i»«)  -  <  Va  (|<  -  i\  +  ||t»  -  w\\x) 

8.  Q(t,w)A(t,w)Q(t,w)~1  =  A(t,w)  for  each  (t,w)  G  [0,T]  x  W. 

Theorem  48  When  the  above  eight  hypotheses  are  satisfied,  the  following  conclusion  holds: 
For  each  u(0)  =  «0  £  W,  there  is  some  T  >  0  and  a  unique  solution  u  to  ut+A(t,  u(t))u  =  0 
with  u(0)  =  u0  such  that  u  G  C([0,T];  W)  n  C1([0, 2s);  -ST). 

The  metric  in  Hypothesis  5  appears  in  [24:pg  275]. 

8.5  Lemmas  for  use  in  the  proof  of  Theorem  48 

The  proof  requires  several  technical  results.  It  will  be  useful  to  begin  with  an  overview 
of  the  i  "eas  involved. 

The  solution  which  is  to  be  obtained  will  be  Y- valued.  This  means  that  only  points 
in  Y  will  be  considered  for  values  of  u(t).  A  candidate  solution  can  be  thought  of  as  a 
curve  in  the  space  Y  with  initial  point  uq.  There  is  no  guarantee  that  solutions  can  be 
propagated  for  long  periods  of  time.  Hence,  consideration  is  given  to  candidate  solutions 
which  are  supposed  to  be  valid  until  some  time  denoted  T  which  is  not  yet  specified. 

The  proof  then  is  roughly  outlined  in  the  following  steps.  First,  a  set  E  is  formed 
which  contains  tlm  '••’ndidate  solutions.  Elements  of  E  are  curves  in  Y  with  initial  point 
uq.  It  will  be  established  that  the  chosen  set  is  a  complete  metric  space  under  an  appro¬ 
priate  metric.  Second,  it  is  shown  for  each  fixed  v  G  E,  that  the  family  {A(t,u(t))} 
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has  an  associated  evolution  system  (which  generates  a  solution  to  the  linear  prrblem 
ut  +  A(t,v(t))u  =  0).  Third,  it  is  shown  that  the  solution  to  the  linear  problem  is,  in 
fact,  in  E.  Also,  the  mapping  from  E  to  E  thus  established  is  a  contraction  mapping. 
Then  its  fixed  point  is  the  desired  solution. 

There  are  many  details  between  this  outline  and  the  completion  of  the  proof.  They 
will  be  presented  now. 

3.5.1  Preliminary  lemmas 


Lemma  49  The  proposed  metric  for  N(X)  in  Hypothesis  5  is  valid. 


Proof:  Validation  is  straightforward  from  the  definition,  eg  [54:pg  45],  [66:pg  27],  or  [73:pg 
4].  In  particular,  if  |j  •  ||M  =  ||  •  ^  then  it  is  clear  that  d(\\  •  ||^,  ||  •  ||„)  =  0  since  log  1  =  0. 
Also,  since  jj^  or  its  reciprocal  will  always  be  greater  than  or  equal  to  one,  it  is  clear 
from  the  properties  of  the  log  function  that  d(\\  •  ||M,  ||  •  ||„ )  >  0.  Similarly,  if  ||  •  ||  ^  ||  •  H,,, 
then  d(\\  ■  ||  ,|[  •  !]„)  >  0.  Symmetry  is  obvious.  The  triangle  inequality  is  a  little  tedious 
but  straightforward.  The  details  follow. 


> 


d(\\  •  IU II  •  llv)  +  ^(11  •  ll«/>  li  •  llo)  =  log  sup  max 


Nk  MA 

\\yV\\y\\A 


+log  sup  t rax 
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fjMk  Mil 
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log  sup 

O^yZY 


(mariMuMk  MiiM  11 

V  1  IlirlU  ibllo  ’  IMU  \\v\\v  ’  IMU  IMlo  ’  IlsrlU  M»  j  / 


> 


log  sup 
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{  1  lli/SU ’  llytl,2  ’MMo'MA) 
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For  use  in  the  proof,  a  collection  of  possible  solutions  is  formulated  as  a  metric  space. 
None  of  the  candidates  is  allowed  to  go  too  far  from  the  initial  condition  vector. 

Let  uo  6  W  be  given  with  M  =  ||no||y.  Choose  p  >  0  with  p  <  oo,  such  that 
p  <  dist(uo,Y\W).  (The  symbol  ‘\’  is  used  for  set  subtraction.)  It  is  advantageous  in 
applications  to  choose  p  as  large  as  possible  within  the  bounds  indicated.  The  set  T\W  is 
nonempty  since  W  is  an  open,  proper  subset  of  Y. 

Definition  50  Let  E  be  the  set  of  all  functions  v  :  [0,T]  — ►  Y  such  that  for  each  t,t  £ 
[0,T],  |KO-«o||y  ^  and  ||v(t)  -  ,y(t)||x  <  L\t  -  Z|  where  L  =  M)XqXqc>'n 

and  the  value  of  T  is  not  yet  specified. 

The  purpose  of  these  requirements  is  to  keep  v(t)  in  W.  The  value  of  p  has  already 
been  chosen.  The  choice  of  T  is  constrained  by  (148)  and  the  discussion  in  the  proof  of 
Lemma  60.  These  are  constraints  which  guarantee  that  a  certain  mapping,  to  be  developed 
in  the  proof  of  the  theorem,  is  a  contraction  mapping.  Details  concerning  the  choice  of  T 
will  be  provided  after  the  mapping  has  been  presented.  But  first,  there  are  several  more 
preliminaries. 

Define  a  distance  function  for  pairs  of  elements  in  E  as  follows.  For  v,  w  £  E 

d(v,w)  =  sup  ||v(<)  —  w(t)|U-  (1°3) 

0<(<f 

Lemma  51  The  set  E  with  the  distance  function  d  is  a  metric  space. 

Proof:  The  verification  that  this  is  a  legitimate  metric  is  completely  straightforward.  It 
is  easy  to  see  that  d(v,v)  =  0.  On  the  other  hand,  suppose  d(v,w )  =  0.  Tiien  for  each 
t  £  [0,T],  || v(t)  -  w(<)IU  =  0.  Since  ||  •  ||x  is  a  norm,  it  follows  that  v(t)  =  w(t)  for  each 

i  C.  fO  *P 1  T'hot  ic  -»t  —  on  \A/hon  oi  -£  on  tlioro  ic  cnme  i  curb  that  0)(i\  rk  V)(f  ^  HThon 

v  v_  ^  5  ^  j  •  ~  — . ~  /  —  j - ^  w - - - - ~  v  /  /  "  \  ‘  - - 

d(v,w)  >  sup  ||v(/)  -  w(t)||x 

<€[0,t) 

>  IKO-MOlU 

>  0  (104) 
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since  ||  •  ||x  is  a  norm.  Symmetry  is  obvious.  The  triangle  inequality  is  the  only  remaining 
issue.  Let  u,  v,  w  E  E  be  given.  Then 

d(u,w)  =  sup  \\u(t)  -  w(t)\\x 
<6{o,  t] 

=  sup  ||u(t)-u(t)  +  u(t)-w(t)||x 
<€[0,71 

<  sup  ||«(t)  -  +  sup  \\v(t)  -  w(t)||x 

<€[0,f]  <6  [0,71 

=  d(u,v)  +  d(v,w).  (105) 

Thus,  E  is  a  metric  space.  □ 

It  will  be  important  for  E  to  be  a  complete  metric  space. 

Lemma  52  The  set  E,  with  the  metric  defined  in  (103),  is  a  complete  metric  space. 

Proof:  Let  {«„}  be  a  Cauchy  sequence  in  E.  Then  for  each  t,  {un(t)}  is  a  Cauchy  sequence 
in  X  as  is  easily  seen  by  examination  of  the  metric  on  E.  Since  X  is  a  complete  space, 
there  is  some  v(t )  E  X  such  that  vn(t)  — >  v(t)  in  the  X  norm.  It  is  easy  to  see  that  v  is 
unique.  It  is  not  at  all  clear  whether  v  lies  in  Y  or  is  an  element  of  E.  These  issues  are 
addressed  next. 

According  to  a  standard  theorem,  (eg  Theorem  7.70  of  [20:pg  204]),  {un(0)  treated 
as  a  sequence  in  Y  has  a  weakly  convergent  subsequence,  say  {2n(*)}-  Then,  (see  Theorem 
7.6b  of  [20:pg  202])  there  is  some  v  E  Y  such  that  zn  — >  v.  If  v  =  v  then  no  other 
subsequence  could  converge  to  any  other  point.  But  this  is  easily  shown.  Note  that 
zn  -»  v  in  the  X  norm.  Thus  ]|u(t)  -  u(t)j|A  =  0.  This  mea  is  that  v  =  v  pointwise.  Thus 
v  is  thu  only  weak  limit  of  {vn}. 

If  v  satisfies  the  two  r .  |  ments  on  elements  of  E,  then  the  argument  for  complete¬ 
ness  wiii  be  finished.  But, 


II v(t)  -  =  |K*)  -  vn(t )  +  vn(t )  -  ?i0||y 

<  \\v(t)  -  vn(t)\\Y  +  |M0  “  «o|| Y 

<  e  +  p  (106) 
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for  any  e  >  0  when  n  is  sufficiently  large.  Thus  the  first  requirement  is  satisfied.  Further¬ 
more, 


=  v(t)  -  Vn (i)  +  »„ (t)  -  V„(i)  +  Vn(t)  -  u(t)|^ 

<  IKO  “  Mt)\\x  +  ||»«(0  “  +  ||v»(*)  “  v(*)\\x 

<  e/2  +  L\t  -  t|  +  e/2 


(107) 


for  any  e  >  0.  Thus  the  second  requirement  is  satisfied. 

This  completes  the  argument  for  the  completeness  of  the  metric  space  E.  □ 

3.5.2  Preparing  for  an  evolution  system 

The  proof  of  the  theorem  will  require  the  solution  of  a  sequence  of  linear  equations. 
For  each  linear  equation  a  solution  is  obtained  from  Theorem  46.  In  this  subsection  the 
linear  equations  are  introduced  and  it  is  verified  that  they  satisfy  the  hypothesis  of  Theorem 
46. 

Some  convenient  notation  is  now  introduced.  For  each  v  G  E,  let  Nv(t)  =  N(t,  v(t))  = 
II '  ll<)  Qv(] 0  =  Q(i*  t>(t)),  and  Au(t)  =  A{t,v(t)).  Throughout  the  remainder  of  Section  3.5, 
v  will  always  represent  an  element  of  E. 

Consider  now,  for  fixed  v,  the  (still)  nonautonomous  (but  now)  linear  problem 


(I'll  A 

—  +  Av(t)u  =  0;  0  <t  <T;  u(0)  =  u0.  (108) 

( XL 

The  next  major  step  is  to  establish  the  existence  of  an  evolution  system  which  solves 
this  linear  problem.  The  lemmas  which  follow  are  necessary  to  establish  the  hypotheses  of 
a  Theorem  46. 


Lemma  53  The  family  (Au(t)}  C  G(X 1»  0)  is  stable  with  stability  constants  Mux  = 
eXN+2ixN(i+L)f  f  and  uj  =  o. 


I1  for  the  given  v  and  arbitrary  t,t  £  [0,  T].  First, 


Proof:  It  is  important  to  bound 
consider  that  Hypothesis  6  and  the  second  property  of  E  reveal  a  bound  for  the  distance 
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between  the  two  norms. 


d(IMUHIt)  ^  VN{\t-i\+\\v(t)-v(t)\\x) 

<  fiN(\  t-i\+L\t-t\) 

=  pN(l  +  L)\t-t\ 


(109) 


From  the  definition  of  the  distance  function  on  the  collection  of  norms  (see  Hypothesis 
5)  it  follows  that 


log  sup  maiB,  jnr\  -  Wv(l  +  L)  \  t  -  i 

o#yex  IIMIt  lll/llt  J 


so  that 


Mi  <  „„J IMk  Mi' 

Mi  -  \M(’  Mi. 


(110) 


and  the  desired  bound  is  c/‘n(i+£)I<-<I. 

A  method  from  the  proof  of  Proposition  3.4  of  [31:pg  245]  is  used  in  the  next  step. 
From  the  definition  of  a  stable  family  (Definition  35)  and  the  equivalence  of  norms,  write 


IlTO+^r1* 

i=i 


<  e^{nL)(f-tk)  JJ  (Av(tj)  +  A)-1  y 


lli=i 

Wk-1 


Ik 


A 


<  l.e^N{^+L)(f-tk)enN(l+L){tk-tk^i) 

A 


n  (A«fo)  +  A)-1  y 
f=i 


tk 


fl(A'’(tj)  +  \)-'y 

3= 1 


/ 1  \  k 


<  j  eVN(l+L)(T-ik)eHH(l+L)(tk-tk- 1)  . .  .  e/itf(l+L)tj  ||y||0 


<  (jJ  ew(1+L)r||i/||o 

<  (i)fce^(1+L)(t+3‘,)||2/||3-, 
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(Ill) 


<  (I)‘<*»<>+«*|»||1, 

Clearly  this  bound  holds  for  any  of  the  norms.  Note  also  that  A  =  1  is  a  valid  choice  in 
the  current  setting. 

This  yields  stability  of  the  family  {^"(t)}  with  stability  constant  M  =  e2ttN^1+^  for 
any  Xt,  meaning  the  set  X  with  any  one  of  the  norms  N(t,v(t)).  It  remains  to  translate 
this  to  the  original  norm  on  X.  But,  from  Hypotheses  5  and  6,  it  is  easy  to  establish 
that  ||  •  <  eAjv ||  •  ||t.  Namely,  from  the  definition  in  Hypothesis  5  and  the  bound  in 

Hypothesis  6  it  follows  that 

wli’ll}  =  ^■|li'IMU) 

<  A  N.  (112) 


Then 


sup  max 

0?y(-X 


(Ml  M*\ 
l MU’  Mti 


<  eXN. 


(113) 


Now  it  is  clear  that  <  eAw  and  hence  ||?/||^  <  eAN||y||t  for  any  nonzero  y  E  X.  Then 
(Av(t)}  is  a  stable  family  in  X  with  stability  constant  Mux  =  eXNM  =  e*w+2Mw(i+£)i’  as 
desired.  □ 


Lemma  54  For  each  v  E  E,  the  mapping  t  ->  Av(t)  E  B(Y,X)  is  Lipschitz  continuous. 


Proof:  The  proof,  which  is  straightforward,  is  outlined.  Let  i,t  E  [0,T],u  E  E  be  given. 
Then 


\\A'’(t)-A'’(i)\\Y-.x 


= 

II  A(t,v(t))~  A(t,  u(f))||y_x 

< 

II  » 

( \i  —  f  1  -L  1 1  ■/Ml 

— 

/~/i 

\r  vi  1  ii  ~  w  ~  wii^  j 

< 

PA 

(\t-i\  +  L\t-i\) 

< 

PA 

(l  +  L)\t-i\.  □ 

Lemma  55  For  each  t  €  [0,T],  v  E  E,  the  bound  ||<2vC0||y-*;r  <  A q  holds. 


(114) 
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3.5.3  The  evolution  system 

The  hypotheses  of  Theorem  46  are  now  satisfied  and  the  theorem  gives  an  evolution 
operator  Uv(t,  s).  This  is  a  good  time  to  review  how  the  proof  all  fits  together.  Based 
on  the  previous  discussion,  an  evolution  operator  can  be  obtained  corresponding  to  If 
the  image  of  elements  of  W ,  under  the  action  of  the  operator  Uu°,  is  again  in  E,  then 
an  iterative  procedure  is  justified.  This  subsection  consists  of  preliminaries  necessary  to 
establish  that  such  images  are,  in  fact,  again  in  E. 

Now  that  the  evolution  operator  has  been  identified,  it  will  be  useful  to  identify 
bounds  for  it.  This  will  require  attention  to  a  portion  of  the  proof  of  Theorem  44.  Bounds 
are  desired  for  ||17v||^  and  ||17t'||y. 

Lemma  58  For  each  v  G  E,  \\UV\\X  <  ex^^1+L)f  . 

Proof:  The  bound  for  ||f/"||x  is  immediate  from  Lemma  53  and  property  1  of  Theorem 
44.  In  particular,  property  1  of  Theorem  44  says  that  ||17i;(i,s)||x  <  Muxew^~sh  Lemma 
53  says  w  =  0  and  Mux  has  the  value  indicated.  □ 

A  bound  for  ||17v||y  is  not  as  simple  to  obtain.  The  bound,  which  is  established  in 
the  next  lemma,  will  be  denoted  by  Muy. 

Lemma  59  For  each  v  G  E, 

m\y  <  AgAgexp  [(A*  +  2/^(1  +  L)T)  +  /ig(l  +  L)\qT e(^+Wi+i)r)] .  (118) 

Proof:  The  strategy  is  to  define  an  intermediate  operator  V"  such  that  Uv  =  (Qv)~lVvQv. 
Notice  that  while  Uv  :  Y  ->  Y  in  this  equation,  the  isomorphism  Qv  allows  Vv  to  be  a 
mapping  from  X  to  X.  A  bound  on  ||V'v||^-  can  then  be  multiplied  by  AqAq  to  give  a 
bound  for  1117%. 

Let  Cv(t )  =  Qv(t)(Qv(t))~1  where  Qv(t )  is  the  same  as  discussed  in  conjunction  with 
Theorem  44  with  the  choice  Q  =  I  +  A.  The  intermediate  operator  Vv  will  be  defined 
in  terms  of  Uv  :  X  -*  X.  The  operator  Vv,  here,  is  given  by  the  V  of  Lemma  43.  From 
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Proof:  This  simply  requires  a  careful  look  at  the  hypotheses.  Let  t  €  [0, 2s],  v  €  E  be 
given.  Then 


||Qu(0l|y-*x  =  ||Q(i>v(0)lli'-or 

<  Aq  (115) 

where  the  inequality  holds,  for  each  t ,  by  Hypothesis  4.  Hence,  since  the  bound  is  uniform, 
it  also  holds  for  the  supremum  over  all  t.  □ 

The  case  for  ||<2w(0‘"1IU-»y  <  A q  is  similar. 

Lemma  56  The  bound  ||Qv(t)  —  Ql'(t)||y_^  <  hq(  1  +  L)  \  t  —  i  |,  holds  for  i,t  €  [0,T] 
and  each  v  £  E. 

Proof:  The  proof  is  a  straightforward  application  of  Hypothesis  4  and  the  definition  of  E. 

mo  -  =  II Q{t,  HO)  -  Q(.i  v(i))\W~x 

<  Mq( I  t~i\  +IM0  -  v(i)\\x) 

<  Mli-il+Llt-il) 

=  »Q(l  +  L)\t-i\  .  O  (116) 

Lemma  57  For  each  t  £  [0,T]  and  v  <E  E,  Qv(t)Av(t)Qv(t)~1  =  Av(t). 

Proof:  This  is  almost  immediate  from  Hypothesis  8,  which  holds  for  each  t. 

Qv(t)Av(t)Qv(t)~1  =  QitAmitxwMt))-1 
=  A(t,v(t)) 

=  Av(t).  □  (117) 

This  completes  the  preliminaries  necessary  to  obtain  an  evolution  system  corresponding 
to  the  linear  equation  108. 
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3.5.3  The  evolution  system 

The  hypotheses  of  Theorem  46  are  now  satisfied  and  the  theorem  gives  an  evolution 
operator  Uv(t,s).  This  is  a  good  time  to  review  how  the  proof  all  fits  together.  Based 
on  the  previous  discussion,  an  evolution  operator  can  be  obtained  corresponding  to  u o.  If 
the  image  of  elements  of  W,  under  the  action  of  the  operator  Uu°,  is  again  in  E,  then 
an  iterative  procedure  is  justified.  This  subsection  consists  of  preliminaries  necessary  to 
establish  that  such  images  are,  in  fact,  again  in  E. 

Now  that  the  evolution  operator  has  been  identified,  it  will  be  useful  to  identify 
bounds  for  it.  This  will  require  attention  to  a  portion  of  the  proof  of  Theorem  44.  Bounds 
are  desired  for  ||I/v||jr  and  ||f/u||y. 

Lemma  58  For  each  v  6  E,  Hit'll*  <  eAw+2*‘N(1+i')?’  . 

Proof:  The  bound  for  ||17,'|]x  is  immediate  from  Lemma  53  and  property  1  of  Theorem 
44.  In  particular,  property  1  of  Theorem  44  says  that  ||Z7v(i,  ^)||x  <  Muxew^~s\  Lemma 
53  says  w  =  0  and  Mux  has  the  value  indicated.  □ 

A  bound  for  ||f/”||y  is  not  as  simple  to  obtain.  The  bound,  which  is  established  in 
the  next  lemma,  will  be  denoted  by  Muy- 

Lemma  59  For  each  v  €  E, 

i|tfv||y  <  AgAgexp  [(A*  +  2W(1  +  L)t)  +  hq(1  +  L)XQTe^+2^+L^] .  (118) 

Proof:  The  strategy  is  to  define  an  intermediate  operator  Vv  such  that  Uv  =  (QV)~1VVQV. 
Notice  that  while  Uv  :  Y  — »■  Y  in  this  equation,  the  isomorphism  Qv  allows  Vv  to  be  a 
mapping  from  X  to  X.  A  bound  on  j|Vv||;f  can  then  be  multiplied  by  XqXq  to  give  a 
bound  for  ||t/vj|y. 

Let  Cv(t)  =  Qv(t)(Qv(t))~1  where  Qv(t )  is  the  same  as  discussed  in  conjunction  with 
Theorem  44  with  the  choice  Q  =  I  +  A.  The  intermediate  operator  Vv  will  be  defined 
in  terms  of  Uv  :  X  -*  X.  The  operator  Vv,  here,  is  given  by  the  V  of  Lemma  43.  From 
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Lemma  43,  if  H(t)  =  C(t )  =  <5(t)(<3(t))-1,  then  for  any  x  e  X,VV  is  the  unique  solution 
of  the  integral  equation 


Vu(f,  s)x  =  Uv(t,  s)x  +  J  Vv(t,r)C(r)Uv(r,s)x  dr. 


(119) 


It  is  shown  in  the  proof  of  Theorem  44  that  Vu,  defined  in  this  way,  satisfies  Uv  = 
(QV)~1VVQV .  For  completeness,  the  relevant  details  are  presented.  (For  convenience  in 
writing,  the  superscripts  V  will  be  omitted.) 

From  the  boundedness  of  Q-1  and  the  differentiability  of  Q  it  follows  that  Q~l  is 
differentiable  a.e..  Note  that 

jt  (< Q(tr **)  =  -Q(trlQ(twrlx-  (120) 

As  a  temporary  notational  convenience,  let  F(t,r)  =  U(t,r)Q(r)~1.  Then 

~F{t,r)x  =  7^U(t,r)Q(r)~1x  +  U(t,r)~Q{r)~lx 

=  -U(t,  r)A(r)Q(r)~lx  -  U ( t ,  r)Q(r)-1Q(r)<2(r)-1a: 

=  -  (F(t,r)A{r)+ F(t,r)Q{r)Q{r)-lS)x  (121) 


where  the  last  step  uses  the  commutativity  of  an  operator  and  its  resolvent.  Let  Un(t,  r )  — + 
U ( t ,  r )  so  that  for  any  y  €  Y 


d_ 

dr 


Un(r,s)y=  An(r)Un(r,s)y. 


(122) 


Now, 


JzU{t,r)Q{r)  1Un(r,s)y 


v/  i 


0 

—F{t,r)Un(r,s)y 

i 

- F(t,r )  (A(r)  +  Q(r)(5(r)-1)  Un(r,s)y 
+F(t,  r)An(r)Un(r,  s)y 

-F(t,r)  (A(r)  +  Q(r)g(r)-J  -  An(r))  Un(r,s)y(m) 
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The  next  step  is  to  integrate  the  last  expression  from  s  to  t  with  respect  to  r.  This  will  be 
displayed. 


U(t,r)Q(r)  1Un(r,s)y\ta  =  -  I  F(t,r)  (A{r)  +  Q{r)Q(r)  1  -  An(r))  Un(r,s)ydr 

J  S 

Q(t)~1Un(t,s)y  -  U(tis)Q(s)~1y  =  -J  U(t,r)Q(r)~1Q(r)Q(r)-'[Un(r,s)ydr 

+  ^  U(r,t)Q(r)~1(An(r)-  A(r))Un(r,s)ydr  (124) 

Note  that  the  second  integral  vanishes  in  the  limit.  Then,  from  the  denseness  of  Y  in  X, 
it  follows  that 

-Q(t)-1U(t,s)x  +  U{t,s)Q(s)~1x  =  Jl  U(t,r)Q(r)-1Q(r)Q(r)~1U(r,s)x  dr  (125) 
or 

U(t,s)Q(s)~1x  =  Q(t)~1U(t,s)x  +  jf  U(t,r)Q(r)-1Q(r)Q{r)-'U(r,s)xdr.  (126) 

It  is  useful  to  compare  this  to  the  form  of  (119).  Specifically,  apply  <5_1(t)  to  the  left  side 
of  each  operator  in  (119).  Then  it  becomes 

Q~l(t)V(t,  s)x  =  s) x  +  J*  Q-  K>)V(t,  r)C{r)U{r , $)x  dr.  (127) 

It  is  now  clear  that  Q~1(t)V(t,s)x,  and  U(t,s)Q(s)~1x  are  each  solutions  to  the  same 
integral  equation,  which  is  known  to  have  a  unique  solution.  Hence 


V  =Q~1VQ 


(128) 


as  desired. 

A  construction  of  V  is  now  gi ,  .  The  purpose  of  this  construction  is  to  obtain  a 

bound  on  || V'Hjf .  Let 

V(o  )(t,r)x  =  U(t,r)x  (129) 
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It  is  easy  to  see  that. 

I|V(0)IU  <  Mux.  (130) 

Next, 

V(i)(i,r)a:  =  V(0)(t,s)C(s)U(s,r)x  ds 

=  J  V(0){t,s)C(s)V(0)(s,r)xds 

llv(i)(*»r)*l|jc  =  v(o){i>s)c(s)v(o)(s>r)xds 

<  j'\\V{o)(t,s)C(s)V{0)(s,r)x\\xds 

<  Jr  l!V(o)(^^)!UI|C(s)V-(0)(s,r)x|UdS 

<  MUxJr  ||C(s)||A’||V(0)(s,r)a:||xds 

<  MUx  IIC'Hoo  J*  ||V(0)(s,  r)|U||x||x  ds 

<  Mu^WCWooWxWx  J*  ds 

<  Mu^WCWMxit-r) 

<  MUx*\\C\Ux\\xf  (131) 


(132) 


(133) 
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The  next  step  is  to  show  that  (133)  holds  for  V(*+1).  Then, 

ll^(fc+i)(*ir)®ll*  =  | Jf  V(k){Us)C(s)V(0)(s,r)xds 

<  Jr  *)C(*)' V(0)(5,r)a:||A-ds 

<  Mc/xfc+1||C||^^  \\C(s)V{0)(s,r)x\\x~^-ds 

rf,k+ 1 

<  Afir**+3||Cr||*+1  H.IIjp (134) 

The  operator  V  obtained  by  the  construction  is 

V(t,s)  =  '£vm(t,s).  (135) 

fc= 0 

It  is  appropriate  to  verify  that  this  expression  for  V  agrees  with  (119)  of  Lemma  43.  The 
issue  is  whether  V  satisfies  (119)  since  (119)  is  known  to  have  a  unique  solution.  So, 
consider 

00  -t  oo 

Ev(*)(*,s)®=  /  J2v{k)(^r)c(r)v{ o)(r,s)xdr  (136) 

k=i  Js  k=0 

where  the  k  =  0  term  from  the  left  hand  side  has  cancelled  with  the  first  term  on  the  right 
hand  side  of  (119).  Consider  the  first  remaining  summand  on  each  side.  By  definition, 

V(i)(t,  s)x  =  jf  V(o )(t,  r)C,(r)V(0)(r,  s)x  dr.  (137) 

Likewise,  for  each  succeeding  pair  of  terms,  the  summand  on  the  left  equals  the  summand 
on  the  right.  Hence,  the  sums  are  the  same. 

A  bound  for  V  is  available  in  terms  of  the  bounds  on  the  Vy.y  Namely, 

nu  <  £iivwiu 

k= 0 

OO  rpk 

<  E^H1M~Tr 

k= 0  K' 
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(138) 


<  MUxeMVx\\oUt 

where  the  standard  series  representation  for  e  is  easily  recognized. 

Everything  is  in  place  to  bound  ||I/]|k.  The  superscripts  will  now  be  used  again. 

l|H"!|y  =  IWT’V’Q’lli' 

<  IIMT'llx-.HIHIxlieli'-x 

<  M«IIHIx 

<  \q\qMUx  llc”ll»:f  (138) 


Recall  that  Cv  •- 

The  final  expression  can  be  made  more  explicit  in  terms  of  fundamental  quantities. 


\\Cv(t)\\oo  =  ess  sup  II Qv{t)Qv(t)  *|U 
<6(0, f) 

<  ess  sup  ||<5v(t)||y_>r 

«e[o,f) 

<  ess  sup  hq{  1  +  L)Xq 

t€[0,T] 

=  1  +  L)Xq  (140) 


The  combination  of  estimates  (139)  and  (140)  gives 


ll^lly  <  AgAgexp  [(A*  +  2^(1  +  L)T )  +  /zg(l  +  L)XqT e(*»+Wi+i,)f) 


(141) 


The  bound  given  by  the  right  hand  side  of  this  equation  will  be  denoted  by  M\jy.  This 
concludes  the  proof  of  the  lemma.  □ 

Tlm<;  hounds  on  the  evolution  onprators  are  obtained  and  t.lip  preliminaries  for  es- 
tablishing  a  contraction  mapping  are  complete.  Note  that  the  linear  problem  (108)  is  all 
that  has  been  solved  at  this  point.  But,  as  it  turns  out,  the  solution  is  in  E  and  a  mapping 
from  E  —>  E  is  established.  Thenit  will  be  established  that  the  mapping  is  a  contraction. 
The  fixed  point  of  the  contraction  map  is  the  desired  solution. 
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3.5.4  A  contraction  mapping 

Attention  now  returns  to  the  solution  of  (108)  in  terms  of  Uv(t,s).  In  particular, 

uv(t)  =  Uv(t,  0)tt0.  (142) 

So,  uv(t)  is  a  continuous  mapping  of  [0,  T]  into  Y  by  properties  4  and  5  of  Theorem  44.  By 

property  2  it  is  differentiable  for  small  t.  Also,  from  the  definition  of  a  Y- valued  solution, 

it  is  continuously  differentiable  into  X. 

It  is  desirable  to  establish  that  uv  £  E.  This  is  the  point  in  the  argument  where  the 
choices  of  L  and  T  occur.  The  conditions  that  uv  must  satisfy  are: 

IK(<)  -  uo||k  <  P  (143) 

\\uv(t)  -  tt*(<)Ux  <  L\t  -  t\.  (144) 


It  will  now  be  shown  that  suitable  choices  are  available  to  make  this  so. 
Lemma  60  For  each  v  £  E  and  t  £  [0,T],  uv(t)  6  E. 


Proof:  Consider  the  inequality  (144)  first.  Note  that  ^uu(t)  =  —Av(t)uv(t).  Recall  that 
M  =  HuoIIk  and  R  is  the  radius  of  a  ball  centered  at  uq.  Then, 


uv(t)  -  uv(£)||  <  sup 

X  tg[o,il 


< 


d  .  . 


I  <-<l 


=  \\-A'>(t)u'’(t)\\x\t-i\ 

<  il«v(t)||y|t-t] 

<  A^||17v(t,0)«o||x|t  - 1| 

<  A^||^||y  ||«o||y|«-i| 

<  xA(R  +  M)\\uv\\Y\t-i\ 

Xa(R  +  M)\qXqgxp  [(Aat  +  2pN(l  +  L)T) 
+pQ{  1  +  L)XQfe^+2fi»(1+tW}  1 1  -  t\. 


(145) 
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The  goal  is  to  have  this  quantity  less  than  or  equal  to  Iu  —  <|.  But, 

L  =  2Xa  (R  +  M) %Q\QeXN  .  (146) 

Clearly  T  can  be  chosen  small  enough  for  the  condition  to  be  satisfied. 

The  inequality  (143)  is  satisfied  for  some  f  >  0  since  uv(0)  =  «o  and  uv  is  continuous, 
hence  uv  £  E.  □ 

Thus,  a  mapping  §  :  E  E  is  identified.  In  particular,  for  any  v  €  22,  §(u)  is 
obtained  by  the  following  steps: 

1.  Obtain  an.  evolution  operator  Uv(t,s)  for  the  problem  $-tu  +  Av(t)u  =  0. 

2.  Let  uv(t)  =  Uv{t,0)u o. 

3.  5>(u)  =  uv. 

The  next  item  to  establish  is  that  $  is  in  fact  a  contraction  map.  This  will  require 
a  prelim5  nary  lemma. 


Lemma  61  Jet  V  and  V  be  evolution  operators  corresponding  to  Au  and  Av  respectively. 
Then 

V{itr)y- J  U(t,s)(Au(s)  -  Av($))V($,r)yds .  (147) 

Proof:  See  Appendix  D.  □ 

The  lemma  will  now  be  used  to  establish  that  $  is  a  contraction  mapping. 

Lemma  62  The  map  $  :  E  ->  E  is  a  contraction. 

Proofs  The  argument  is  straightforward.  Let  u,  v  G  E  '>e  given.  Let  u  =  v  =  $v. 


d(u,v)  =  sup  jju(t)  -  6(i)|U 
o  <t<f 


=  sup  ||f7u(t.,0)uo  -  ^v(t,  0)«olU 
o  <t<T 
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< 

< 

< 


sup 

0  <t<T 


Uv(t,s)(Au(s)  -  Av(s))Uu(s,[})u0ds 


x 


sup  [  \\Uv(t,s)(Au(s)~  Av(s))Uu(s,0)uo\\xds 
o  <t<f  Jo 

sup  MUx2 /iAd(u,v)\\u0\\xt 
o  <t<T 

Mux2HAd(u,v)\\uQ\\xt 


(148) 


where  Myx  is  taken  from  Lemma  53. 

It  is  now  clear  that  T  can  he  reduced  if  necessary  to  ensure  that  $  is  a  contraction 
mapping.  □ 


3.6  The  proof  of  the  theorem 

The  Contraction  Mapping  Theorem,  eg  [54:pg  126],  [13:pg  181],  or  [46:pp  40-42], 
applies  to  give  the  desired  soluti.  and  the  theorem  follows  as  outlined  at  the  beginning 
of  Section  3.5.  □ 


5.7  Applicability  of  Theorem  f  8  to  an  a (t)  case 


Consider  utt  +  a(t)uxxxx  =  0  with  the  abstract  formulation  ut  +  Aju  =  0;  u(0)  =  uo 


where 


An  = 


0 

a(t)D4 


(149) 


As  before,  a(t)  >  am;n  >  0  for  all  t  G  [0,T].  Furthermore,  require  a(t)  to  be  continuously 
differentiable.  Thus  there  is  some  amax  such  that  a(t)  <  amax  for  all  t  £  [0,T].  Similarly 
there  is  some  a'max  such  that  a'{t)  <  o!max  for  all  t  £  [0,T].  For  this  equation,  it  has  simply 
been  assumed  that  a  could  vary  with  t  in  a  model  originally  derived  for  constant  a. 


Notice  that,  in  this  particular  application,  the  operator  does  not  have  any  dependence 
on  the  solution.  Thus  the  full  power  of  the  theorem  is  not  exercised.  The  formal  presence 
of  w  in  A.r(t,  w)  will  nevertheless  be  retained  in  the  following  exposition. 

The  space  X  is  almost  the  same  as  before.  The  point  set  is  unchanged  from  (30) 
and  (3l).  However,  since  the  inner  product  on  X  uses  the  value  of  a  in  its  definition,  it 


73 


would  appear  to  have  a  variable  norm.  Though  it  will  be  useful  to  treat  a  family  of  1101ms 
corresponding  to  the  various  values  of  a  this  is  not  acceptable  for  the  definition  of  the 
Banach  space  X.  Since  a  is  a  continuous  function  of  t  on  a  closed  and  bounded  set,  there 
is  some  tmm  €  [0,T]  such  that  a  attains  its  absolute  minimum,  say  amin.  The  value  amm 
replaces  a  in  the  definition  of  the  inner  product  on  X  as  given  in  Theorem  10.  The  space 
Y  is  D(A$')  with  the  norm  defined  as  follows.  The  norm  on  Y  is  a  graph  norm  which 
depends  on  the  operator  A.  Since  the  operator  is  now  allowed  to  vary  with  >  and  xv  it  is 
necessary  to  be  careful  so  that  the  norm  is  well  defined.  To  this  end,  choose  A-(lm ii0) 
as  the  operator  in  the  norm  for  Y, 

\\v\W  =  Ill/ll  A'  +  p7(fmm>«o)2/||A-  (150) 

The  set  W  is  a  ball  in  Y  with  radius  R  and  center  at  the  initial  condition  «o-  Other 
symbols  with  the  same  meaning  as  before  are  M  =  ||«o||y .  and  Q  =  I  +  A. 

It  is  easily  seen  that  Hypotheses  1-3,  5,  and  S  of  Theorem  48  are  satisfied  from  the 
same  arguments  that  were  used  in  the  proof  of  the  theorem.  For  example.  Hypothesis  8 
is  satisfied  since  Q~l  is  the  resolvent  of  A  and  resolvents  commute  with  their  genciatois. 
The  other  hypotheses  must  be  addressed  individually. 

Consider  Hypothesis  4.  An  upper  bound  on  ||Q(t,  tu)||y— a  is  desired.  Let  values  of 
i  and  xv  be  given.  Note  that  w  does  not  have  any  role  in  this  particular  application.  Then. 

||Q(M«)||v-A  =  ||(/  +  .47(/,w))||y_.Y 

||(/  + A7(f,w))y|lA 

lllfllr 

_  [j(/  +  A~(t,w))y\\x 

MIX  "b  II  A~(tm{n,  Wo).7/||a 

<  Mi +  imtu.wMIa 

-  O&r  Wx  +  Ur (Ixnin  >«o)j/||a 
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f  -y2  \ 

IUIU  +  * 

V  Q!(i)-D  2/1  /  y 

—  SUD  x _ _ _ ±L 

.  OjSygy  llvlU  +  ||A7(imtn,«o)2/|  X 

llsllx  +  (ctm;„  fo  lx")2  dx  +  a(tfti(D4y{f  <fe)'/2 
■  0*61-  llvk  +  M7(<mi».«o)l/|U  (  ' 


The  remainder  of  the  argument  is  carried  out  in  two  separate  cases.  First,  suppose  that 
<*mm  <  1.  Then 


||Q(t,u;)||y_X 


Mix  +  a!?.  fo'  Wf  dx  +  «(!)''  /„'  (!x)'r‘ 

-  o SSr  llvllx  +  l|AT(imi„,«0)j,||A- 

IlylU  +  IS  (iftjs  /o  (ya")2  ^  +  a2min  jj  (D4yi)2  dxj 
jgr  l|y|U  +  IIMW«o)y|U 

\\y\\x  +  (flmin  fo  (y  2"?  dx  +  a2min  /o1  {D*yxf  dx)1'2 
~  oJSey  !MU  +  ||A7(tm,-„,«o)y|U 

IlylU  +  («m,n  fo  (V2")2  dx  +  a2min Si  {D4yx)2  dx)l/2 

-  oSfy  IlylU  +  IIM^fn.ttoJylU 

IlylU  +  tellA7(tm«n>Ho)y|U 

-  SUD  - mm  — 

OfyeY  l|y|U  +  p7^mtn,«o)y|U 

<  222*.  (152) 

"min 


For  the  second  case,  suppose  am,-n  >  1.  Then 

IIQ(<»w)lli'-Jr  £  SUP 

ofyev 

jjyjU  +  "(*)  (gffi/p  fo")2  dx  +  fo  ( DAy\  f  dx )1/2 
O&ZY  l|y|U  +  M7(<mi»,tto)y|U 

IlylU  +  "(0  ("min  /o  fo")2  dx  +  "min  fo  (^Vlf  dx) 

-  oSfy  IlylU  +  P7(imin,«o)y|U 

IlylU  +  "mas  ("min  fo  fo")2  dX  +  a2min  Jp  (£>fo)2  dxf^ 

~  oSey  ||y|U  +  ll^min,«o)y|U 


IlylU  +  ( j§)  (a^  fo  W?  dx  +  Q^)2  fo  W2  dxf/2 
||y|U  +  ll^min,«o)y|U 
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toll*  +  (ami,  fj  (i Iff  dx + <&.  /o1  fe)1/2 

_  amai0|“Py  toll*  +  M7((mta,«o)!/|U 

=  ®mtt  r*  (153) 


Therefore,  choose  Aq  =  max{^f^,amax}. 

A  lower  bound  on  ||<5(t,  w)j|y_^  is  useful  to  help  identify  an  upper  bound  for 
||Q(t,u;)~1||^_+y.  Since  —A7(t,w)  is  dissipative, 


||9(<»w)tflU  =  ([/+  A7(t,vi)]y,[I  +  A7{t,w)]y)1/2 


=  ((y,  y )  +  2  (y,  A7(t,  w)y)  +  ( A7(t,  w),  A7(ty  w))f/2 

>  (lMlx  +  M7(i,w)j/||^)1/2 

>  2-1/2(\\y\\x  +  \\A7(tyw)y\\x) 

=  2~1/2(^\\y\\x+ (aminJ^  (y2")2  dx  +  aitf  (D4yi)2  dx^j  j 

>  2_1/2  (^\y\\x  +  (*rninj  (V2 "f  dx  +  (X2min  jf  (T»41/l)2  j 

=  2-1/2m\x  +  \\A7(tmin,w)\\x) 

=  2-l'2(\\y\\x  +  \\A7(tm{n,u0)\\x) 


=  2-1/2||*,||y. 


(154) 


It  follows  from  Theorem  5.7.1  [54:pg  244]  that  an  acceptable  choice  is  A q  =  21/2. 
The  next  task  is  to  calculate  hq. 


\\Q(t,w)-Q(i,w)\\Y^x 


II  -  Q(t,w))  y\\x 

sup  - ii — fi - - - 

o &&  WvWy 


IIQ(<iw)y-Q(f,w)y|U 

!I?/||v 

II  (I  +  A7(t,  w))  y  -  (/  +  A7(t,  w))  y|U 

sup  - - ' - 

OtvGY  My 
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Hi47(*,w)y-  A7(t,w)y\\x 
o?£y  II* 


.  |j  (A7(t,w)-A7(t,w))y\\x 

sup  — - j t-j. - - - 

OjiyGY  Pl|y 


sup 
0  *y& 


'  0  0 

>  (a(t)  -  a(i))D4  0 


IMIy 


sup 

o^j/ey 


\  (a(t)  “  <*(t))D4p i 


X 


|a(i)  -  a(f)|  sup 

OySj/gy 


0  \ 
D4yi  ) 


IMIy 

/ 


|a(i)  —  a({)j  sup 


DAy i  j 


0#{ ,ZY  ||y|U  +  p7(^m»n,w)T/||x 


|a(i)  -  a(£)|  sup 


o#yey 


V 

U/ 

/ 


/ 


0 

D4yi 


\ 


/ 


\ 


-y  2 

^  &minD  t/i  J 


X 


ja(i)  -  a(£)j  sup 
0  #ygy 


0 

D4y  i 


/ 


0 


\ 


&mtnD  I/I 


* 
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1 


=  |a(t)  -  a(i)  | -  sup 

ttmin  O^y&Y 


l  o  ' 

D4y, 


X 


X 


1 

amin 


=  (155) 

CXmin 

So,  choose  hq  =  This  completes  the  arguments  to  satisfy  Hypothesis  5.  Notice  that 
in  this  application  the  isomorphism  does  not  depend  on  w. 

Consider  Hypothesis  6.  Let  N(t,  w)y  =  a(t)  (yi")2  dx  +  y 22  dx.  Then 


loq  sup  max  i  dx  +  &  ^  dx  Qmi'n &  (yi"?  dx  +  £- y 2 dx  } 

o^yex  \  amin  Jj  (y\")2  dx  +  jg  y22  dx  ’  a(t )  jg  (yi")2  dx  +  jg  y22  dz  J 

log  sup  i— 

1.  a min  J 


<  /oy^£.  (156) 

«min 

So,  choose  A/y  =  log |j^-.  It  is  easy  to  see  that  /z/y  —  Aw  is  suitable  for  the  current 
problem. 

The  next  argument  requires  a  preliminary  lemma.  The  lemma  is  elementary  but  is 
included  to  help  clarify  the  argument. 


Lemma  63  For  positive  real  numbers  a, b, and  c  with  b>  c,  it  follows  that 


Proof:  See  Appendix  D. 


□ 
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Consider  Hypothesis  7.  Claim:  A^t,w)  G  G(X^tiW^,  1,0).  Recall  that  XN(t,w)  is 
simply  the  set  X  with  the  norm  that  uses  a(t).  For  each  (t,w)  this  reduces  to  the  constant 
coefficient  case.  In  particular,  for  each  fixed  value  of  t,  Lemmas  21  and  26  and  Theorem 
27  apply  to  A^t,w)  just  as  they  did  to  A\. 

Claim:  ||A7(f,  <  A^.  This  is  straightforward,  as  follows. 


\\A7(t,w)\\Y^X  = 


p7(t,w)yl[x 

<¥yey  My 


=  sup 


\\M^  W)2/IU 


\\v\\x  +  m7(tmi»»i»)»IU 
\\M^w)y\\x 


<  sup 


0 SZr  WMt  min  5 


I  -»  ' 


=  sup 

QtydY 


a(t)D4yi 


J 


( 


-2/2 

<XminD4yi 


< 


sup  ■ 
Ojij i<-Y 

(&min  Jo 

(y2")2 

dx  + 

o(0 /o1  (n'vtf 

dx'j 

1/2 

(amin  Jo 

(V2")2 

dx  -f 

Ofmin  Jo  (D*yi)2 

dx^ 

,1/2 

sup  ■ 
o;£y€y 

((Xmin  Jo 

(y2"f 

dx  + 

«««./<>  (D*yi)2 

dx 

j  1/2 

(&min  Jo 

(»")* 

dx  + 

(XminjJ  (D4Vl)2 

\ 

dx 

* 

j  1/2 

&max 

a 


(157) 


mtn 


Choose  Xa  =  75m?t. 

wmm 

It  is  suitable  to  choose  ha  =  Hq-  Continuity  holds  as  argued  in  the  /3(t )  case. 

This  completes  the  verification  of  hypotheses  lor  this  case.  The  theorem  applies  and 
guarantees  the  existence  of  a  unique  solution. 
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8.8  Chapter  summary 

A  theorem  for  existence  and  uniqueness  of  solutions  to  a  broad  class  of  abstract 
Cauchy  problems  has  been  presented,  along  with  its  proof.  Also,  an  application  has  been 
described,  formulated,  and  shown  to  satisfy  the  hypotheses.  The  existence  of  a  unique 
solution  is  guaranteed. 
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IV.  A  nonlinear  damping  term 


A  =  A(u)  = 


\ 


) 


(158) 


Consider  utt  +  D2(P(u)D2ut)  +  aD4u  =  0.  Formulate  this  as  an  abstract  system  as 
before  with 

/ 

0  -1 
aD4  D2(P(u)D2) 

where,  in  the  abstract  system,  u  is  a  vector  with  components  ui  and  u 2.  The  space  X  and 
D(A)  are  the  same  as  described  in  the  previous  chapter. 

It  is  convenient  and  informative  to  treat  a  specific  case  and  demonstrate  how  the 
conditions  are  verified. 

Consider  utt  +  D2((Po  +  P\u)D2ut)  +  aD4u  =  0  where  a,  Po,  and  Pi  are  constants. 
Formulate  this  as  an  abstract  system  as  before  with 


As  =  A(u)  =-. 


0  -1 
aD4  D2((p0  +  PiUi)D2) 


\ 


(159) 


and  D(As)  =  D(A3).  Certain  restrictions  will  be  placed  on  Po  in  terms  of  other  constants 
in  the  problem.  Portions  of  the  problem  formulation  which  are  different  from  before  are 
described  next. 


4-1  Preparation  for  application  of  the  existence  theorem 


Let  Y  be  the  set  of  points  (pairs  of  functions)  in  D(As)  endowed  with  the  norm 


ll*  =  WU  +  IP"n|[tl  +  ||£>J»|k. 


(160) 


Notice  that  does  not  depend  on  i.  The  linear  space  Y  is  complete  as  shown  before. 
The  graph  norm  has  not  been  used  this  time.  Use  of  the  graph  norm,  according  to  the 
previous  pattern,  would  require  W  to  be  closed  and  bounded.  The  current  strategy  uses 
the  norm  (160). 
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Some  earlier  terminology  is  now  reviewed.  Let  W  =  B(uq]  R)  be  an  open  ball  in 
Y  with  center  «o  (the  initial  condition  vector)  and  radius  R  (not  yet  specified).  Let 
M  =  |K||y- 

Several  lemmas  are  now  presented.  They  will  be  used  to  identify  the  restrictions  on 
0o  and  in  other  computations.  The  purpose  of  these  lemmas  is  to  identify  bounds  on  y2" 
and  y2 in  terms  of  D4y2. 

Lemma  64  For  any  c  G  [0, 1]  and  y2  G  H 3,  ||2/2"||oo  <  |y2W(c)|  +  ||2/2W||l2- 

Proof:  This  follows  immediately  from  the  same  line  of  argument  as  in  the  proof  of  Lemma 

12.  □ 

Lemma  65  For  any  y2  G  H3  with  y2(0)  =  2/2,(l)  =  0  there  is  some  c  G  [0, 1],  such  that 
y2n(c)  =  0. 

Proof:  Note  the  continuity  of  y2"  as  an  element  of  H 3  and  the  boundary  conditions 
2/2'(0)  “  ya'(l)  =  0.  The  existence  of  the  desired  c  G  [0,1],  such  that  y2"(c)  =  0  is 
immediate  from  Rolle’s  Theorem.  □ 

It  is  useful  to  note  that  there  are  at  least  two  such  values,  say  c\  and  c2.  This  is 
because  an  application  of  Rolle’s  Theorem  to  the  continuous  function  y2  and  the  boundary 
conditions  on  y2  gives  a  point  c  G  (0, 1)  such  that  y2'(c)  =  0.  Then  the  argument  of  Lemma 
65  can  be  applied  to  each  of  the  intervals  (0,c)  and  (c,  1). 

Corollary  66  For  y2  G  H3  and  satisfying  the  boundary  conditions, 

Moo  <11^1^.  (161) 

Lemma  67  For  any  c  G  [0,  l]  and  y2  G  H4, 

Mloo  <1^)1+11^211^.  (162) 

Proof:  This  follows  immediately  from  the  same  line  of  argument  as  in  the  proof  of  Lemma 

12.  □ 
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But,  again  by  Rolle’s  theorem,  there  is  some  c  £  [ci,^]  such  that  1/2 w(c)  =  0.  Thus, 
the  following  corollary  is  immediate. 

Corollary  68  For  t/2  G  H4  and  satisfying  the  boundary  conditions, 

ll»"IU  <  ll^wlU,.  (W3) 

Lemma  69  For  w  £  B(uo\R )  =  W ,  |rt>i |  <  where  M  =  ||tto||y  and  a  is  the 

constant  in  (159)  which  is  also  used  in  the  definition  of  ||  •  ||y . 

Proof:  Since  w  £  W,  it  is  clear  that  ||w||y  <  R  +  M.  Recall  that  components  of  elements 
of  Y  satisfy  the  boundary  conditions  in  the  specification  of  D(A).  From  the  definitions  it 
is  not  difficult  to  see  that 


M 


—  IMloo 

<  ||Wl"||i, 

1 

“  a»/2 

< 

< 

< 


al/2 
1 

M  +  R 
0(1/2 


INI* 


Mv 


(164) 


The  restriction  that  fio  will  be  required  to  satisfy  is  now  given.  It  is  required  that, 
for  every  x  £  [0, 1],  +  >  0  for  all  w  £  W.  Since  ||wi||oo  <  this  can  be  satisfied 

with  a  finite  choice  for  fio-  Or,  more  to  the  point,  it  is  desirable  that  >  -wi(x)  for  all 
x  £  [0,  l].  Note  that  the  choice  of  small  values  for  R  and  M  will  allow  more  flexibility  in 
the  choice  of  acceptable  fi0  and 


J,.'2  Applying  Theorem  J,S 
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Consider  whether  the  nonlinear  damping  problem  satisfies  the  hypotheses  of  Theorem 
48.  Hypotheses  1  and  2  are  clearly  satisfied.  Hypothesis  3  is  satisfied  with  Q  =  I  +  A  as 
before. 

Bounds  on  Q  and  (Q-1  for  Hypothesis  4  may  be  obtained  with  some  effort.  In  order 
to  obtain  such  bounds  most  easily,  several  elementary  lemmas  are  presented. 

Lemma  70  For  a,  b  >  0,  (a2  +  i2)1^2  <  a  +  b. 

Proof:  See  Appendix  D.  □ 

Lemma  7x  For  any  real  numbers  a  and  b,  (a  +  6)2  <  2 (a2  +  62). 

Proof:  See  Appendix  D.  □ 

Corollary  72  For  positive  real  numbers  a  and  b,  2 _1/2(a1/2  -f  61/2)  <  (a  +  6)1/2. 

Lemma  73  For  positive  real  numbers  a,b,  and  c  it  holds  that  (a+6+c)2  <  3(a2  +  62  +  c2). 


Proof:  See  Appendix  D.  □ 

Lemma  74  For  any  positive,  real  a,  b,  andc,  it  holds  that  (a+b+c)1/2  <  (a1/2-}-!)1/2-}-^1/2). 

Proof:  See  Appendix  D.  □ 

Consider  a  bound  for  Q.  Equation  (40)  will  be  used  several  times  in  the  following 
string  of  inequalities. 


||<2(ti>)||y_Of  =  \\T  +  As(w)||y-or 
il(J  +  A8(w))y\\x 

oSfi'  II* 

oim  II  y  +  M™)y\\x 

HtflU  +  \\d4Vi\\l2  + 
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<  _ Mix  + 

“  oJ”&IMI*  +  ll^4»l||La+IP4»llfa 


0 lliiiu  +  PVP  +  II D4»IU,  t  Ml;  +  IP^iP  +  P*»IU, 


llsllx 


+  sup 


MsMylU 


<  1 +  sup 


p8(n>)y|U 


o^yMx  +  WD^  +  WD^Wl, 


<  1 +  sup 


PsMi/lk 


P4J.  life  +  IP4!/2llt, 


=  1 +  sup 
0  PjZY 


(a  fo  (y2")2  dx  +  fj  (otD4y1  +  p0D4y2  +  PiD2  (toa D2y2))2  dx ) 


1/2 


||-D42/i||l2  +  |P42/2||l2 


<*1/2  (fo  W')2  dx )1/2 

-  ‘^'IrPSCTPa 

( /o  (aD4yi  +  +  A-D2  (wi&yi))2  dx)'1* 

+  oSfy  PVIk  +  PVh 


< 


(/o1  {aD\  +  PoD4y2  +  pxD2  (w^y,))2  dx)1/2 
1  +  “  +  &r - »1 - 


<  1  +  a1/2  +  21/2  sup 


<  1  +  a1/2  +  21/2  sup 


(fo  (<xD4yi)2  dx  +  /o  (p0D4y2  +  PiD2  (wxD2y2))2  dx 


\2  j_\i/a 

wiwy2)x 

WnWi*  +  IP4^I|l2 

\2  J  \l/2  .  /  /*1  /  i-v/4  ^  /  *%o  \\  2  *  \l/2 


<  1  +  a1/2  }-  21,?a  +  21/2  sup 

ofyeY 


(fo  MV)2  dx)  +  (Jo1  (ft,  fl  V  +  PiD2  (w1D2y2))2  dx) 

ll-DVIV  +  IPVIV 

21/2  (/J  (PoD4y2)2  dx  +  $  (pxD2  (wxD2y2))2  dx) 
\\D4yi\\L2  +  \\D4y2\\L2 


1/2 


< 


1  +  a1/2  +  2l/2a  +  2 


SUD  WPqD^Wl* 

ofrer  \\D4yi\\L2  +  ||-D42/2||l2 
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( fo  {PiD2  {w\D2y2))2  da;) 1/2 
+2oSSy  ||-£>42/i!|z,2  +  !I-C>42/2||z,2 


=  l  +  a1/2  +  21/2a  +  2|/?0|  +  2  sup 

<¥veY 


\Pi\  (fo  {wiD4y2  +  2wi'D3y2  +  w\D2y2)2  dxf12 
\\DAn\\h  +  \\D*V2\\i* 


<  1  +  a1/2  +  21/2a  +  2\p0\ 

31/2|/?il(/o  Wi2(D4y2)2  d *  +  2/q1  (wi1)2 (D3y2)2  dx  +fx(w1")2(D2y2)2  dxf12 

+2oSfy  IP42/iIU2  +  IP42/2||l2 


<  1  +  a1/2  +  21/2a  +  2|Aj| 

(. R  +  M )  (Jo1  (^2)2  dx  +  2  /J  (D32/2)2  dx  +  fj  (P2^)2  dx) 

+2 . 3./2|A|4upy - i - VKJL  +  VUh. - 


1/2 


<  1  +  a1/2  +  21/2a  +  2j/30| 

+2  •  31/,2|/?j|(i2  +  M)  sup  2^fr+}}D^ 

0  \\D4yi\\l2  +  IP ‘2/2 1, ^ 


<  l  +  a1/2  +  21/2a  +  2|/30|  +  8-31/2|/?1p  +  M)  (165) 

So,  choose  Aq  =  1  +  a1/2  +  2l^,a  +  2|/?o|  +  8  •  Zll2\^i\{R  +  M). 

Now  that  Q  is  bounded,  recall  that  it  is  also  linear,  one-to-one,  and  onto.  Then,  by 
a  standard  corollary  (see  [22:pg  47],  or  [73:pg  77])  of  the  Open  Mapping  Theorem,  Q~x  is 
also  bounded,  say  by  Xq. 

A  suitable  choice  for  hq  is  ha  which  is  determined  in  (168). 

Hypotheses  5  and  6  are  trivially  satisfied  since  the  variable  norms  are  not  used  in 
this  application. 

The  first  part  of  Hypothesis  7  is  satisfied  for  each  w  €  W  as  argued  in  the  (3(x ) 
case.  For  clarity,  the  dissipative  argument,  which  appears  on  the  surface  to  be  different,  is 
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presented.  First,  some  notation  is  clarified. 


^8  = 


As(u)u  = 


/ 

\ 

/ 

V 

/ 


\  /„  \ 


As(u>)u  = 


0  -1 

aD4  i?2((/?o  +  /3itii)D2) 

0  -1 

oD4  D2((/Jo  +  /3i«i)D2) 

-«2 

aD4ui  +  D2  ((/5o  +  /?i«i)jD2«2) 
-«a 

aD4ui  +  jD2  ((/?0  +  PiWi)D2u2) 


* 


(A8{w)v,u) 


-a  J  u2"u\"  dx  +  J  u2aD4uidx  +  J  u2D2  ((/3o  +  P\’W\)D2u2j  dx 

u2D  ((/Jo  +  /?iWi)D2u2)  \o  ~  Jq  Du^D  ((A>  +  P\W\)D2u2j  dx 

-Du2  ((/?o  +  PiW\)D2u2j  |o  +  J  D2u2(P0  +  PiWi )D2u2  dx 

[  (Po  +  PiWi)(D2u2)2  dx  (166) 

Jo 


where,  in  the  line  with  three  integrals,  the  first  two  cancel  out  after  integration  by  parts. 
Also,  the  boundary  terms,  arising  from  integrations  by  parts,  are  all  zero.  If  it  happens  that 
the  quantity  Po  +  P\W\  is  greater  than  or  equal  to  zero  for  all  x  (which  is  required  above) 
then  -  A8{w)  is  dissipative.  The  remainder  of  the  argument  to  establish  A$(w)  £  G(X,  1, 0) 
is  the  same  as  for  the  p(x)  case.  This  completes  the  discussion  of  the  first  part  of  Hypothesis 
7. 


For  the  second  part  of  Hypothesis  7  note  that,  from  the  work  for  Xq  leading  to  (165), 
it  follows  that 

Xa  =  <*1/2  +  21/2a  +  2|A>|  +  8  •  31/2! px\{R  +  M)  (167) 
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is  suitable. 


* 


Next,  a  suitable  ha  must  be  identified  for  the  third  part  of  Hypothesis  7.  Notice  that 
in  the  current  application  there  is  no  explicit  time  dependence  in  .4g.  For  the  first  inequal¬ 
ity,  in  the  following  sequence  of  equalities  and  inequalities,  a  comment  is  appropriate.  The 
first  two  terms  in  the  integrand  are  bounded  in  terms  of  the  supremum  norm  on  the  factoi 
involving  w\.  The  third  term  in  the  integrand  is  bounded  in  terms  of  the  supremum  norm 
on  the  D2yi  factor  which  leads  to  an  Li  bound  on  DAyi-  Now,  by  definition 


-  ^s(w)||r-A' 


IPs(u>)  -  Mw)Mx 
o*y£y  IMb' 


II (Mw)  -  i4s(w))ylU’ 

o#y4V  \\y\\x  +  ll^'yilUa  +  \\d4V2\\l2 


=  sup 


sup 
0  £yeY 


0  0 
0  PiD2  ((wx  -  wi)D2) 


V  r 

y2 


IMU  +  \\D4yi\\Li  +  \\DAy2\\L2 


X 


sup 


01 D2  {{wi  -  tdi)D2y2) 


Wv\\x  +  \\D*VilL»  +  \\DAy2\\Li 


01  ( So  (^2  ((wi  -  wi)D2y2))2  dx)1/2 
oJIlv  h\\x  +  IPVIU2  + 


=  IAI  sup 
0  *y<:Y 


( So  ((i»i  ~  t»i )D4y2  +  2(wi  -  w\)'Dzy2  +  (t»i  -  wi )"D2y2)2  dxf12 

MU  +  WD4y\\\L2  + 


<  \0x\  sup 


otueY  \\D4y2\\L2 

,  -  tf.")ll| ,+  |K’ -  tf.ll,)  IP'fclli,) 

<  3  |y#i|  sup  - iFnirii - 

0#yer  \\D  y2\\L2 

=  31/2|/?1|61/2|K,-uJ1"|k2 


1/2 


1/2 


=  3  •  2l!2\0\  \  ||w  -  w||a' 


(16*) 
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so  it  is  acceptable  to  choose  [ia  =  3  •  2 1^2|/?i|. 

Hypothesis  8  follows  immediately,  as  before,  since  Q-1  is  the  resolvent  of  /ls. 

Identification  of  an  acceptable  T  is  a  several  step  process.  Since  this  process  has 
been  demonstrated  in  detail  for  a  previous  application,  it  will  not  be  repeated  here. 

The  first  step  of  the  solution  algorithm  is  now  described.  Apply  the  arguments  of  the 
0{x)  case  to  /ls(«o)-  This  yields  a  solution  u(t )  based  on  the  single  operator  .4s(t/o).  Next 
apply  the  iterative  scheme,  as  in  the  0(t,x)  case,  to  the  operator  family  .4g(w(/)).  This 
gives  an  evolution  system  U(t,s )  which  generates  a  solution  u(t).  This  is.  .just  one  slop  of 
the  iteration.  The  fixed  point  algorithm  (from  the  proof  of  Theorem  -IS)  guaiantocs  that 
the  iterates  will  converge. 

The  continuous  dependence  result,  Theorem  3.2  of  (35:pp  170-171],  is  applicable  to 
this  problem. 


./.5  .4  numerical  example 

In  this  section  a  specific  example  is  presented.  The  equation  is  the  same  as  earlier  in 
the  chapter,  but  now  the  constants  take  on  specific  values.  Also,  a  specific  initial  condition 
vector  is  given.  The  example  and  the  software  to  propagate  its  solution  aie  from  [60]. 

Choose  the  following  values  for  the  constants. 


a  =  .008 

0o  =  -01 

0i  =  .001  (169) 


For  the  first  component  of  vq  choose 


si  nil  u  -  sm//  sinli// -  sm//  ,  .  ,  . 

til  =  Sill  VX  +  - : —  COS  VX - ; —  COSll  V  -  Slllll  V.  ( 1  /  0 ) 

cos  v  -  cosh  //  cos  //  -  cosh  v 


This  is  shown  in  Figure  4  with  /'  =  4.73.  The  choice  of  //  satisfies  cof-(u)co.'-li(n)  =  J.  This 
is  required  for  wq  to  satisfy  the  boundary  conditions  necessary  for  it  to  In*  in  D(A$). 


.  $9 


0.600 


Choose  0  for  the  second  component  of  wo¬ 
lf  fli  were  zero,  then  this  problem  could  be  solved  by  separation  of  variables.  The 
fundamental  frequency  corresponds,  in  this  case,  to  i /  =  4.73  approximately.  Choose  this 
value  for  v  and  test  the  basic  algorithm  for  convergence. 

The  solution  algorithm  is  outlined  now. 

1.  Establish  values  for  the  constants  in  the  problem. 

2.  Establish  a  grid  and  assign  an  initial  estimate  to  the  solution  at  each  grid  point.  As 
the  first  estimate  suppose  that  the  initial  condition  vector  is  a  constant  solution. 

3.  Initialize  a  counter  for  time  increments,  say  j  =  0. 

4.  Evaluate  the  quantity  30  +  /3ju  for  each  grid  point,  based  on  the  cuirent  estimate  for 

u. 

5.  Increment  the  counter:  j  =  j  +  1. 

(5.  Propogate  the  approximate  solution  from  /  =  [j  -  1)A/  to  t  =  ,/A/.  (This  step  is 
done  with  the  program  DGEAR  from  the  standard  Fortran  package  known  as  IMSL.) 

7.  If  the  desired  final  time  has  not  been  reached,  go  to  step  5. 

8.  If  the  solution  lias  not  converged,  go  to  step  3. 

9.  This  completes  the  algorithm. 

The  results  are  presented  in  Figure  5.  The  plot  is  for  the  vertical  disp.accment  of 
the  midpoint  of  the  beam.  The  line  across  the  top  represents  the  initial  estimate,  which 
is  the  initial  condition  as  a  constant  solution.  The  next  curve  down  is  the  estimate  after 
one  iteration  of  the  algorithm.  The  third  curve  is  the  estimate  after  two  iterations.  The 
third  iteration  lies  on  top  of  the  second  one  and  cannot  be  distinguished,  though  it  has 
been  included. 

A  second  example,  with  0\  =  .008.  is  shown  in  Figure  6.  Notice  dial  while  the 
spacing  is  different,  the  basic  character  is  the  same. 
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4-4  Chapter  summary 

The  nonlinearly  damped  beam  vibration  problem  has  been  shown  to  satisfy  the 
hypotheses  of  the  theory  developed  in  the  previous  chapter.  An  example  has  been  presented 
which  demonstrates  rapid  convergence  of  the  algorithm  developed. 
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V.  Conclusions  and  Recommendations 


This  chapter  provides  a  brief  summary  of  results  which  have  been  demonstrated.  It 
also  includes  recommendations  for  further  study. 

5.1  Conclusions 

A  standard  mathematical  model  for  approximating  the  transverse  vibration  of  a  beam 
has  been  generalized.  The  existence  and  uniqueness  of  solutions  have  been  established. 
Certain  continuous  dependence  results  have  been  presented. 

A  nonlinear  problem  which  was  not  previously  known  to  have  a  solution,  which  could 
be  obtained  by  convergent  sequences  of  approximations,  has  now  been  shown  to  have  such 
a  solution. 

5.2  Future  work 

Several  extensions  are  quite  logical  for  the  work  contained  herein.  Certainly  it  would 
be  desirable  to  extend  the  work  to  more  general  boundary  conditions.  This,  howevet,  may 
not  be  at  all  straightforward.  See  Appendix  B. 

Use  of  the  equations  in  a  parameter  identification  scheme  is  certainly  appropriate. 
Indeed,  it  was  with  such  use  in  mind  that  this  project  was  undertaken.  This  document 
provides  the  theoretical  basis  to  undertake  a  nonlinear  identification  scheme  along  the  lines 
of  the  linear  equation  based  schemes  in  [12].  The  implementation  of  mimeiical  methods 
for  this  problem  should  be  very  interesting. 

Of  course,  issues  relating  to  the  stability,  long  term  behavior,  stability  of  iletative 
schemes,  and  continuous  dependence  on  other  elements  of  the  equation  aie  all  of  intciest. 
This  would  include  time-dependent  Trotter-Katr  results,  such  as  in  [57],  [48:pp  17,  19], 
[32],  and  [41].  Nonlinear  Trotter-Kato  type  results  may  be  found  in  several  references.  See 
[47:pp  469-476],  [41],  [50:pp  223-224],  [51:pp  403-404],  [52:pp  24-25],  and  [3-1]  for  some  of 
the  early  work.  A  good  summary  is  in  [IS].  A  product  formula  version  is  addicsscd  in  [62] 
and  [61]. 
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Application  of  the  style  of  analysis  presented  in  tlr  work  to  other  equations  is 
also  of  interest.  For  example,  if  the  original  equation  had  included  rotary  inertia,  then 
the  nonlinear  version  of  the  equation  would  have  beer,  different.  A  separate  analysis  is 
necessary.  Allowing  time  dependence  of  the  fundamental  parameters  in  the  derivation  is 
another  variation  that  would  be  interesting  to  pursue. 
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Appendix  A.  Derivation  of  the  Euler-Bernoulli  model 


The  Euler-Bernoulli  model  is  derived  as  follows.  The  total  kinetic  energy  of  the  beam 
is  assumed  to  be 

r=W.  m(x)vf(to,x)dx.  (171) 

The  elastic  potential  energy  due  to  bending  (assumed  to  be  the  total  potential  energy 
since  changes  in  gravitational  potential  will  be  ignored)  is  assumed  to  be 


V  =  1/2  J  EI(x)ulx(to,x)dx. 


(172) 


For  convenience,  the  product  EI{ x)  is  treated  as  a  single  entity  in  the  equations. 

The  total  mechanical  energy  in  the  system  for  the  Euler-Bernoulli  model  is  assumed 


to  be 


S  =  T- \-V  ss  1/2  J  [ro(x)ti?(t0,  *)  +  EI(x)ulx(t0,  x)]da 
and  the  equations  of  motion  are  obtained  by  setting 


(173) 


dS/dt  =  0. 


(174) 


Consider  the  expression  for  the  time  rate  of  change  of  the  total  mechanical  energy. 
The  arguments  of  u  are  suppressed  for  convenience.  Leibnitz  rule  (see  [72:pg  5]  or  [69:pp 
163,170])  is  applied  to  take  the  differentiation  inside  the  integration. 


d£ 

dt 


=  -  -lL 

dt  \  2  Jo 


,  .  fdu\ 2  .  .  / d2u\ 

”WU)  +B7W(w) 


dx 


If  dx  +  l  f  &!(*)§ dx 

‘  Gt,  \Gt>  /  6  j o  Gv  \G+'  J 


&  J  0 


-  !L  —  fLpT(r\d2u  ^  Ar 

~  Jo  O  dt  dt 2  d  +  Jo  EI^  dx 2  dtdx 2  d 


(175) 
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Now  expand  the  second  integral  in  the  last  expression  using  integration  by  parts.  (It  is 
assumed  that  u  is  smooth  enough  to  justify  changes  in  the  order  of  integration.) 


d2u  dzu 

d^didi *dx 


EI{x) 


d2u  d2u 
dx i  dtdx 


L 

o 


rL  d2u  d 
Jo  dtdx  dx 


dx 


=  EI(x) 


d2u  d2u 
dx 2  dtdx 


L 

o 


d_ 

dx 


du 

dt 


L 

0 


Hence, 


(176) 


d_ 

dx 


(177) 


It  is  not  difficult  to  see  that  for  clamped,  pinned,  or  free  boundaries  the  boundary  terms 
in  this  expression  are  zero.  So,  the  equations  of  motion  must  come  from 


,  ,d2u  d2 

ro(l^+w 


dx  —  0. 


(178) 


The  quantity  ||  will  not  be  zero  over  any  interval  for  any  interesting  beam.  In  fact, 
as  time  goes  by,  will  take  on  a  wide  variety  of  values.  The  only  way  to  guarantee 
d£/dt  =  0,  is  to  require 


roW§?+&('B/(l)5?)=0  (179) 

for  all  x  €  [0,  L)  and  t  >  0.  This  is  all  made  precise  in  the  calculus  of  variations.  The 
terminology  admissible  function  is  used  for  possible  values  of  |^.  See,  for  example  [72:pp 
217-220]. 
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In  the  case  of  constant  density,  constant  cross  section,  and  constant  modulus  of 
elasticity,  the  equation  reduces  to 


82u  „r84u  n 

mdfi  +EISx*  _0' 


(180) 


This  is  usually  written  as 


utt  +  ( auxx)xx  =  0  <  x  <  L]  0  <  t  (181) 

where,  consistent  with  the  physical  situation,  it  is  assumed  that  a  >  0.  Since  a  is  a 
constant  the  parentheses  are  not  really  necessary  but  they  are  suggestive  of  a  more  general 
equation. 
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Appendix  B.  No  simple  extension  to  other  boundarxj  conditions 


Consider  an  example  which  demonstrates  that  the  type  of  analysis  done  herein  will 
not  carry  over  to  spaces  similar  to  X  but  with  different  boundary  conditions.  In  particular, 
if  the  boundary  conditions  on  the  first  component  required  a  zero  value  on  the  second 
derivative  but  not  on  the  first,  then  the  space  would  not  be  complete. 

Consider  the  sequence  of  functions 


*«(*)  = 


x  0  <  x  < 

"  T— 2"  ®  <X  <  22^T1 

1-®  1 


(182) 


The  limit  of  this  sequence  is  given  by 


z{x)  =  lim  zn(x )  =  < 


0 


x  =  0 


-ar  +  x  0  <  x  <  1 


0 


x  =  1 


(183) 


It  is  easy  to  see  that  each  element  of  the  sequence  has  second  derivative  zero  at  the 
end  points  but  the  limit  does  not.  Slight  modifications  of  this  example  demonstrate  that 
boundary  conditions  other  than  those  cited  originally  do  not  lead  to  a  complete  space  in 
the  given  norm. 


One  possible  approach  to  get  the  second  derivative  under  control  would  be  to  base 
X  on  II3  x  H°  instead  of  if2  x  H°.  But,  when  a  sequence  is  given  as  Cauchy  in  this  new 
X  (with  the  old  norm),  it  is  not  possible  to  show  that  its  first  components  form  a  Cauchy 
sequence  in  H3.  So,  there  is  no  simple  answer  here. 


Another  approach  is  to  let  X  be  the  completion  of  the  subset  of  II'2  X  JJ°  determined 
by  the  boundary  conditions.  This  approach  has  not  been  fully  analyzed. 


Another  difficulty  arising  from  other  boundary  conditions  is  in  establishing  the  dense¬ 
ness  of  D(A)  in  X.  In  particular,  Theorem  7  and  its  generalization  in  Appendix  D  will 
not  apply.  But,  see  [I2:pp  16-18, 42ff]  and  [58:pg  8,  2.5]. 
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Appendix  C.  Not  the  generator  of  an  analytic  semigroup 

While  X  has  been  defined  as  a  real  Banach  space,  complex  valued  functions  could 
be  allowed  and  the  inner  product  modified  (with  conjugation  of  the  second  factor  in  the 
integrand)  to  make  everything  proper  for  the  possibility  of  extension  to  an  analytic  semi¬ 
group. 

Theorem  75  The  operator  —A\  given  by  (26)  is  not  the  generator  of  an  analytic  semi¬ 
group. 

Proof:  The  proof  of  this  theorem  will  be  given  after  several  preliminaries  have  been 
established.  □ 

The  key  to  the  argument  is  Theorem  2.5.2.a,c  of  [58:pp  61-63].  For  convenience  the 
relevant  portion  of  the  theorem  is  stated  next.  Its  proof  is  in  the  reference  as  cited. 

Theorem  76  Let  S(t)  be  a  uniformly  bounded  Co  semigroup.  Let  A  be  the  infinitesimal 
generator  of  S(t)  and  assume  0  G  p(A).  The  following  statements  are  equivalent: 

a  S(t)  can  be  extended  to  an  analytic  semigroup  in  a  sector  A  $  =  {z  :|  argz  |<  5}  and 
||S'(z)||  is  uniformly  bounded  in  every  closed  subsector  As>,  S'  <  5,  of  As- 

c  There  exist  0  <  6  <  7r/2  and  M  >  0  such  that 

p(A)  D  S  =  |a  :|  arg  A  |<  ^  -f  6  j  U  {0}  (184) 

and 

M 

||IE(A  :  A)||  <  —for  A  e  E,  A  ?  0.  (185) 

rrooi:  oee  joo:pp  6j.-6oj. 

From  Theorem  76  it  is  clear  that  the  resolvent  of -A  must  include  the  entire  imagi¬ 
nary  axis  if  -  A  is  to  be  the  generator  of  an  analytic  semigroup.  But,  as  will  be  established 
for  -Ai  (see  (26)),  the  point  spectrum  includes  infinitely  many  well  spaced  points  on  the 
imaginary  axis.  Hence  the  conclusion. 
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Lemma  77  The  point  spectrum  of  —A\  includes  arbitrarily  large  values  on  the  positive 
and  negative  imaginary  axis. 

Proof:  First  recall  that  -A\  is  skew-adjoint.  This  leads  to  the  necessity  of  any  eigenvalue 
being  on  the  imaginary  axis.  In  particular,  if  —A\u  -  Xu  then  (~Aiu,u)  =  (Ati,u).  But, 
as  was  shown  earlier,  -A\  is  dissipative  and  hence  Re{-A\u,u)  =  -Re(A\U,u)  <  0. 
From  skew-adjointness,  (-A\u,u)  =  (u,(-Ai)*u)  =  (u,.4i ti).  Thus  Re(Aiu,u)  <  0. 
Hence,  Re(A\U,u)  =  0.  On  the  other  hand,  j?e(Ati,u)  =  Re(X(u,u))  =  jRe(A||u||£2)  = 
jRe(A)||u||£  ,  which  must  be  zero  for  arbitrary  u.  Hence  Re( A)  =  0. 

This  has  only  established  that  any  eigenvalues  for  -A\  are  constrained  to  the  imag¬ 
inary  axis,  it  remains  to  show  that  there  are  eigenvalues  there.  This  will  be  done  in  a  very 
direct  constructive  fashion.  Solve 


for  A.  This  is  the  same  as  the  system 

u-i  =  A«i 

-aDAu  i  =  Xu2-  (187) 

This  can  be  written,  by  simple  substitution,  as 

-  aDAUi  =  X2ui  (188) 
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which  can  be  written 

ft2 

r4  -  —  =  0,  (191) 

a 

it  is  clear  that  the  roots  are 

r  =  itfWa-V*.  (192) 

For  convenience,  let  7  =  /31/2a-1/4.  Then 

u  =  cie71  +  C2e~'ix  +  C3COS7S  +  C4sin7a: 

u'  =  7Cie'YX  -  7C2e-72:  -  7C3sin7«  +  7C4COS7®  (193) 

Now  use  boundary  values  to  solve  for  the  coefficients. 

u(0)  =  0  =$>  Cl  +  C2  +  C3  =  0 

ti'(0)  =  0  =*►  Ci  -  C2  +  C4  =  0 

u(l)  =  0  =i>  e yC'  +  e-7C2  +  C3 COS7  +  c4 sin7  =  0 

u'(l)  =  0  =$>  e7ci  —  e-7C2  -  C3sin7  +  C4COS7  =  0  (194) 

The  important  issue  here  is  to  find  7  so 
requires 

1  1 
1  -1 
e7  e"7 
e7  -e"7 

This  is  reduced  as  follows. 

1 
0 

0  -»,7  +  e-7  -e7  +  cos7  sin7 
0  -e7  -  e-7  -e7-sin7  COS7 


that  this  system  has  a  nontrivial  solution.  This 


1  0 

0  1 

cos  7  sin  7 
-sin 7  cos 7 


=  0 


(195) 


-2 


1 

-1 


=  0 
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-2 


-1 


1 


-e7  +  e-7  -e7  +  cos7  sin7  =  0 
-e7  —  e-7  -e7  -  sin7  cos7 
— 2[— e7  cos  7  +  cos2  7  +  e7  sin  7  +  sin2  7] 

+[-e7  cos  7  +  e-7  cos  7  +  e7  sin  7  +  e-7  sin  7] 

+[e27  +  e7sin7  -  1  -  e-7sin7  -  c27  +  e7  C0S7  -  1  +  e-7  COS7]  =  0 


-4  +  2e7  cos  7  +  2e-7  cos  7  =  0 
e7  4.  p — If 

-1  + - - - cos  7  =  0 

cosh  7  cos  7  =  1 


cosh  7  =  sec  7  (196) 


The  problem  has  now  been  reduced  to  finding  real  values  of  7  such  that  cosh  7  =  sec  7.  But 
cosh  7  is  defined  for  all  real  7  and  is  always  greater  than  or  equal  to  1.  On  the  other  hand 
sec 7  has  asymptotes  at  7  =  for  n  =  •  •  •  —  2,  — 1, 0, 1,2,  •  •  • .  Except  for  n  =  0,  each 

asymptote  is  approached  by  a  unique  branch  of  sec7  which  intersects  cosh 7.  (For  n  =  0, 
two  branches  intersect  at  the  same  point.)  Thus  an  infinite  number  of  arbitrarily  large 
positive  and  negative  values  of  7  are  obtained.  Now  recall  7  =  f31^2a~1/4,  so  (3  =  72a1/2 
where  a  >  0  is  a  constant.  It  is  clear  that  (3  can  be  arbitrarily  large  and  hence  so  are  the 
A  values. 

This  establishes  the  lemma.  □ 

Inconveniently,  zero  is  in  the  spectrum  of  -A  and  hence  0  £  p(-A).  Hence,  an 
adjustment  is  necessary  to  apply  Theorem  76.  Since  -A  is  the  generator  of  a  Co  semigroup 
its  resolvent  set  contains  the  positive  real  axis.  A  small  shift  of  the  problem  will  put  0  in 
the  resolvent  set. 


Lemma  78  The  point  spectrum  of  -A  -  el  includes  arbitrarily  large  values  on  the  line 
Re(  A)  =  -e. 

Proof:  The  argument  is  exactly  the  same  as  in  the  previous  lemma  with  A  replaced  by 

A  +  e.  □ 
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Lemma  79  -A  -  d  is  not  the  generator  of  an  analytic  semigroup. 


Proof:  By  Theorem  6.1  of  [19:pg  38],  -A  -  d  is  the  generator  of  a  Cq  semigroup  of 
contractions.  And  hence  the  generator  of  a  uniformly  bounded  semigroup.  (Or,  see  [30:pg 
499].)  To  justify  application  of  this  theorem,  note  that  -d  is  a  dissipative  operator  for 
any  e  >  0.  In  particular 


( du,u )  =  -e(u,u)  - 

=  -«2 

<  0  (197) 

There  is  one  more  preliminary  before  the  proof  of  the  lemma.  It  must  be  established  that 
0  e  p(-A-el).  That  is,  show  that  A  +  d  has  a  bounded  inverse.  But,  e6/)(-A)-^6  +  A 
has  a  bounded  inverse. 

Now  Theorem  76  applies.  Furthermore,  any  wedge  will  cross  the  line  x  =  -e  and 
hence  the  sector  will  include  elements  of  the  point  spectrum.  Thus,  the  wedge  does  not  lie 
intirely  in  the  resolvent  set  and  the  lemma  is  established.  □ 

This  concludes  the  preliminaries.  The  proof  of  Theorem  75  is  structured  as  a  proof 
by  contradiction. 

Proof:  Suppose  -A  is  the  generator  of  an  analytic  semigroup.  Then  by  Corollary 
3.2.2  of  [58:pg  81],  -A-d  is  the  generator  of  an  analytic  semigroup.  This  is  a  contradiction 
to  the  lemma  just  established. 

The  theorem  follows.  □ 
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Appendix  D.  Miscellaneous  proofs 


D.l  Proof  of  Lemma  7 

Lemma:  If  y  G  H2  and  y(0)  =  y(l)  =  y'( 0)  =  y'(  1)  =  0,  then  y  G  Hq. 

Proof:  When  y  G  J72,  it  follows  that  y>y\y"  G  L 2.  Since  y"  G  X2,  there  is  some 
sequence  {y'f}  C  Cq°  such  that  y"  -4  y". 

Let  y„(a:)  =  /q  y([(®)  di.  Note  that 


<  [  bn(*)l  d® 

Jo 

<  [  ll/»(*)l  d® 

Jo 

<  /  (j/n(*))2d®  (198) 

Jo 

which  is  bounded  since  y"  G  £2-  (Note  that  the  last  inequality  need  not  hold  for  y”  <  1, 
but  in  this  case  boundedness  is  obvious.)  Hence  y'n  G  L^. 

Claim:  y'n  -4  y' 


y'n(x)  j  y"(x)dx  =  y\0)  +  y\x)  =  y'(x). 
Jo 


(199) 


Therefore,  the  claim  holds. 

Similarly,  let  yn(®)  =  foVn(x)^x>  and  no*;e  that  Vn  £  -£2*  It  is  easy  to  see  that 
yn  G  H$.  If  yn  -4  y,  then  y  G  Hq  by  the  completeness  of  H$.  But  this  follows  by  the  same 
line  of  argument  used  to  show  y'n  -»  y'.  □ 

A  more  general  version  of  this  lemma  is  also  known,  as  stated  below. 

Lemma:  Let  u  G  Bm  be  given.  Then  u  G  //q1  if  and  only  if  u^( 0)  =  uM(l)  =  0  for 
every  k  <m  -  1. 

Proof:  This  is  essentially  a  corollary  to  Theorem  3.3  of  [56:pg  67].  It  is  taken  from 
[56:pg  91].  The  proof  requires  several  other  theorems  which  are  also  in  the  cited  text.  It 
will  not  be  repeated  here.  □ 
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D.2  Proof  of  Lemma  24 

Let  (•,  •)  denote  an  inner  product  on  a  linear  space  X .  For  fixed  x  6  X,  let  f(y)  = 
(x,y).  Then  /  is  continuous. 

Proof:  Let  t/o  €  X  and  e  >  0  be  given.  It  is  sufficient  to  identify  6  >  0  such  that 
||y  ~  3/o||  <  8  =>  |  (*,  y)  -  (x,  Vo)  |<  €•  If  x  =  0  the  result  is  trivial.  Assume  x  0.  Then 

I(*»i0-(s»ifo)|  =  |(*,V-|fo)| 

^  11*11  lb-  !to||* 

(See  [54:pg  273],  [63:pg  41],  or  [66:pg  25l]  for  example,  for  the  Schwarz  inequality.)  But 
||a:||  is  known  since  x  is  fixed,  so  choose  6  <  e/||x||.  □ 

D.3  Comments  on  the  proof  of  Lemma  42 

The  cited  proof  omits  some  of  the  detail  concerning  the  boundedness  of  ||^Q||.  The 
Uniform  Boundedness  Principle  must  be  applied.  See  the  proof  of  Theorem  46  (below)  for 
an  example  of  such  an  application. 

D.4  Comments  on  the  proof  of  Theorem  46 

Theorem:  Let  {A(f)}te[0|r]  be  a  stable  family  of  infinitesimal  generators  of  Co  semi¬ 
groups  on  AT.  If  j D(A(t))  =  D  is  independent  of  t  and  for  vgj),  A(t)v  is  continuously 
differentiable  in  X  then  there  exists  a  unique  evolution  system  U(t,s),  0  <  s  <  t  <  T, 
satisfying  the  5  results  of  Theorem  44  where  Y  is  the  set  D  equipped  with  the  norm 

Iblb  =  Mx  +  m(o)v|U- 

Proof:  See  [58:pp  145-146].  This  appendix  is  to  expand  on  one  portion  of  the  proof 

cited. 

The  proof  includes  a  claim  that  Q(t)  =  Ao I  -  A(t)  is  an  isomorphism  of  Y  onto  X. 
It  is  appropriate  to  comment  on  the  validation  that  Q(t )  is  in  fact  such  an  isomorphism. 
For  Aq  >  w,  Theorem  1.5.3  of  [58:pg  20]  gives  existence  of  an  inverse  for  Q(t)  which  means 
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it  must  be  1-1.  Theorem  2.16.3  of  [22:pg  55]  gives  onto.  (See  also  [19:pg  13].)  Linearity  is 
clear  and  an  algebraic  isomorphism  is  established. 

However,  for  a  topological  isomorphism,  Q  and  Q~x  must  each  be  bounded.  This 
can  be  a  difficult  issue.  First,  consider  the  boundedness  of  Q. 

For  any  v  6  F,  Q(t)v  is  continuous  from  the  hypothesis  on  .4(f).  Hence,  there  is 
some  Mv  such  that  [|Q(i)u||y_.x  <  Mv  for  all  t.  Now  the  Uniform  Boundedness  Principle 
(eg  [65:pg  196])  applies  to  give  a  uniform  bound,  say  Mq  for  ||Q||. 

Next,  recall  that  bounded,  linear,  one-to-one,  and  onto  operators  have  bounded  in¬ 
verses  (eg  [73:pg  70]  or  [22:pg  47]).  Thus,  ||<2-1||  is  bounded. 

Now  it  is  clear  that  the  proposed  isomorphism  is  legitimate. 

It  is  interesting  in  the  above  argument  that  ||Q||y_+x  is  taken  to  be  bounded  by 
hypothesis,  yet  the  norm  is  specified  separately.  It  is  appropriate  to  verify  that  these  are 
consistent,  ie,  that  Q  really  is  bounded  when  considered  as  a  mapping  from  the  given  Y 
to  the  given  X. 


||Q(0lly-*  = 


< 


sup 

o^yey 


llQ(*)y|U 

l|y||y 


snn  U(i-M(*))yi|x 

o*ylv  IM|y 

ll(I-M(0)-/l(0)-M(*))ylU 
0  3?y  ||y||y 


^  ll(J+i4(0))y|U  +  ||(A(Q-i4(0))yHx 
o£Iy  ■  l|y||y 


1  +  sup 

OitySY 


H(A(Q- A(0))y|U 
l|y||y 


(200) 


But  the  continuous  differentiability  of  A(t ),  on  a  closed  and  bounded  interval,  gives  Lips- 
chitz  continuity  to  A.  That  is,  for  some  K,  \\A(t)y  -  -A(0)y||x  <  Then 


||Q(t)||y_x  <  1  +  KT. 


(201) 
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D.5  Proof  of  Lemma  61 

Lemma  61  Let  U  and  V  be  evolution  operators.  Then 

V(t,  r)y  -  U(t,  r)y  =  jf  U(t ,  s)(Au(s)  -  Aw(s))F(s,  r)y  ds .  (202) 


Proof:  The  proof  consists  of  taking  the  derivative  of  a  cleverly  chosen  quantity  and  then 
integrating  the  result.  Recall 


ftV(t,s)  =  -A(t)U(t,s) 


(203) 


Now  the  differentiation  is  presented. 

IpMvMi]  =  lim 

ds 1  v  1  v  ’  ’  1  As-fO  As 


,.  U(t,  s  +  As)V (s  +  As,  r)y  -U(i,s  +  As)F (s,  r)y 

=  lim  - - - 

As-*0  A  s 


,  U(t,s  +  As)y (s,  r)y  -  U (f,  s)F (s,  r)y 

+  lim  - - - 

As-»o  As 


v  tt/  a  Sv(s  +  As,r)-V(s,r) 

lim  U(t,s  +  As) - t“ - -y 

As-tO  '  As 


As— >0  As 


=  lf(t,5)^y(s,r)?y+ ^(t,5)y(s,r)y 


=  U  (t,  s)(— Av(s)V’(s,  r)y)  +  U(t,s)Au(s)V(s,r)y 


=  ~U(i,s)Av(s)V($,r)y+  U(t,s)Au(s)V(s,r)y 
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(204) 


Now,  integrate  both  sides  from  r  to  t. 


U(t,s)V(t,s)y\3)Zir  =  [ tU(t,s)(Au(s)-A'’(s))V(s,r)yds 

Jr 

U(t,t)V(t,r)y»U(t,r)V(r,r)y  =  jf  U(t,s)(Au(s)  -  Av(s))V(s,r)y  ds 

V(t»r)y-U(t,r)y  =  jf  U(t,s)(Au(s)  -  Av(s))V(s,r)y  ds  □  (205) 

A  similar  argument  is  available  in  [53:pg  552]. 


D.6  Proof  of  Lemma  63 

Lemma:  For  positive  real  numbers  a,  6,  and  c  with  b  >  c,  it  follows  that  <  |. 

Proof:  Let  ri  =  ^  and  r2  =  Notice  that  rj  =  |  +  1  and  r2  =  |  +  1,  so  b  >  c 
implies  r2  >  r2  and  ^  >  1.  Then,  r\b  =  a +  b  and  r2c  =  a  +  c.  Hence 

a  +  b  _  rib 
a  +  c  r2c 

b 

~  ri  r2  c 

=  □  (206) 
c 

D.  7  Proofs  of  Lemmas  70  -  72 

Lemma:  For  a,  b  >  0,  (a2  4-  62)^2  <  a  +  b. 

Proof:  For  positive  values  of  a  and  b  it  is  clear  that 

(o  +  b)2  =  a 2  +  2  ab  + 

>  a2  +  £>2.  (207) 

The  desired  result  follows  by  taking  square  roots.  □ 
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Lemma:  For  any  real  numbers  a  and  6,  (a  +  h  <  2(a2  +  62). 
Proof:  Note  that 


Now, 


(a -by 

> 

0 

a2  -  2a6  +  b 2 

> 

0 

a2  +  62 

> 

2ab. 

(208) 

(a  +  6)2  =  a2  +  2a6  +  b 2 

<  2(a2  + 12).  □  (209) 


Lemma:  For  positive  real  numbers  a,  b,  and  c  it  holds  that  (a+6+c)2  <  3(a2+62+c2). 
Proof:  This  is  a  straightforward  calculation. 

(a  +  b  +  c)2  =  a2  +  2  ab  +  2oc  +  b2  +  2  be  +  c2 

<  a2  +  (a2  +  &2)  +  (a2  +  c2)  +  62  +  (62  +  c2)  +  c2 
=  3(a2  +  62  +  c2)  □  (210) 

Lemma:  For  any  positive,  real  a,  b,  and  c,  it  holds  that  (a  +  fc  +  c)1/2  <  (a1/2  +  61/2  + 

c1/2). 

Proof:  This  is  the  same  as 

a  +  6  +  c  <  (a1/2  +  61/2  +  c1/2)2.  (211) 


(ai/2  +  6i/2  +  ci/2)2  =  a  +  6  +  c  +  2a1/261/2  +  2a1/2c1/2  +  261/V/2 


and  the  desired  result  follows  by  taking  square  roots. 


□ 
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or 


ii(i+ AA)-»*n<  mi. 


(216) 


But  || (I  +  AA)_1||  =  sup||.||_1{||(J  +  AA)-1;?||}  and  the  original  form  is  obtained  from  the 
last  equation  above.  For  the  nonlinear  case  the  expanded  version  is  used  for  the  notion  of 
accretive. 

It  is  also  useful,  in  the  nonlinear  case,  to  adopt  a  set  theoretic  representation  of  the 
operator  A.  That  is,  for  A  C  X  x  X,  define 

Ax  =  {y  :  [z,y]  g  A} 

D(A)  =  {a  :  Ax  ^  <f>} 

R(A)  =  U{Ax:x€D(A)}  (217) 

and  note  that  this  allows  A  to  be  a  multivalued  function. 

For  A,  B  C  X  x  X  and  A  g  make  the  natural  definitions 

A  +  jE?  =  {[x,y+ z]:y  e  Ax,  z  e  Bx} 

A  A  =  {[a:,  Ay] :  y  g  Ax} 

A'1  =  {[2/>®]:[s,y]e  A).  (218) 

In  this  setting,  call  B  C  X  X  X  accretive  if  (/  +  A B)~l  is  a  function  for  all  A  >  0 
and  ||(J  +  \B)~lx  -  (I  +  Af?)-1y||  <  ||a;  -  y||  for  all  x,y  g  D  ((/  +  A B)~l). 

Theorem  80  (Crandall-Liggett)  Let  X  be  a  Banach  space.  Let  A  C  X  x  X  and  u>  g 
satisfying  A  +  wl  is  accretive,  be  given.  If  R(I  +  A  A)  D  D(A)  for  all  sufficiently  small 
positive  A,  then 

lim  (/+  —A)~nx  (219) 

n-+  oo N  n 

exists  for  x  g  D(A)  and  t  >  0.  Moreover,  if  S{t)x  is  defined  by  this  limit,  then  S(t )•  g 
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Note  that  A  is  not  required  to  be  closed  or  densely  defined  or  linear  or  accretive.  Also, 
note  that  u>  and  A  are  now  fixed. 


The  proof  is  best  approached  by  first  establishing  a  host  of  preliminary  results.  This 
will,  of  course,  require  more  notation.  Let 

J\  =  (I+XA)-1 
Dx  =  D(J\) 

\Ax\  =  inf  {||y||}.  (220) 

yGAx 

Lemma  1.2.1:  Choose  A  >  0  such  that  A cj  <  1.  Then  J\  is  a  function,  and  for  x,y  € 

\\J\x  -  Ja»||  <  (1  -  Aw)_1||a:  -  y\\.  (221) 

Proof:  The  proof  of  this  lemma  begins  with  a  claim. 

Claim:  (/  +  pA^A)-1  has  Lipschitz  constant  |1  +  tu\  for  t  >  0,  1  +  tu  ^  0. 

This  is  shown  in  the  following. 

(222) 

=  (1  +1^)  [((1  +  tv)I  +  M)_1  *  -  ((1  +  tu)I  +  M)"1  y 

=  |1  +  tu>\  ||  (J  +  t(ul  +  A))-1  x  -  ( I  +  t(ul  +  A))-1  r/|| 

<  jl  +  tw|  ||®  -  y\\  (223) 

where  the  ineouality  follows  by  the  accretiveness  of  A+uC  for  t  >  0,  and  directly  for  t  =  0, 
and  establishes  the  claim. 

Now,  choose  t  =  then  A  =  where  t  >  0.  Hence,  only  those  t  such 
that  tu  >  -1,  need  be  considered.  The  claim  gives  |1  +  M  as  a  Lipschitz  constant  for 
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(I  +  AA)-1,  i.e.  for  J\.  But 


|l  +  to>|  =  1  +  tw 

i 

A 

=  (1  -  Aw)-1  (224) 

and  the  Lipschitz  constant  has  been  validated. 

Next,  establish  that  J\  is  a  function.  The  key  is  that  A  +  u>I  is  accretive  and  hence 
(7  +  t(A  +  w/))-1  is  a  function  for  all  t  >  0.  In  particular, 

(7  +  t(A  +  w7))  -1  =  (7  +  tA  +  tijjl)  1 
=  ((l  +  tu))I  +  tA)  1 

=  (1+,w)  (l  +  fw)~1((l+iw)/  +  M)"1 
=  +  (225) 

For  1  +  tu>  0  (legitimate  as  above)  it  follows  that  (7  +  jq—yl)-1  is  a  function  for  all  t  >  0. 
Choose  t  such  that  A  =  as  above  and  J\  is  a  function.  This  completes  the  lemma.  □ 

As  a  corollary  to  Lemma  1.2.i 

IW*  -  <  (1  -  Awm« y\\  (226) 

for  all  positive  integers  n. 

This  is  established  by  induction.  The  lemma  has  already  established  the  base  case. 
Suppose  now  that  ||J”-1*  -  J”-1i/||  <  (1  —  Aw)“(n-1)||*  -  y\\.  Let  *  =  J"-1*,  y  = 

Then 

=  px*-m 

<  (1  -  Aw)-1  II*  -  2/|| 

<  (1  -  Aw)“n||*  -  2/||  (227) 
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and  the  corollary  is  established. 

Note  that,  in  the  above  argument,  it  is  assumed  that  J^x  G  D ^  for  k  =  1,2,  •  •  • ,  n—  1. 
Also,  the  fact  that  J\  is  a  function  was  used  to  show  tb-t  J\x  is  uniquely  determined.  This 
can  be  taken  as  an  additional  hypothesis  without  preventing  its  desired  use  later. 

As  a  specific  application  of  the  corollary,  choose  x  =  J\y.  Then 

K+1S  -  JiJ/ll  <  (1  -  Ao.)-”||JaS,  -  A  (228) 

Lemma  1.2.ii:  Choose  A  >  0  such  that  Aw  <  1.  Then 

||J^a:  -  .r|j  <  A(1  -  Aw)_1|A*|  (229) 


for  x  G  D ^  n  D(A). 

Proof:  Let  [x,j/]  G  A  be  given  with  x  G  D\  H  D{A).  Let  Xi  =  (I  +  AA)_1a:.  Then 

(I  +  AA)xi  =  x 

x  =  xi  +  AAxi 

=  si  +  Ayx.  (230) 

Recall  that  J\x  =  is  an  alternate  notation  for  (/  +  AA)-1x  =  x\.  Now 

||Jax-x|!  =  ||a:i-x|| 

=  H*/a(»i  +  Ayj)  —  Ja(«  +  Aj/)|| 

<  (1- Aw)-1||(*1  +  Aj/1)-(x  +  Ay)|| 

=  A(1  -  Aw)-1||y||  (231) 


which  is  valid  for  all  y*  HcncC) 


\\J\X  -  x]| 


<  A(l  -  Aw)  1  inf  {\\y\\} 

y£Ax 

—  A(1  -  Aw)-1|Ax| 


(232) 
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as  desired.  Note,  in  the  above  arguments,  that  (I  +  AA)x  =  x  +  XAx  =  x  +  Xy  means  that 
x  =  J\{x  +  Xy)  for  all  x  £  D(A).  Also,  x\  +  Xyi  -  x  =  0.  □ 

Lemma  1.2.iii:  Choose  A  >  0  such  that  Aw  <  1.  Let  n  be  a  positive  integer, 
x  £  and  A|w|  <  1.  Then, 

||  -  x||  <  n(l  -  A[w|)-n+1||Jx*  -  ®||.  (233) 

Proof:  Observe  that 


Jxx  -  x  =  J%x  -  J*~lx  +  Jx  rx  —  •  -  •  +  J\x  -  x 

=  £  (•nr**  -  jr(<+1>)  <234) 

3= 0 

Claim:  x  £  D(J\)  =>  x  £  j D(J"-1). 

Since  x  £  D(J%),  (7+AA)“nx  is  well  defined.  But,  (J+AA)-nx  means  [(/  + AA)-1x]n 
and  its  well  definedness  implies  that  of  (J  +  XA)~*x.  Finally, 

(I  +  AA)"nx  =  (I+XAy^il+XA)-^ 

=  (I+\A)-(n~Vjxx 

=  (235) 


and  the  well  definedness  of  all  the  other  terms  gives  that  of  J"-1 .  Inductively,  the  claim 
holds  for  all  terms  of  interest.  So,  all  of  the  newly  introduced  terms  make  sense.  Thus, 

k*-*i  = 

t=0 

<  £(1  -  Aw)-n+(^)||  JAx  -  *||  (236) 

t=0 

by  Lemma  1.2.i.a. 

There  are  three  cases  to  consider.  Case  1:  Aw  =  0.  This  case  is  trivial. 
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Case  2:  Aw  <  0.  It  is  desired  to  show  that 


n-l 

2(1  “  Aw)-n+(*'+1)  <  n(  1  -  A|w|)“n+1.  (237) 

i=0 

But,  Aw  <  0  =t>  0  <  -Aw,  1  <  1  -  Aw.  Hence,  the  terms  in  the  summation  are  increasing. 
The  largest  term  is  (1  —  Aw)-n+(n-1+1)  =  1  and  the  sum  is  less  than  or  equal  to  n.  This 
case  will  be  complete  if  (1  —  A|w|)  <  1  (since  then  the  inverse  will  be  greater  than  or  equal 
to  1)  which  the. hypothesis  A|w|  <  1  clearly  assures. 

Case  3: .  Aw  >  0.  This  time 


-Aw  <  0 

1  -  Aw  <  1  (238) 


Hence,  the  largest  term  in  the  sum  is  the  first  one,  (1  -  Aw)-n+1.  The  sum  is  less  than 
or  equal  to  n(l  -  Aw)~n+1.  The  proof  will  be  complete  if  (1  -  Aw)-n+1  <  (1  -  A|w|)-n+1. 
But  Aw  >  0  =»  w  >  0,  i.e.  |w|  =  w  and  the  equality  holds.  □ 

Lemma  1.2.iv:  Choose  A  >  0  such  that  Aw  <  1.  If  A  >  0,  n  6  9i,  and  x  G  D^,  then 

jx  +  Xx  G  (239) 

and 

Ja*€  +  (240) 

Proof:  Since  x  €  D\,  J\x  makes  sense.  Say 


J\x  —  (J  +  AA)-1a: 

=  ®0,  (241) 

whence 

x  =  (I  +  AA)xo 
=  xq  +  Ai4o:o 
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=  s0  +  Xyo 


(242) 


where  [20,2/0]  £  A.  Now  write 

a  \  —  u  -  a,  .  .  A  -  u 

-s  H - ^ — J\x  =  -(®0  +  Aj/o)  4 — ^ — ®o 

u  a 

=  -20  +  MO  +  So  -  ^®0 

=  So  +  Mo  (243) 

It  remains  to  show  that  so  +  Mo  G  But,  [so,yo]  G  A ,  so  so  €  D(A)  and 
(7  +  nA)xo  =  so  +  Mo-  This  means  (7  +  /x.4)-1(so  +  Mo)  =  *0  and  the  first  part  of  the 
lemma  follows.  But,  recall  J\x  —  20  and  so  G  7M(so  +  Mo)-  Hence,  J\x  G  JM(s 0  +  Mo)  as 
required.  □ 

As  a  corollary:  If  s  €  D(J%),  then  J^z  G  (^7*-1s  +  for  all  integers  k. 

The  lemma  establishes  the  base  case.  It  remains  to  show 

4*  <=  4  +  ^4*)  *  4"*  6  J„  (2jj[*  +  (244) 

Let  s  =  J*s.  Then,  s  G  7^  ^7*-1s  +  ^j^s)  and  the  desired  result  is 

7A*G7^s  +  ^7As).  (245) 

Since  s  G  D(7A),  let  J\x  =  (7  +  AA)-1s  =  sx.  Then  s  =  (7  +  AA)si  =  si  +  AAsi  = 
*1  +  A y\.  Note  that  [sj,yi]  G  A.  Write 

M-  .  X-n  T  .  H,  lN  N  ,  A- n 
jx  +  -y-J\x  =  ^(sj  +  Ayi)+ — si 

=  Si  +  Ml  (246) 


Consider 


(7  +  nA)  si  =  Si  +  Mi 

si  =  (7  +  /iA)_1(si  +  wi) 


(247) 
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i.e.  (®i  +  Mi)  £  Dp.  But,  J\x  =  a:i  and  ®i  €  1  +  Mi)  so  Ja»  €  JV(»i  +  Ml)  = 

(jr*  + 

Before  presenting  the  next  lemma,  a  brief  review  of  some  identities  among  the  bino¬ 
mial  coefficients  is  appropriate. 


B(n,  m) 
B(n,n) 
B(n,0) 

B(n,  m)  +  B(n,  m  +  1) 


n! 

m!(n  -  m)\ 

=  1 

=  1 

=  B(n  +  l,m+l) 


(248) 


Lemma  1.3:  Let  A  >  fi  >  0,  Aw  <  1,  and  x  6  D(J ™)  n  D(J£),  where  m  and  n  are 
positive  integers  with  n>m.  Let  a  =  (3  =  Then 


m-l 

<  (i-u»rnY:*iPn-jB(n,j)\\jrj*-*\\ 

3=0 

+  (1  -  w n)-’ampj-mB{j  -  l,m  -  l)!!^*  -  x\\.  (249) 

j=m 

An  abbreviation  will  be  useful.  For  integers  j  and  k  satisfying  0  <  j  <  n  and  0  <  k  <  m, 
let 

ak,j  =  \\4x-Jkxx\\.  (250) 

Also  introduce  ai  =  (1  -  //w)-1  §  and  Pi  ■■=  (1  -  /iw)-1  In  this  notation  the  result  of 
the  lemma  takes  the  form 


m-l  n 

am,n  <  al3Pin~3S(nJ)am-j,0  +  JZ  «1  ™ P\3~m B{j  -  1,  TO  -  l)a0,n-j.  (251) 

i=0  j=m 

Proof:  The  proof  begins  with  a  basic  inequality  relating  the  ak,r  This  will  require  several 
steps  to  establish.  In  the  definition  of  a^j  it  is  desirable  to  make  the  replacement 

4x  =  J„  fat'*  +  ^4*)  •  (252) 
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This  follows  from  the  corollary  to  Lemma  1.2.iv  and  Lemma  1.2.i.  Now, 

ak,j  —  ||  —  ^a^II 

=  |4*-^(f4~'*+^p4*)| 

<  (i-/*o‘1||4'I*-(f4"1*+^4*)| 

<  (l  -  -  4*11 

=  aiak-i,j-i  +  Pidkj-i  (253) 

where  the  first  inequality  follows  from  Lemma  1.2.i,  the  second  inequality  uses  ^  +  =  1 

and  the  triangle  inequality,  and  the  last  step  is  simply  a  change  of  notation.  This  is  a  basic 
inequality  to  be  used  in  the  proof. 

The  proof  is  by  induction  on  n.  The  inductive  proposition  P(n )  is:  For  all  m  <  n 

m- 1  n 

am,n  <  Y  <*\}  P*'3  B(n,j)am-jfi  +  Y  axmPi3~mB{j  -  l,m  -  l)a0,n-j.  (254) 

j—0  j—m 

Let  n  =  1,  for  the  base  case.  Note  that  m  =  1  is  the  only  possible  value  for  m.  Claim: 

The  following  inequality  holds. 

o  1 

oi,i  <  £«iJ/5i1--'5(l,j)a1_i)0  +  2ai1/3i:'”lj5(i- 1>°)ao,i-i 

i= o  i=i 

=  ftoi.o  +  ao,o-  (255) 

The  basic  inequality  says  a i(i  <  oiao.o  +  Pi^ijo  and  the  base  case  is  complete.  Note 
that  ao,o  =  0. 

Now  suppose  that 

to-1  n 

am,n  <  Y  a13p1n~3B(n,  j)am-j<0  +  Y  <*\mPi3~mB{j  -  l,m-  l)o0,„_j  (256) 

i=0  j=m 
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for  all  m  <  n  and  seek  to  establish 


m-l  n+1 

flm.n+1  <  E  c*i3 Pin+1~3 +  +  E  aim  Pij~m  B(j  -  l,ro- l)a0,„+i-j  (257) 

i=0  1  j=m 

for  all  m  <  n  +  1.  Start  with  the  basic  inequality, 


®m,n+l  —  &lGm— l,n  d*  Pi ®m,n- 


(258) 


For  each  m  <n  the  induction  hypothesis  applies  to  yield 

(m— 2 

E  Gl3 Pin~3 j)am-l-jfl 

J=o 

+  E  aim-1/V-(m-1)j^(j  -  l,m  -  2)a0|„-j 

i=m-l 

(m-l  n 

E  aij/?in-J'5(n,  j)am_i>0  +  E  «i mPij~mB(j  -  1,  m  -  l)aot„-j 

j-0  j-m 

m-2  n 

=  E  «ij+Vin"i5(n,i)am_(i+a)>0+  E  aim/5iJ+1_m5(y-l,m-2)a0ln-j 

i=0  j=m-l 

m-l  n 

+  E  «ijAn-i+15(n,  j)am_jl0  +  E  <*\m P\**l~m B{j  -  1,  m  -  l)a0,n_(259) 

j=0  j'=m 

Note  that  m  =  n  + 1  and  m  =  1  must  be  treated  as  separate  cases.  This  will  be  done 
later.  For  now,  focus  on  the  first  and  third  summations  in  the  preceeding  inequality. 


m— 1 


m-l 


E  oti3 Pin  3+1  B(n,j  -  l)am_j,0  +  E  aiiPin~3+1B{n,j)am-j,o  +  ^n+1am, o 

i=i  j=i 

m-l 

=  E  <*i3Pin~3+1  [B(n,j  -  1)  +  B(n,j)\am-j,o  +  /?in+1am>0 
i=i 

m-l 

=  E  ai 3Pin~}+1B(n  +  l,j>m-j,o  +  Pin+1am,o 
i=i 


m-l 

=  Ea^i  n"i+15(n  +  1,j>m_ji0 
j=o 


(260) 
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Next  consider  the  second  and  fourth  summations. 

n 

E  ai”W+1_m  [B(j  -  1,  m  -  2)  +  B(j  -  1,  m  -  1)]  a0,n-j  +  aima0,„_(m_i) 

j=m 

n 

=  E  ai”W +1~mB(j,  m  -  1  )a0>n-j  +  aima0,n-(m-i) 

j=m 

=  E  «1W+1“m5(i>m-l)a0l„_j 

j=m-l 

n*fl 

=  E  <*™(hj-mB(j  -  1,  m  -  l)a0l»+i-i  (261) 

j=m 

It  is  easy  to  see,  when  the  terms  are  recombined,  that  the  desired  result  is  obtained. 
Return  now  to  the  case  m  =  i.  It  is  necessary  to  show 
0  n+1 

«i,n+i  <  E  °» jPin~j+1B(n  +  +  E  -  1, 0)ao,„+i_j 

i~o  i=i 

n+l 

=  /?in+1fli.o  +  E  Q?i/?iJ~1flQ.w+i-j  (262) 

i=i 

This  must  be  established  for  all  n  so  another  induction  is  in  order.  For  n  =  1,  it  is  required 
to  establish 

2 

al,2  5:  /?l2fll,0  +  E  qi/V~1q0,2-7 

i=i 

=  /5i2ctiio  +  aiao.i  +  ai/?iao,o>  (263) 

But,  from  the  basic  inequality 


al,2  <  »la0,l  +  PlO-1,1 

=  o;iao,i  +  /?i(QJao,o  +  /?i«i,o)  (264) 


as  required.  (Note  again  that  ao,o  =  0.) 
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Claim:  The  following  holds. 


n+2 

ai,n+2  <  /?in+2ai,o  +  ^2  ai/V-1ao,n+2-j  (265) 

3=1 

Prom  the  basic  inequality  and  the  induction  hypothesis  it  follows  that 


®l,n+2 


.<  CiiOo.n+l +/?lfll,n+l 

^  <*iao,n+i  +Pi  ^in+1ai,o  +  ^  «i/^iJ_1ao,n+i-i j 

n+1 

=  «iao,n+l  +  /?in+2ai,0  +  ^2  al Pi*  aO,n+l-j 

3=1 

n+l 

=  Pin+2aifl  +  ^2  &Otn+l-j 

3=0 
n+1 

=  /5in+2ai,0  +  ]> 2  alPl*~1<l0,n+2-j 

3=1 


(266) 


as  required. 

Now  consider  the  case  m  =  n  +  1.  This  time  it  is  required  to  show 


«n+l,n+l  <  *  Bin  +  ^JW+l-j.O 

3=0 

n+1 

+  J2  ain+1AJ-(n+1)5(y  -  l,n)flo,„+i-i 

j=n+l 

=  Y2  <*i*  Pin+1~*  B(n  +  1>  i)«n+i-i,o  +  «in+1a0)o>  (267) 

3=0 

This  can  best  be  handled  by  establishing  the  following  intermediate  result.  Let  all  con¬ 
ditions  be  as  stated  in  the  Lemma,  except,  consider  n  <  m.  The  result  then  becomes 
am,n  <  Dj=0  (*i* Pin~*  B(n,j)am-jfl-  Notice  that  once  this  is  established,  the  desired  re¬ 
sult  follows  by  replacing  m  and  n  with  n+1  and  noting  that  the  n  +  1  term  is  zero. 
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With  attention  on  m  >  n  it  is  appropriate  to  represent  m  by  n  +  k  where  k  is  any 
appropriate  nonnegative  integer.  Then,  proceed  with  induction  on  n  for  the  proposition 

IX 

P(n)  =  For  every  k,  an+k,n  <  ')T,a13Pin-3B(n,j)an+k-j,o .  (268) 

i=o 

Let  n  =  i.  Then 

ai+Jt,i  ^  (269) 

by  the  basic  inequality.  This  clearly  has  the  desired  form  and  the  base  case  is  complete. 
Suppose  P(n )  in  order  to  establish  P(n  +  1). 


fln+l+fc.n+l  <  ai''!n+fc,n  +  /?lan+l+fc,n 

n  n 

<  <*i  £  ai3/3in~JB(n>i)an+k-jto  +  A  £  aiJ/?in"JP(n,i)a„+it+i-i,o 
j= o  i=o 

=  £ai,+Vin_J5(n,  j)a„+fc-;,o  +  £  ocij  Pin+1~j  B(nJ)an+k+1-jfi 

3- 0  3=0 

n+1  _  n 

=  £  ®ij0in+1~jB(n,j  -  1  )an+k+i-jfl  +  £  ai: An+1-J B(n,j)an+k+i-jfi 
j=i  j=o 

n 

=  ain+1a*,o  +  £ai,0in+1*vli?(»,  j  ~  1)  +  B{n,j)]an+ jt+W.o 

3=1 

+Pin+1  an+k+i,o 

n 

=  ocin+1ak,o  +  £  <*i^ Pin+1~i B(n  +  l,i)an+Ar+i-i,o  +  /?in+1«n+Jt+i,o 
3=1 

n+l 

=  £  aiJ'/3in+1-J'P(n  +  l,i)an+Hl_ii0  (270) 

3=0 

and  the  induction  is  complete. 

This  completes  the  proof  of  the  lemma.  □ 

Before  the  next  lemma,  it  is  appropriate  to  recall  the  Schwartz  inequality  in  the  form 


EM.l<(ENl2),/J(Ew2) 


,1/2 


(271) 
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where  the  generic  sums  suggest  that  the  result  holds  for  both  finite  and  infinite  sums,  eg 
[54:pg  548]. 

Several  other  preliminaries  are  also  in  order.  The  quantity  (a  +  /?)"  will  appear. 
Wiite  this  in  its  binomial  expansion  and  take  derivatives  with  respect  to  a  to  obtain  some 
useful  identities. 


E5(«.iK'rj 

3=0 

3=0 

j2jB(nJ)a?pn~j 

3=0 

3=0 

3=0 


(a  +  p)n 
n(a  +  /?)n_1 
an(a  +  /?)n-1 


an(n  -  l)(a  +  P)n~2  +  n(a  +  P)n  1 


a2n(n  -  l)(a  +  P)n  2  +  an(a  + /3)n  1 


(272) 


Recall  now,  that  a  +  P  =  1.  The  above  then  simplify  to 

E5(«,iVri  =  i 

3=0 

n 

]> ZjB(n,j)oPpn~j  =  an 

3=0 

n 

^y'ij2B(n,j)o?prx~i  =  a2n(n-l)  +  cm  (273) 

3=0 


Some  similar  identities  are  obtained  from  the  McLaurin  series  for  (1  -  p)  m.  The 
series  looks  like 

P2 

(1-P)~m  =  1  +  mP  +  m(m+  l)y  +  •••.  (274) 

But 


U  ~  1)! 

(m  -  l)!(j  -  1  -  (m  -  1))! 
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(i  -  i)i 


(m  -  l)!(j  -  m)! 


(275) 


and  it  is  clear  that 


E  -  1,  m  -  l)F~m  =  (1  -  p)~m  (276) 

j=m 

The  next  step  is  to  take  derivatives  with  respect  to  /3  on  both  sides.  This  yields 

OO 

j=m 

oo 

=  m/3(l  -  (3)~m~l 

j=m 

oo 

E  O'  -  m)2B(j  -  1,  m  -  l)/33-m_1  =  m(m  +  1)/3(1  -  /3)-m"2  +  m(l  -  /S)-”1"1 

j=m 
oo 

^(i-m)25(i-l,m-l)j83_m  =  m(m  +  1)/32(1  -  /3)" 


\-m-2 


j=m 


+m/3(l  -  /?)’ 


-m-1 


(277) 


Recall  that  1  -  (3  =  a  and  reduce  these  to 

E  B(j  -l,m-  l)am(3i 

j=m 

oo 

E  O’  “  m)BU  -  1> m  -  l)«m/?J 

j=m 

oo 

E  O’  -  m)2£0  -  1,  W  -  l)am/3; 


— m  _ 


=  1 


-m  _ 


—m  _ 


m/3 
a 

m(m  + l)/32  m/3 


0—m 


of 


+  — 
a 


(278) 


Some  additional  modifications  to  the  form  of  these  equations  is  desired.  To  prepare  for 
them,  note  that 


<*  +  /?  =  1 

ml3a  +  m/32  =  m/3 
m/3  m/32  _  m/3 

a  a2  a2 
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(279) 


m/32  _  to/?  to/? 

a2  a2  a 


Now  write 


—m  _ 


j=m 

oo 


TO/? 

a 


+  TO 


2  -  l.nr  -  l)am/?J 

j=m 


— m  __ 


to(to  + 1)/?2  to/? 


cr 


H - 1-  2to 

a 


(*+»)- 


TO 


_  to(to  +  l)/?2  to/?(1  +  2 to)  ^  ^2 


a' 


a 


+  to 


(280) 


It  is  finally  time  to  state  and  prove  the  lemma  for  which  all  these  preparations  have 
been  made. 

Lemma  1.4:  Let  n  >  to  >  0  be  integers  and  a,  /?  be  positive  numbers  such  that 
«  +  /?=!.  Then 


m 

53  B(nJ)qi(3n~*(m  -  j )  <  ((na  -  to)2  +  na/?) 
i- o 


1/2 


(281) 


and 


B(j  -  l,m  -  l)am/?J-m(n  -y)  <  (— (—  +  to  -  n) 
i=m  V a  '  a  ' 


2\l/2 


(282) 


Proof:  Consider  the  first  inequality.  With  n  >  to,  apply  the  Schwarz  inequality  to  write 


YjB{n,j)oPpn  3{m-j )  <  j|TO-j| 

j=o  i= o 

(  \  1/2  /  \  1/2 

£>(»,  j)A 

<  £*(»,  j)o>7?"-* )  53P(n,y)aJ'/?n-jTO2 

\i=o  /  \i=o 

i/2 

-  53  B(n,j)oP f3n~j(2mj)  +  '£B(n,j)c^ (3n~>j2  ' 
j=o  j= o 
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=  ( 'm 2  -  2 man  +  a2n(n  -  1)  +  anj  ^ 

=  ^(na  -  m)2  +  an(l  —  ^ 

=  ^(na  —  m)2  -f  anfi'j  ^  (283) 


This  completes  the  proof  of  Lemma  1.4.  □ 

Some  final  preliminaries  will  now  be  presented.  These  are  some  small  details  that 
are  best  established  outside  of  the  main  line  of  argument. 

Claim:  If  fj,  >  0  and  n\u\  <  1  then  1  -  fi\w\  <  1  -  nu. 

This  is  equivalent  to  nu  <  n\u>\  and  is  clear.  Thus  (1  -  /z|w|)-1  >  (l  -  juw)-1. 

Claim:  For  sufficiently  small  positive  A,  (1  -  A|w|)-m  >  (1  —  A|w|)-(m-J)  as  long  as 
j  <  m. 

For  small  positive  values  of  A,  it  is  clear  that  1  -  A|w|  <  1  and  hence  (1  -  A|u>|)_1  >  1. 
Then  the  result  is  clear. 
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Claim:  Given  e  >  0,  there  is  some  positive  integer  N  such  that 

— —  —  <  e  whenever  n.m>  N  (285) 

m  n 

Consider 

1__  1 
m  n 

< 

< 

Choose  N  >  \  and  the  result  is  clear. 

Claim:  If  n  is  a  positive  integer  and  t 
The  proof  is  by  induction.  Let  n  =  1. 

(1  -  t)-1  <  e2t  (287) 

Notice  that  for  t  =  0  the  equality  holds.  Also  note  that  for  t  =  |  the  inequality  holds 

strictly.  This  is  also  true  in  the  form  1  <  (1  -  t)e2t.  The  base  case  will  be  completed  by 

showing  that  the  right  hand  side  is  strictly  increasing  on  t  G  [0,|). 

Let  f(t)  =  (1  -  t)e2t.  Then 

f\t)  =  2(1  -  t)e2t  -  e2t 
=  (2-2t-  l)e2< 

=  (1-2  t)e2t  (288) 


6  [0,  |],  then  (1  -  t)  n  <  e2nt. 
It  is  required  to  show 


1  1 

—  H — 
m  n 

2 


min{n ,  rri) 
2_ 

N 


(286) 


S\nro  sforivptiva  ic  plwrpve  r>nc»t*Y°  Kaco  rpco  ic  -noY/ 

AAft  WV  VAftW  U  W4  «  t  Cv  V  A  •  W  AW  W.A  «  t  Wl  W  W  WWA  V*  •  VAAW  w  MW  w  vWw  w  AW  AtW  M  w*  vW*  A 


Suppose  (1  - 1)  n  <  e2n<,  then 


(1  -t)"(n+1)  =  (l-t)"1^-/)-" 

<  (l  -t)-Vn< 
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(289) 


<  e2te2nl 

_  e2(n+l)t 


and  the  claim  is  established. 

This  concludes  the  preliminaries  and  the  proof  of  the  theorem  will  now  be  presented. 
The  first  order  of  business  is  to  establish  the  existence  of  the  limit. 

Let  *  €  D(A).  Assume  A  >  n  >  0,  n  >  m,  A|u>|  <  5.  Assume  x  €  D(J ™)  n  D(J£). 
The  plan  is  to  establish  that  Jt”n*  is  a  Cauchy  sequence  for  rational  values  of  /z  and  A, 
then  an  appeal  to  continuity  will  complete  the  problem.  Most  of  the  previous  lemmas  are 
used  in  the  algebra. 


rri — 1 

I \J>  ~  J?* II  <  (1  -  w/*)""  E  j)\\Jrj*  ~  *11 

1= 0 

+  E  (1 "  w/0 -  1,  m  -  l)\\j;~jx  -  *|| 

j—m 

m 

<  a  -  -  *11 

1=0 

+  E  (1  -  mM )- -  l,m  -  1)11  J^X  -  *|| 

j—m 

m 

<  (1  -  mM)-"  E  <*3  P*'3  B{n,  j)(m  -  j)(  1  -  A|w|)_(m"j)+1||  Jxx  -  x\\ 

j=o 

j—m 


< 


A(1  -  ^|)"n(l  -  A|w|)-m  E  aj(3n-jB(n,j)(m  -  j) 

i=o 


+/x(l  -  /z|w|)-n  E  <*mPj~mB(j  -  1,  m  -  1)(»  -  j) 


3—m 


\Ax\ 


Ae^Hg^AH  ^na  _  m^2  _j_  napj 


\  1/2 


+^e 


2n/i|w| 


(  a'  \  a 


2\  1/2' 


+  m  -  n 


\Ax\ 
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_j_e2M n/1  ^m<\2  _  ^  |_42;| 


=  e2H(m.+mA)  (^n{l  _  mXy  +  n/i(A  _  M))1/2 

^g2|o)|nu  q_  (m,X  _  n/i)2^  ^  j  |/lx|  (290) 

This  establishes  a  slightly  tighter  bound  than  the  paper  gives  for  its  equation  1.9.  However, 
no  additional  strength  is  given  to  the  theorem. 

For  any  given  t  and  sufficiently  large  n  and  m,  legitimate  values  of  /i  and  A  are  given 
by  and  ~  respectively.  Some  preparation  is  needed  to  substitute  this  in  to  the  inequality 
above.  The  purpose  for  all  this  is  to  establish  a  result  for  rational  values  and  then  appeal 
to  continuity.  The  following  are  useful  for  substitution. 

nn  =  mA  =  t 
nn  +  mA  =  2 1 
nn  ~ =  0 


=  mA  -  nn 


=  mA(A  -  n)  (291) 
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Substitution  yields 


IWfci*  -  J?^n*ll 


< 

< 


(292) 


Recall  that  u>  is  fixed  (from  the  very  beginning)  and  that  consideration  is  currently 
being  given  to  a  particular  x.  The  limit 


lim  JP/-X  (293) 

n-*oo 

exists  since  the  terms  have  just  been  shown  to  form  a  Cauchy  sequence  in  a  Banach  space. 
From  the  corollary  to  Lemma  1.2.1 

Pt/nx  ~  Jt/nv\\  <  (l  “  II*  "  2/11  (294) 

Recall  that 

lim  (l-~ \  =ewt.  (295) 

n-»oo  \  n  / 

Let 

S(t)x  =  Jihn  J?/nx.  (296) 

Then  S(t)  has  ewt  as  a  Lipschitz  constant.  Thus,  if  5  is  a  semigroup,  then  S  G  Qw. 

The  definition  of  S  for  elements  of  D(A )  has  been  given.  But  S  must  be  defined  on 
all  of  X.  The  extension  will  rely  on  the  denseness  of  D{A)  in  X.  Let  x  £  JJ(A)\D(A)  be 
given.  There  is  a  sequence  {a:,}  of  elements  in  D{A)  which  converges  to  x.  The  issue  to 
be  resolved  is  the  existence  of 

lim  J?/nx.  (297) 

n-*oo  ‘v7* 
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Let  e  >  0  be  given.  Choose  N\  such  that  i  >  N\  implies  Hz;  -  z||  <  jfzr.  Choose  such 
that  for  any  i,j  >  N\  and  n  >  it  follows  that 


\\^t/nxi  ^i/n^ill  5:  2ew  ||z,'  —  Zj|| 


Let  N  =  maz{iVi,iV2}.  Then 


ll^t/n2*  ~  J?/n*ill  < 


and  the  sequence  is  Cauchy  as  desired.  Also, 


<  2eut\\xi~Xj\ 


=  2ea’i||z,-  -  x  +  x  -  Xj\ 


<  2e' 


2eut 


No v;  it  makes  sense  to  define 


S(t) 


x  =  _lim_  lim  . 


TI—+OO  J— >00 


(298) 


(299) 


(300) 


(301) 


The  Lipschitz  continuity  of  S(t)  will  now  be  established.  Let  r  >  t  >  0  be  given.  For 
x  6  D(A),  legitimate  values  of  fi  and  X  are  given  by  £  and  -  respectively.  Choose  n  =  m 
in  equation  1.9  and  write 


||  J?/nx  -  Jt%z||  < 


((i_r)2  +  ,Z_!)1/2e2H(HT) 


+  ^T~~  +  (r  - 1)2) 1/2  e2|u,|< 

Jim  || JnT,nx  -  J?/nx\\  <  [|i  -  r|e2I-K*+-)  +  (r  -  <)e2M‘]  |Az| 

=  (r  -  t)|Az|  (e2|w|(‘+T)  +  e2^) 


|Az| 


(302) 
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But, 


=  \\S(t)x  -  5(i)*||  (303) 

Thus  S(t)x  is  Lipschitz  continuous  in  t  for  x  £  D(A)  and  bounded  t  intervals. 

Claim:  For  a:  G  D(A),  S(t)x  is  continuous  as  a  function  of  t. 

Let  x  6  D(A)\D{A)  be  given.  There  is  a  sequence  {*,-}  ->  x  where  each  X{  is 

in  the  domain  of  A.  Let  e  >  0  be  given.  It  is  required  to  find  some  6  >  0  such  that 

||5(t)x  -  S(r)a;||  <  e  whenever  ||r  -  t\\  <  6.  But,  for  any  x,- 

||S(7>;  -  S(t)xi\\  <  |i  -  r|  \Ax\  (e*H(«+r)  +  e2W‘)  (304) 

Furthermore,  ||S(i£)xt-  -  S(t)x\\  <  e/3  for  sufficiently  large  i.  Then 

||S(r)x  -  S(i)x||  =  ||S(r)x  -  S(r)xt-  +  S(t)x{  -  S(t)xi  +  S(t)x{  -  S(t)x|| 

<  \\S{r)x  -  S{t)x{\\  +  ||S(t)x{  -  S(t)x{ ||  +  || S(t)x{  -  S(*)x|| 

<  e/3  +  \t  -  r|  |Ax|  (e2M(‘+r)  +  +  e/3.  (305) 

All  that  remains  is  to  choose  6  such  that 

1 1  -  t\ \Ax\  (e2M(t+T)  +  e2^)  <  e/3  (306) 

That  is,  it  is  required  that 

S  \Ax\  (e2M2t‘  +  e2H‘)  <  e/3  '  (307) 

where  i  is  the  maximum  allowable  value  for  t. 

Choose 

S  ~-= - t — - r  (308) 

4|Ax|  (e4M«  +  e2M*J 
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and  the  claim  follows. 


It  is  clear,  from  the  definition  of  the  J\  terms,  that  S’(O)  is  the  identity.  The  previous 
claim,  applied  for  t  =  0  gives  the  continuity  required  for  5  G  Qw.  The  proof  is  concluded 
with  the  establishment  of  S(t  +  r)  =  S(t)S(r).  As  has  been  the  pattern,  there  are  a  few 
preliminaries.  Namely,  by  applying  the  definition  and  previously  established  continuity 
results,  it  follows  that 


w  =  fay.%]” 

=  feh/J” 

-  jaw 

(309) 

Also,  it  now  follows  that 

S{mt)  =  Km  J”t/n 
=  jgn  Jmt,mk 

=  faN" 

=  [S(<)]m  (3W) 


It  will  now  be  established  that  the  semigroup  property  holds  for  rational  values  of  t 
and  r.  Let  l,  k,  r,  and  s  be  positive  integers.  Then 


(311) 


But,  since  S(t)  is  continuous  as  a  function  of  t  and  Lipschitz  as  an  operator  on  X,  it 
follows  that  the  relationship  holds  for  all  real  numbers.  □ 
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